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The equations for the energy spectra of electrons and photons in a cascade shower are 
written in a form suitable for numerical applications with accurate cross sections. Trial 
calculations were carried out to check the feasibility of a step-by-step numerical integration 
procedure similar to that used successfully at low energies in noncascade problems. Numerical 
results were obtained for lead by using a monoenergetic source of electrons at 360 mec? and 
extreme relativistic cross-sectional forms, which assume complete screening. The method 
was designed to permit the determination of cascade-shower spectra, spatial distributions, 
and directional distributions that are exact in the sense that the main limiting factor is the 
accuracy with which the cross sections are known. The numerical results of the trial caleula- 
tions were found to agree with previously published results based on more restrictive methods 
of analysis 


1. Introduction 


High-energy electrons and photons dissipate their energy by generating secondary electrons 
and photons. The complex radiation arising in this manner from a high-energy electron or 
photon source is called a cascade shower. Most attempts to describe this phenomenon have 
used approximations that are valid only at very high energies, i. e., in the extreme relativistic 
range; and the results have only limited applicability in a broad range of lower energies where 
the cascade phenomenon is still of importance. At these lower energies realistic calculations 
must be largely numerical because the accurate cross-sectional forms are too complicated to 
permit analytic transformations. 

In this paper a calculational procedure which applies without restrictions on cross-sectional 
forms, and which has been used in noneascade problems, is adapted to the cascade situation. 
It should provide an adequate way of describing cascade-shower spectra in any energy range. 
Because the equations are part of a linked set of moment equations, which describe also spatial 
and directional distributions, and possibly fluctuations, the procedure in this paper may 
eventually form the basis for a more complete description of many aspects of the cascade 
problem. 

The equations are applied, in a trial calculation, to the special case of a monoenergétic 
electron source. Calculations for an electron source are more important than for a photon 
source for the following reason: If a photon source is assumed, it is possible to calculate the 
resulting photon spectrum at all energies. This type of calculation has been done many 
times and offers no difficulty (for an example, see reference {1]).2. One can then determine by 
simple quadrature the electron spectrum generated by this photon spectrum. The whole 
cascade spectrum of electrons and photons, except for the initial photon component already 
determined, can then be determined by appropriate superposition of solutions for mono- 
energetic electron sources. Thus, calculations for a monoenergetic electron source demonstrate 
that the main difficulty blocking a complete numerical analysis of the cascade problems is 
removed. 

For the initial calculation, extreme relativistic cross-sectional forms that assume complete 
screening are used. They include the same types of singularity as the exact cross sections 
and therefore provide a meaningful trial run. Also they are quite similar to the cross-sectional 
forms used in earlier applications of conventional cascade theory; so they make possible a direct 
check with previously published results. 

1 "Phis work © is supported by the Office of Naval Research and the Atomic Energy Commission 


2 Submitted to the University of Maryland in partial fulfillment of the requirements for the M. 8. degree 
Ficures in brackets indicate theli ture references at the end of this paper. 
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2. Interactions Involved in Cascade Showers 


The two tvpes of interactions that predominate in the development of cascade showers 
are bremsstrahlung and pair production. In a bremsstrahlung interaction a high-energy 
electron or positron, under the influence of a nuclear or electronic field, generates a photon. 
In pair production a high-energy photon in the electromagnetic field of a nucleus or electron 
gives up all its energy to generate an electron and a positron. The combination of these two 
processes results in the creation of photons, electrons, and positrons in increasing numbers at 
lower energies. 

At all energies below the extreme relativistic range, there are additional effects that influ- 
ence the development of the showers. There are energy losses due to atomic interactions that 
slow down the electrons more or less continuously along their path. Some of these interactions 
result in secondary electrons with high energies, 1. e., knock-on electrons. Compton scattering 
also contributes some high-energy electrons to the shower, but at low energies it is of primary 
importance through its influence on the photon spectrum. The photoelectric effect and 
photonuclear interactions have cross sections that are much smaller than the other interactions 
at all cascade-shower energies, though the former is of importance below shower energies. 


3. Discussion of Earlier Work 


An excellent discussion of the lines along which the mathematical description of the cascade 
phenomenon has been developed can be found in a book by B. Rossi [2]. Because the proce- 
dural methods of this paper are directed mainly toward a solution of the cascade problem at low 
energies where accurate cross-sectional forms are needed, comments here on earlier work will 
be restricted to developments pertinent to these lower energies. 

The first attempt to describe a low-energy cascade spectrum was made by Arley [3]. He 
neglected interactions that slow down the electrons, i. e., ionization losses, above the critical 
energy;* and he further neglected bremsstrahlung and the generation of secondary electrons 
by all processes at energies below the critical energy. His results as compared with more 
recent calculations are low by a factor of 3 or 4. 

Tamm and Belenky [4] have obtained a much more accurate solution by an analytic 
method, using asymptotic cross-sectional forms for bremsstrahlung and pair production, assum- 
ing continuous ionization loss and neglecting the Compton effect and knock-on electrons. 
Their solution has proved to be more accurate than was at first expected, considering the 
approximations. Friedman [5] has discussed their solution at length and computed a first 
correction term. This correction term never exceeded 4.5 percent of the original solution in 
the case that Friedman considered. 

Perhaps the most satisfactory treatment of the energy-spectra problem for low energies 
has been given by Richards and Nordheim [6]. They have used reasonably accurate repre- 
sentations for all important interactions, but they treat small radiation and ionization losses 
as continuous, thus incurring an error of up to 8 percent, according to their own estimates. 
They obtain results by numerical integration for energies from 4 to 10 times the critical energy. 
Their results are not absolute and do not describe the spectra near the initial energy. To 
obtain more complete spectra, it is necessary to normalize their equations to a solution valid 
at high energies, and they have used the results of Rossi and Greisen [7] for this purpose. 

All these calculations are difficult to generalize to a more accurate treatment of all the 
different interactions, and hence there is a need for a new approach that will permit use of 


arbitrarily accurate cross sections. 
4. Analytic Treatment of the Problem 


In this section equations will be derived for the energy spectra of electrons and photons. 
No distinction is made between positrons and electrons, the two being considered as having 


4‘ The critical energy in a given material the energy at which the radiation losses of an electron equ il the ionization losse 
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identical cross sections and combining to make up an “electron” spectrum. The differentia] 
energy spectrum for electrons will be represented by /(/2)dE, which is the track length traveled 
by all electrons in the energy range £ to E-+-dE. I(a)da represents the photon spectrum simi- 
larly.® In all cases energies are in units of me’, and the electron energies include the rest energy. 

To set up equations in terms of /(//) and //(a), a time-independent source is assumed. 
Physically, this means a continuous cascade in which the average number of electrons, positrons, 
and photons at any given energy remains constant with time. The equations will be derived in 
a general form, and considerations of the interactions and cross-sectional forms will be reserved 
for a separate section. 

Continuity equations for /(/2) and //(a@) can be written as follows: 


dE’k* (E’ FE) Ib’) =s*(E) + da’ k®® (a! 2) H(a’) (1) 
I K 


| da’ k?® (a! a) H(a’) = 8?* (a) 4 dE'k**(E’ a)I(E’), (2) 
a Ja 


where 

ke*( 2’, F) is the total probability per unit path length traveled for an electron of energy EF’ 
to undergo an interaction that removes it from the energy range above EF, 

k**(o’,F) is the total probability per unit path length traveled for a photon of energy a’ 
to generate an electron of energy greater than £, 

s®(#-) is the total source strength for generating electrons with energies greater than EF, and 

k?°( a! a), k?®( Ee’), and s?*(a) are similar quantities with the photon and electron roles 
reversed. 

Note that the superscripts refer to “electron scattering,” “photon generating,” etc., al- 
though strictly speaking the interactions include types other than scattering. Physically, eq 
(1) states that for a time-dependent source the rate at which electrons are removed from the 
energy range above F is equal to the rate at which electrons come into the energy range above 
E. Equation (2) has the same meaning for photons. 

Equations (1) and (2) completely define 7(/) and H/(a) for specified sources; however, 
these equations as they stand are not very suitable for numerical work. It is advantageous to 
make further transformations suggested by the nature of the variations of spectra and cross 
sections. The particular choice given here is based largely on numerical and analytic ex- 
ploratory calculations that are not sufficiently pertinent to the main discussion to warrant a 
detailed description, but some general arguments to justify this choice are given below. 

I(f2) may vary rapidly if s**(/2) is discontinuous (i. e., if the source is concentrated at some 
energy), as for example, in the case of figure 1 in reference [8]. A more slowly varying function 
on which to base a step-by-step numerical integration Is 


oo 


F(E) dE’ I(E’). 
a E 


which gives the track length traveled by all electrons while their energy is above E. The fune- 
tion /7(a@) may also be singular if the source function s?*(@) is discontinuous, but photon cross 
sections permit the removal of such singularities by iteration (see, for example, reference [9]). 
It is therefore not necessary to utilize an integral form of the photon spectrum. 

Equation (1) could be rewritten in terms of F(/) by partial integration, but this results in 
the differentiation of k@(/’,F), which is already singular as -’ approaches &. <A further trans- 
formation is needed to make the equation more tractable, namely integration of all terms of eq 

1). With regard to eq (2), k?*(a’,a) and k?*(’,a) are integrals of nonsingular cross sections. 
It is therefore possible, and in fact advantageous, to accompany the partial integration of the 
last term in eq (2) with a differentiation of all terms in the equation. 


EF) and Ha) arein the unit en Mc Sec These are track-length units referring to a source of one electron, The dimensions of 
li} ind Illa n etotl e of flux if the source l electron per gram per second 
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Applying these transformations we obtain for the left side of eq (1) 


sa . . nd 
| dk’ | dE’ k*(E”’,E’) T(E”) dk’ Ik’ dE” k*"*(E’ ,E” 
e E e gE’ e k e E 
. *; , 
| dE’F(E") 5, |  dE"k(E',E"), 
JE ok Jk 
The other term in the integration by parts, —/ EN] | . dE’ k**(h’ E” ] » is zero at both 


limits. 


Equation (1) thus becomes 


dE’K®(E’ »E) F(R’ S*(f)-4 | da’ K®(a’ FE )H(a’), 2) 
eI a Ff 
where 


K*(E’,E)—~2 | dE" k*(E',E”), 


ke a’. FE) dE’h® a’ 5’ 
Jf 


E 


Equation (2) becomes 


lar Hia)— da’ K*(a’ a)H (a’ ay 8 a) T dk’ F (k’)/ Kk EF’ a), 4 


Ja Ja 


Note that the other term on the right side, 


O ’ 2 ad 
2 FUE’ )ke(E 0) ree), 


obtained after integration by parts, Is zero at both limits. 
The various quantities in (4) are defined as follows: 


2 
fa Oya 


is the total photon interaction probability per unit path length for photons at energy a 


, 


™ 9 ol” a _ P 
Ae a .a@ = : NZeo @ ,a@), (4b 
a 


where o(a’,a) is the photon differential scattering cross section per electron, and N is the 


number of atoms per gram. 


- : O7k?* (Ee a) O , ’ 
K”*(E’ a) =———,, —— No(E’,a), : 
s ' OE’ da on’ *** ; . 


where o(/2’,a) is the cross section per atom for producing radiation. 


OS?* (a 
td 


NS? aji--- 
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s the photon source strength ata 
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5. Interactions and Cross Sections 


Equations (3) and (4) as derived in section 4 make no reference to specific interactions or 
cross-sectional forms, and consequently this choice can be adapted to the specific nature of each 
f problem. In terms of the interactions discussed in section 2, the functions A may be written 


s(L » > 
kK E’ .E) Kor msstrahlung losses ~ Ktomic interactions, 


_ ae Pag > 
AK "(a .E) K pair production T Ke. mpton electrons» 


oth 
Ks a’ a K compton photons: 
Kr kh’ 1a Fi teemests ihlung photons» 
7 sla) M\ @) pair production tT M\@) Compton photons > \ @) absorption- 


At this point accurate cross-sectional forms could be introduced in the equations, but an 
examination of their suitability can be made with simplified forms similar to whose used by 
earlier investigators. Specifically, (a) extreme relativistic cross sections are used, (b) com- 
plete screening is assumed, and (c) the rest mass of the electron has been neglected insofar as 
it influences the limits of the integrals. Detailed cross sections are given in appendix 1. The 


K funetions derived from these cross-sectional forms are 


} K*(E.E vae{ | —}—j na wo athe 
4 | = n 
, }<) ~~ 
—jg79 IN A—a +t U- ip +eey 4°? i “| a | 


where n= Ek EI’: 


eae — 7 7 De Ste Ed a | n? 
K*(a’ EE) = N8&ea {EE g Its jtty 7 |in 1I83Z eee =4 


j r n' ung? * i—9, Pm i 1, ! 3__1 
oa YNZ ger (1 n) (InE+In —"+In2 5 +5 tg vz b 
where n Ka’: 
K**(a@’ a)=NZ = n 
ax 


ih) , 
where n7—a a 


. ” a 4 a rr 
K"®* (Ee a) = N4¢ | n°? +2n*) In 1832 or *} 
a o Y 
where n=a@ i’; and 


ula NAG E In 18327 3— ai |+ NZ , r+ | In2a+5 } 


6. Numerical Treatment of the Problem 


In determining a convenient system for the numerical integration, we note that the func- 
tions A depend primarily on ratios between the two energies involved. We wish to take 
advantage of this sealing property by choosing integration steps of equal length in the logarithm 
of this energy ratio. This choice of integration mesh has the additional advantages that one 

can integrate over an enormous energy range with comparatively few steps, and that a cross- 
sectional tabulation can be made that is applicable to widely different energy ranges. 
The numerical integration is performed by Simpson’s rule, using modifications of the 


familiar (2, 4, 3,4, 4)A set of weights to make the rule apply also to an odd number of intervals. 
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The integrals were computed essentially as follows: 


dh’ KR (E’ BE) FUE’ dn BE’) E’R*(E’ EF) FUE’ 
n 
— me 
=>) w',,L KO (E,,E,) FE 
n . 
2S) w,, ASN, EF) FU 
Note that the integration weights ,, decrease in size as the energy decreases. More details 
are given in appendix 2. 

Despite the transformations discussed in section 4, the function A*(/’,/2) in eq (3) is still 
singular. However, this singularity causes no difficulty if one writes the integral in the fol- 
lowing manner: 

dh’ KR ( bh’ ».E) FE dh’ KR’ (BE .E)[FPUE) —F(E)|+FUE dh’ Ke (EE 
Jf a i ae i 
. . . . . | 
Note that the first term on the right is no longer singular and can be integrated numerically, 
The last term can be evaluated without difficulty by integrating the singular parts of A®(2’,E) 
analytically. Notice also that as /’ >, the first integrand is finite, approaching the limit 
AlOF( EF) OF), where A=lim( kh’ — FE) A’ (E’,F). This term is very small and needs to be evalu- ’ 
E’ >} 
ated only approximately. A three-point numerical differentiation formula, described in appen- 
dix 2, was adopted for this purpose, the corresponding weights being denoted by a. By } 
utilizing these devices, we represent the first integral of eq (3) as the sum of terms 
FUE) [Ae CE >) w, KEE Aa,]+>3lw;,A° (EE) + Aa] P(E 
where 4° (EK dE’ RK (E’ FE). Note that a; is zero except for i=n, 0 lL, n—Z. 
« fk 
} 


If the funetional form of A@( Ee’ -) is that given in section 5, 


’ —- 1 4 } l l a l 
(kh V4 ok {[ a9 I—y)In(l—y rs l—n ren ta fn IS3Z IS 


where n=) ko, and £y is the highest energy of the source of electrons and photons. For a 
discussion of the choice of the cutoff value Q, see reference [8]. In this paper it has been ey aluated 


as follows: 


Q i je =4.9214<10-", 


where ] 13.5 ev. 
The function A?*(/2’,a) in eq (4) is not singular, but the exact expression varies rapidly as 


I’->a. This can easily be taken into account by using addition-subtraction as before: 


dE'K* Be’ a) FUE’ dE’ KE’ a) [FP UE"’) —F CE) peal +F(E) peg | dE’ RK" ( KE’ a). 


e x . 
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Written as a sum, this term becomes 
n—1 n—l 
kik KP? a, \— Sw, AX?" hk yt, t >a, (Ej ,a,) FE, ’ 
where 4 "(a dE’ KE’ a). Tf the functional form of A’*(/’ a) is that given in section 5, 


,—l| (4 4 , l 
KS a A Tal er Oe mn) In 1832-834 9 "| 


where n=a/ ky. 
Equations (3) and (4) can now be written in their final form for numerical evaluation: 


n l 7 l 
Fike 7% - | SOW LE) +> WE (a +S¢(E,) | (5) 
n l " l 
Hie ; , | F E,,,) D?®(a,) + SSW FE) + >We (ay) +88 (a, | (6) 
where 
iW w,, Ao (E, la , 





W Ww Kk 5(q fle 
«) 

n l 

C*(E,) = 4°(E,)— 2 onK"(E,,E,) +Aa,, 
? uo 

(? (a ula w, A * (a, a 
n—l 

it - VI" pe 

D?*(a,) =. 4 (a, )— Son pe. 

i= J 


These seven functions depend only on cross-sectional forms, and they are evaluated separately 
prior to the main calculation. Note that eq (5) depends on F(/;) and H(a,) but not on H(a,), 
while eq (6) depends on F(2,) and on F(£;) and H(a,;). Therefore, eq (6) must be determined 
after eq (5) at each step in the integration. The problem of how to start the calculation is 


taken up in the next section. 
7. Special Solution Near E, 


Fike) and F(a) must vanish at the highest energy, &). It would be possible to use eqs 
5) and (6), with the cross-sectional forms that have been listed, to solve for the spectra near 
k,. But it was found more convenient and accurate to do a separate calculation because the 
solution near /) has certain special features. In particular, we note that (1) for small energy 
losses it is always possible to use extremely simple cross-sectional forms, (2) the linkage be- 
tween the electron and photon equations can be simplified because the higher cascades have 
negligible influence in this energy region, and (3) it is possible, though not necessarily desirable, 
to seek analytic solutions. 

Jecause even an analytic development leads to a numerical integration, it appeared most 
expeditious to employ eq (3), with the simplifications just mentioned and a finer integration 
spacing, to develop a solution in this region. In particular, for an electron-generated cascade, 
the term for generating electrons from photons and the term A[OF(/)/OF] are very small and 
can be neglected. The remaining functions A can be written in the approximate form: 
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A : k’—E 
kk Bln BE +, | 


4B) =A In *_ BE) in 4+ C(E—E 
i) eE, 


c 


where A, B, and C' are constants (see appendix 1). 
- uy " ; ' 111] — 
lhe quantities (7) were tabulated at Ay 2 . Where m 9’ 4’ 9’ 1. This spacing required 
a special set of weights; but otherwise the evaluation in this special case is carried out by 


using the appropriate form of eq (5) in a manner similar to the main calculation. Results 


of this special solution not only give the solution near /, but also permit a more accurate 
evaluation of the integrals in eq (3) for a’,k’>k,x2-'°. In this energy region the use of 
Simpson’s rule with 27! spacing is not really a good enough approximation because F(E) 


has large nonquadratic variations 
8. Calculations 
The calculations in this paper are based on a monoenergetic source of 1 electron per second 
at 360 mc? in lead. Results were first obtained for the special solution near /y. The main 
calculation, extending down te 5.6 mec*, was carried out in two parts, using the mesh /)>2-" 


to the energy }/,) and the mesh /,<2~-"* from that energy to h=5.6 me’. | 
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ki RE | Arrangement of tabulated functions 


The calculation was facilitated by arranging a triangular tabulation of the four W’s 
and a linear tabulation of the C’s and J) in a form shown in figure 1. Most of the work went 
into tabulating A functions and Ws, the final computation being carried out in a few hours. 
The results are given in table 1. 

Any error had the effect of introducing an oscillation in the final computation, with alternate 
The integration procedure appeared stable in that such 


points being too high and too low 
By plotting the results it was 


oscillations tended to damp out slowly rather than to magnify. 
possible to use these oscillations to find errors. The final sequence of numbers is smooth, 
i. e., nonoscillatory, to better than | percent, which may indicate that the integration procedure 
introduced errors no larger than this, although the possibility that errors may have piled up 


at low energy is not excluded 
TABLE 1 Results 


E Ful Illa I FUE Ilia I Ful Illa 
“.0 wo. 4 7) Wo oN 4 2 4 1 44 
4. & aol ' ’ ve Wend ¢ 2 oS “ 
19 8 mi) 2 Hu 17.9 y, 

us 1. 35 4)4SS j l aul 14 1; \ 

$20.7 1. 7 M2 ; 115 11. 3 4] O.S 
285.7 2.47 Ooos5 m0 & 4! 14s &Y ; Ww] 
2h. ¢ a 13 i l 235 7.04 s 13 
TH 8 $ 1 Im) 7 13. 4 is4 ‘ ; 14.2 
202 14 124s 4 ne “il 
180. 0 71 0314 Ix. 9 O44 


9. Discussion of Results 


The results in table 1 are plotted in figures 2 and 3, where they are compared with the 
results of Richards and Nordheim |6] and calculations based on the analytic solution given 
by Rossi and Greisen |7]. (See appendix 3 for a discussion of the normalization.) For the 
comparatively low initial energy considered here, the equations of Rossi and Greisen are only 
accurate in a very limited energy region around 130 mec*. At lower energies, as pointed out 
by Richards and Nordheim |6], a more accurate representation of the interaction probabilities 
is necessary. The Richards and Nordheim numerical calculation was based on more accurate 
cross-sectional forms, but their results are normalized to those of Rossi and Greisen at 131 me’. 
Thus, near that energy, the significant comparison is with Rossi and Greisen. At lower 
energies, due to differences in the cross-sectional functions involved, the three curves go apart. 
At energies near /y no comparable calculations have been found. Thus the only energy 
region Where close agreement between the spectra calculated in this paper and those given by 
the other authors could be expected is from } to 4&5. In this region the agreement is shown 
in figures 2 and 3 to be quite good. 

Through a consideration of the specific cross-sectional forms used in the three calcula- 
tions, some further comparisons are possible at lower energies. For example, the spectra of 
Rossi and Greisen indicate a larger number of electrons and a smaller number of photons than 
reported in this paper. This is caused mainly by the cross section for pair production, which 
is significantly lower under the complete screening approximation than the asymptotic form 
used by Rossi and Greisen. Thus, in the calculations in this paper, low-energy photons do 
not generate as many electrons, and there results a decrease in F'(/2) and an increase in /1(q@) 
relative to Rossi-Greisen. The spectra of Richards and Nordheim indicate a much larger 
number of electrons and photons at low energies. Although this must be traceable similarly 


to a difference in cross sections, the reasoning is not as straightforward in this case. 
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FIGURE 2. Comparison of F(E) with previous re sults FIGURE 3. Comparison of H (a) with previous re sulta. 
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10. Final Comments 


The agreement as to general trends and as to absolute magnitude near }/) demonstrates 
the feasibility of the calculational procedure developed in this paper. Logically, the next 
step is to introduce accurate cross-sectional forms into the equations to obtain more accurate 
spectra. Alternatively, one could move immediately toward a solution of the spatial and 
directional distributions. This would involve computing the higher moments in the linked 
set of equations of which these spectra form the zeroth moment (for an example of this, see 
reference |10]). 


11. Appendixes 


Appendix 1. Cross-Sectional Forms 


The following detailed cross-sectional forms were used to derive the functions A’ discussed 
in section 5. 
Bremsstrahlung losses [11 


aes eS It+n? 2 as | 
o(E’.E rr a7 | 7 = | ( “sy —; )In IS3Z-!*4 i} 


where n= FE’. 
Electron atomic interactions [S 


where n= Fk EL’, 
Pair production (either electron or positron) {11 


: a ‘ l 
o,,(a’, Fb Sey : ( . 7. n+] )In 1I83Z—*- q 7 l—y } 


where n Ka’. 
Compton electrons: 


Compton photons 


where Ui) aa’. 
Bremsstrahlung pbotons [11 


" ~1[/4 4 z-441—9 
Op, ke , a ar | (5-3 Ul Ua ) In 183Z Tt 9 } 


The cross-sectional constants were evaluated us follows: 


where n=a lk’, 


© “ir Z- 3.9435 barns, Z=82; 


137) Z 


and 
7 = 4 0.6652 barn. 
The constants A. B. and (in the special solution near ke (section 7 ure obtained from the 


expression for A“(2’,F£); they are coefficients of the terms proportional to /7-(l—y 
In(l1—n), and 1, respectively, in the limit £,h’--+k): 
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Appendix 2. Integration and Differentiation Weights 


To obtain the integration weights @,, used in connection with eqs (5) and (6), the following 


formulas were used: 


daf r) == [flay tA) +-4f (a9) +f(ay—A)|] 
Jz,-4 


| ‘azt J a [Sf In7+A S/ La f Da A)], 


| def I a | flay +2 A) +t 3f(2)], 
a/ Tt, A 


where f(z) is a function that varies slowly for 7 in the region of integration. 
For differentiation the weights a, were determined by using the following formula: 


a ) n(E,—E)) 
=" be D gS! EO) T1(EL—E) 


Appendix 3. Normalization of Previous Results 


The following results were obtained by using a table of calculations given by Richards 
and Nordheim [6 


ke F(E HU S6F pyl Fk Ilia 5.329a(K 
131 6. 45 0. OO839 
91.8 OS _ a 
65. 6 12. 6 384 
5 + es 21.0 .. oa 
26. 2 30.0 2? 61 
19. 7 a4. 4 1. 74 
13. 1 19. 6 10. 7 
6. 56 tL 2 34. 1 


The normalization factors arise in the following manner: 


69.86 0.437 *-— density of lead (11.005) radiation length (0.529) 


5.329=—69.86 ~< enerev conversion factor 


13.11 


| 


For a discussion of the factor in square brackets, see reference [6], 
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The following results were obtained by using the equations of Rossi and Greisen [7] under 
“approximation B”’: 


I F(l ». 8227 ,( Ko, I Illa 2.9752, (Eo, W 
180 1. 77 0. O348 
ag 9. 02 . 132 
1 16, 4 {S82 
22. 5 IS. 5 1. 70 


The normalization factors arise in the following manner: 
5.822 dlensity of lead (11.005) radiation length (0.529), 


9 975 5.822 me~ (0.511) 
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Observer Differences in Color-Mixture Functions Studied 


by Means of a Pair of Metameric Grays 
Kenneth L. Kelly 


The Granville metameric gray panels have served to characterize the color vision of 

} observers in an approximate but useful way. Thus, the 39 observers studied, and the four 
color-mixture funetions were classed into five groups aecording to their indicated amounts of 

These panels were also used to measure the correlation of lens and 
macular pigmentations with observer age, sex, and eye and hair colors. Of the three sets of 

color-mixture funetions intended to refer to 2° field observation, the 1931 Commission 

Internationale de L’Eclairage standard observer agrees best. These data thus afford no 

| basis for supplanting the 1931 Commission Internationale de L’Eclairage standard observer 
with either the Judd “i” or the 1955 Stiles’ 2° color-mixture tunctions. The 1955 Stiles’ 


ocular pigmentation. 


10° functions, however, agree fairly well. 


1. Introduction 


The standard observer and coordinate system, 
now widely used in the interpretation of spectro- 
shotometric and colorimetric data, were recom- 
mended in 1931 by the CIE (the Commission Inter- 
stionale de L’Eeclairage) for this purpose so that 
ll subsequent data would be expressed in the same 
) wistimulus svstem and would therefore be immedi- 

ively comparable. This svstem was based on the 

work done by Guild and by Wright in England and 

til recently has proved quite satisfactory. How- 

eer, it was reported to fail in several instances to 
> «count correctly for the differences in color between 
two samples whose spectrophotometric curves differ 
nthe short wavelength end of the spectrum, espe- 
ially below 420 mu. 

These reported failures resulted in much discussion 
boking toward a possible revision of this standard 
bserver and led Judd to suggest in 1951 a modified 
vt of color-mixture functions based on the data of 
Wright and Guild combined with the standard lumi- 
nosity function modified below 460 my according to 
luminous-efficiency data by Gibson-Tyndall, Wald, 
Weaver, Thomson, and Ishak. Now Dr. W. 38. 
Stiles of the National Physical Laboratory in Eng- 
land has undertaken the first careful direct determi- 
nation of the color-mixture functions of average 
normal observers which the standard observer should 
represent. 


| 


black pigments called the simplex gray (No. 1), and 
a selective gray produced by a mixture of yellow, 
green, purple, and white pigments called the com- 
plex grav (No. 8). Figure 1 and table 1 show the 
spectral directional reflectances of these two grays 
obtained on the General Electric recording spectro- 
photometer. 

According to the 1931 standard observer, the com- 
plex gray will appear greener in daylight (source C) 
than the simplex, and the reverse will be true when 
source A (color temperature 2,854°K, representative 
of incandescent lamp light) is used. The same effect 
results from changing the angular subtense of the 
sample from 10° to 2°; the retinal area stimulated 
first will include the macula and also a portion of the 
surrounding retina, and second will lie inside the 
macula only. This change from 10° to 2° will have 
an effect at least partly analogous to placing a 
vellow filter (macular pigment) in front of the eye 
for the 10° condition and, thus should have some- 
what the same effect as reducing the color tempera- 
ture of the source. 

The purpose of the present study is to establish a 
criterion by means of which it will be possible to 
determine whether any set of color-mixture functions 
gives predictions of the character of the color differ- 
ence between the Granville grays in accord with 
observers known to have normal vision by the ac- 














In these discussions, use of a field size larger than 
the 2° field used for the 1931 standard observer has eee 
been advanced to accord more closely with viewing [ « 
eonditions in industry. Dr. Stiles is therefore mak- ° a . 
ing his measurements with both 2° and 10° fields ad 2 7° . “er 
and has already reported color-mixture data for a | ¢ ° Psy 
plot group of 10 observers whose average age Is 40 m . ° 
about 31 vr [8}.! - ce . 
In 1949, Walter Granville, then of the Container | &30F****eeees “eae 
Corp. of America, painted several metameric gray SIMPLEX GRAY 
| panels to illustrate the possible effect of angular sub- L © COMPLEX GRAY 
lense on color matching [2]. Two of these were l i ! 1 
4 500 600 700 


| chosen for the work herein reported; a nearly non- 
_ lective gray produced by a mixture of white and 


in brackets indicate the literature references at the end of this paper, 


Figures 


WAVELENGTH, my 


Spectrophotometric curves of the two metameric 


Granville grays. 


FiGuRE 1. 
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TaBi kb l, Spectral directional re flectance s of the two Granville 
grays 
Wavelength Simplex (No.1 Complex (No. 8 
3s0 0. S17 0. 30 
390 $24 $07 
we) $31 Lt) 
10 33s $10 
0 37 snd 
450) $30) ut 
40 $4 1) 
1) +4 uM 
1) $44 sit 
470 $32 ‘734 
iM) 331 18S 
wo 340) Js 
“N) 0) 1u7 
510 $31 HI 
2 $31 $41 
i) 42 $57 
‘40 44 $37 
SM) $58) sil 
vn) tim ™U 
570 42 | 
wt) $43 s7 
YO 4 ine 
How sit) He 
ole 347 tSU 
620 M47 2 
i) 347 52 
Ho “7 i~) 
OM) 47 “M4 
twa) 44s ‘17 
O70 447 l 
Oso) $47 “iH 
oun 44s is4 
Tuo $48 iS3 
71 44s 7 Le 
72 47 2s 
7a) $40) it 
7) $46) su 
7) 4 tt 
7 s44 27 
770 $43 27 
FE xtrapolated, 
cepted tests. In particular, this criterion is to be 


applied to the four sets of color-mixture functions 


already mentioned. 


2. Experimental Method 


The experiment consisted of placing the two gray 
panels side by side in a nearly vertical position illu- 
minated by CIE source A and viewing them through 
a Davis-Gibson (CIE) C filter prepared i February 
1955. The spectral transmittances of this filter de- 
termined in July 1956 at the start of the observations 
are given in table The observer 
stand at a line that would cause the retinal image of 
the two panels to subtend a 10° angle and to dis- 
regard the Maxwell spot |7] if present. First he was 
to describe the color of the simplex gray with respect 
to that of the complex gray and then, as the operator 
reduced the voltage on the lamp and thereby its 
color temperature, note when neither panel ap- 
peared redder or greener than the other. The volt- 
age of the source at this point was recorded along 
with the observer’s age, sex, and hair and eve colors 
Observers not passing the 5th Edition Ishihara Tests 
for Colour Blindness were eliminated 


was asked to 


to 
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PaBLE 2. The spectral fransmittances and change in rec 


, ‘proe 
color te m pe rature for the da jlight filters used ' - 
W ive Dav s-Gibson Filter G90A Davis-Gib- Filter No Davis-Gip. 
length ( 2.13 mn son plus S331 
son plus 
C7404 No 8331 
Mpa 
~~) on 0. 751 ( Is. )47 0 > 
400) 71 774 49 4 r 
we 72S ba) “ OA 40% 
10 TOL sit Hom 712 534 
120 7475 TSO S830) Hon 190 
4130 7IS0 7585 dit ol7 443 
140 HH23 vu S21 : $i 
ir SS Hus. 75 {UN 9] 
feu) “a0 Wz {384 40) 4 
170) $557 Has 2801 38S ba 
iN) pone HR 244 $44 140 
140) S532 77 2040 wy 108 
TLL a Al 14} wt) ONIN 
510 2479 5323 1320 254 0630 
520 2178 OSS 1108 22s 497 
530) AV27 $S0)5 vu74 1vy 0403 
40 1961 102 000 179 0351 
Mw) ISU] 1547 O60) 172 032 
tM) 1765 1600) OS12 17 0311 
570 1600) $500) 0730 171 0274 
cSt) 1420) wi mol ™ 0213 
yuu) mM) Sum) O44 114 0149 
wn 11145 $805 O424 14 0129 
ole Ol 37™) USS] 103 O105 
H20 a40 SHO (344 (4 ONS 
630 OSS0) 4448 O30 OSS OO73 
mo os. $318 (274 O72 oso 
OM) o774 $4) O20 (Wit) OO5! 
we O724 $217 (233 OS 0049 
670 OOs4 S440 (228 72 0049 
ont) wis4 145 1 O7S (O49 
nun O5s84 tis ri OSS (O4N 
TOO) O534 ee ola ev O44 
7 (484 wit oly] m™ O89 
TA 444 {0G O1SS OTH O03 
~34) 404 00 0142 O73 (20 
740) 0364 $403 OLS OHS 0025 
7 O344 +418 OLDS Live La 
i ma {37 i l (Ww) (wre 
7 Tk ; OLO7 (Wid O00 
Sywra 22 713.9 st) 4 IS7 300. 0 
Extrapol 
The observer was then asked to step back to 
another line from which the retinal image would 
subtend an angle of 2° and the experiment was 
repeated. The combination of source and Davis- 


Gibson filter resulted in : 
ture of 50° A. 


uw source of higher 


1 maximum color tempera- 
but for those observers requiring 
temperature, more 
measured blue daylight glasses (table 2) were used 
with the Davis-Gibson filter. Therefore, by varving 
the voltage on the lamp and by using combinations 
of these filters, color temperatures could be obtained 
from 1,800° K up and, indeed, it was found that some 
sources bluer than a source corresponding to infinite 
color temperature were required The voltage and 
filter combinations were noted for each red-green 
balance point for each observer and the correspond- 
ing reciprocal color temperatures were determined 
from table 3 by subtracting the change in micro- 
reciprocal degrees Kelvin produced by the 
filter combination from the reciprocal color tempera- 
ture of the bare lamp. 
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Color temperature (0) in K of the lamp, and the 
temperature in urd of the lamp and the lamp 
more of the daylight filters against voltage on the 


TABLE 3. 


afi 
reciprocal color 


with one or 


lamp 
108 @ of 10°08 of 1060 of 10/@ of 
6 of imo lamp and lamp and lamp and lamp and 
Volt lamp of Davis Javis- davis- Javis- 
slome lamp Gibson Gibson Gibson Gibson 
” filter with G90A | with 8331 | with G9OA 
and 8331 
36. 0 1800 a1) rs | 27 16 SA 
2.7 1900) 12H 324 24 137 os 
49 9 O00) 500 oO HM) 110 3 
6 2100 176 274 16 St) s 
i 2200) 4155 ae 174 Hs i3 
=5 2300 455 233 14 45 —33 
a4 2400 417 14 137 27 51 
“4 2500 100 IVS 119 10 Os 
4 2600 385 IS2 104 ) AS 
1S 2700 370 16s Ww 20 us 
1% 280K) 357 15 vi 33 111 
133 284 350. 4 148 70 10 118 


3, Red-Green Balance Points for the Four 
Color-Mixture Functions 


The red-green balance points for the four color- 
mixture functions were determined by first com- 
puting the chromaticity coordinates for the two gray 
panels for each of the functions for three Planckian 
sources (277, 221, and 191 urd), for nine Planckian 
sources (666, 571, 500, 488, 486, 444, 400, 350, and 
combined with the Davis-Gibson “C”’ 
filter (A10°/@—202 urd), and for the same Planckian 
sources combined with the double filter Davis- 
Gibson “C”’ plus G9OA (AL0°/O=280 urd). These 
data were plotted ov the (x,v)-chromaticity diagram 
and fall close to the Planckian locus. 

It was noted that the direction of the straight line 
connecting each pair of points as well as the distance 
between them (indicative of color difference) varies 
regularly with the correlated color temperature of 
the energy reflected by the simplex gray regardless 
of whether the computation is based on a Planckian 
source or on a Planckian source modified by the 
single or the double filter. Figure 2 shows the 
Planckian locus and 10 of the 21 pairs of computed 
points plotted on the (x,v)-diagram for the CIE 
standard observer. Note that near 300 urd the 
simplex gray is seen by this observer as greener than 
(148 urd), the 


308 urd) 


the complex, whereas at source C 
reverse Is true. 

The red-green balance points for the two Stiles 
color-mixture functions were easily determined be- 
cause they each predicted that under a certain one 
of the illuminants the two gray panels would be a 
chromaticity match; that is, the chromaticity points 
for them were found not to be significantly different. 


The reciprocal color temperature of the source 
satisfying this condition for the Stiles’ 2° color- 
mixture functions is 210 urd, while that for the 


Stiles’ 10° color-mixture function is 267 urd. 

The other two color-mixture functions predict 
that at the red-green balance point there will be a 
residual yvellow-blue difference. Therefore it was 
Necessary to establish some criterion to determine 
when the red-green difference equals zero. The 


99 

















+ 
a 
’ 
400PF * 
° 
° 
“. 
¢ . 
. 
o “ 
@ ° 
*e, 
“€, cod ° 
a 
> °.“s 
° 
. 
y ssof _# 
6 * 
| €o 
¢ 
. 
<m 
2/ .& ” 
4 
. +, 
* Ss 
+ te 
* 
“ 
‘e 
300 9 SIMPLEX GRaY 
. 
© COMPLEX GRay 
% 
= = — —— i 
30 35 400 
x 


Ficure 2. Jilustration of the change in relative position and 
distance apart of the CIE (x,y)-points of the Granville grays as 
the color temperature of the source is changed. 


first criterion used in determining the corresponding 
reciprocal color temperature of source was that the 
line joioing the chromaticity points of the two grays 
be parallel to the yellow-blue axis on the (x,y)- 
diagram. This axis, having an inclination of 45°, 
was drawn between the 476 and 578 my points on 
the spectrum locus. These are the proposed centers 
of the yellow and blue segments of the boundary of 
real colors on this diagram [6]. The corresponding 
values of reciprocal color temperature of source 
satisfying this criterion for the 1931 CLE and 1951 
Judd “i” color-mixture functions are 192 and 202 
urd, respectively. 

The second criterion used in determining the 
corresponding reciprocal color temperature of source 
was that the line joining the chromaticity points 
of the grays be perpendicular to the red-green axis 
on a uniform chromaticity scale diagram [4]. This 
axis was drawn from the 508 my point on the spec- 
trum locus to the intersection of the line drawn 
through the 493 mu point and the “C” illuminant 
point with the line connecting the red and violet 
ends of the spectrum locus. These are the proposed 
centers of the red and green segments of the boundary 
of real colors on the (x,v)-diagram [6]. 

A perpendicular to this axis was drawn through 
the “C” illuminant point and was found to intersect 
the spectrum locus at 567 mg. This perpendicular 
was transferred to the (x,y)-diagram by drawing it 
also through the “C” illuminant point and the 567 
my point on the spectrum locus. The angle of in- 
clination of this perpendicular is found to be 66.5° 
and the corresponding values of the reciprocal color 
temperature of source satisfying this criterion for 
the two color-mixture functions are 206 and 213 urd, 
respectively. Table 4 summarizes the values of 
reciprocal color temperature computed by these 
criteria to correspond to these balance points for 
the four sets of color-mixture functions. 








TABLE 4. 


halance points for the four color-mirture functions by diffe rent 


Values of rec procal colo le m pe rature of red-qreen 


crileria 


Reciprocal color ry ire of 
Color-mixture functions source in grd yielding red-green 
balance 
By chromatic h 
Stiles’ 2 21 
Stiles’ 10 207 
Parallel to Y-B Perpendicular 
xis on UCS triangk 
to R-G axis 
1931 CIE Ww 206) 
1951 Judd “i” 202 213 


4. Results and Discussion 


Table 5 contains the individual data for the 39 


observers studied with their descriptions of the two 
panels and the reciprocal color temperature of their 
red-green balance points. Note that observer 39 


TABI E 5 Observer initials, age, ser, eue, and haw color, de scripl 


g positions, and rec iproc al color te mi pe rature of source requir 


di flere nee hel we en the nei 


Sex and 
Observer ize Eve ce Hair « 1 
DBI NESS blue brow! 

2 KLK M45 hazel brown 

; GWH P56 brown brown 

4 MRI M27 hazel lark brown 

CAD M44 blue brown 

th MMB F39 blue lark brow1 
rs LEB MO5l dark brown black 

5 rHe Mal blue blond 

’ IN M3s blue red 

10 Ccpc M21 hazel suburr 

11 F50 hazel brown 

12 HKH M39 blue blond 

13 FCB M62 blue brown 

14 RT\ M34 blue brown 

i RWC M46 blue brown 

16 WAH M29 dark brown black 

17 GHL M34 dark brown black 

18 WHA M47 blue brown 

1 RLM F44 brown dark brown 
20 WRD M42 brown brown 
21 ECW M48 brown brow! 
22 BPG F37 dark brown dark brown 
3 WBE F23 brown brown 
24 SAC Fig brown brown 
25 IJ F29 blue brow! 
2h CAL F25 hazel light brown 
77 RPI M53 blue brown 
28 HJK M52 brown brown 
20 VIB F42 dark brown dark brow1 
$0) KSG M66 blue brow! 

l RD \I68 blue brown 
$2 WLH M21 hazel blond 
$4 WEM M36 blue blond 
44 PI FIs hazel brown 
$5 ICs M28 green blond 
st rol M3 brown light brow: 
$7 Ks NI56 blue browT 
38 IWL M31 dark brown black 
3Y CB M77 blue brow! 

A verages 41.3 








was able to obtain a match at the 10° position by 
that it was not possible to run the lamp at a suff. 
ciently high color temperature for him to obtain A 
match at the 2° position, presumably due to his ia 
heavy ocular pigmentation. ' 

Because the image of a 2° field talls wholly Withi: 
the macula, the 2°-field results depend both on the 
macula and on the lens pigmentations. 
since the image of a 10° field covers the macul 
a large region of the surrounding retina of which the 
area of the macula represents about 10 percent, the 
10°-field results depend only on the lens pigmentation 
provided the observers disregard the Maxwell] spot 
as instructed. Also, the difference between the two 
may be taken as a measure of the macular pigmenta. 
tion on the assumption that the spectral sensitivities 
of the receptors apart from pigmentation are the 
same throughout the 10° field {7}. 

The degree of pigmentation of an observer is 
indicated by his red-green balance point. By this 
interpretation, Stiles’ 10° average observer (balance 
point at 267 urd, see table 4) would be the least 
pigmented of the four. 


ons of colors Huminated by source C 


al the 10 and 
ed by each observer for red-qree n balance at both positions and the 


of aray panels ? 


Color of simplex relative t mplex Reciprocal color temperature of 
match point 
10 - 1 2 Difference 
red green 24 v1 113 
red green 23S 103 135 
red green 179 17 132 
pink pink Ts 14 x2 
red green 1 “4 112 
red green 213 132 SI] 
slightly pink green 213 7s 135 
pink green ith 144 107 
pink green an) 110 wo 
pink green 238 “4 14 
pink slightly pink O76 ISU 67 
pink pink 250) 179 71 
pink green 172 11S 4 
red green 222 135 Ny 
pink green 213 71 142 
pink green 2s IS2 aD | 
pink pink ON} 21 73 
match green ISY 65 124 
pink green 2358 159 79 
pink match 227 149 78 
pink green 2292 17 47 
lavender green 227 127 100) 
pink pink Jt) Itt v2 
lavender lavender 270) 21s ti2 
red match 278 154 124 
lavender match >) 170 80 
match green a5 14 109 
match green 170 2 168 
pink slightly pink 25) 175 81 
green green 4] 1) 131 
green green 130 2s 15s 
red green 233 110 123 
red green 106 2 167 
pink pink 278 170) 108 
pink pink 270 161 109 
pink pink 2s 192 71 
red green 146) M 112 
pink pink 2 182 74 
green green 118 11S 
214 114 100 
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Re ciprocal colo te m pe ature of source required im 


require lower reciprocal color temper 


the reciprocal color temperature 
red-green balance points for the 39 
observers studied table 5) are plotted with 
observer number against observer age. Since no 
we has been given for the 1931 CLE and 1951 Judd 
“ color-mixture functions, the spread in reciprocal 
color temperatures required by them for red-green 
} balance as given in table 4, are indicated by a pair of 
lines which should extend across the graph but have 
been shown as short sections only to avoid compli- 
eating the figure. Stiles gives the average age of 
his pilot group of 10 observers as 31. It may be 
seen from figure 3 that only four of the 39 observers 
studied seem to have less pigmentation than the 1931 
CIE standard observer, only two less than the 1951 
Judd “i? (both of these on the basis of vellow-blue 
criterion) and only one less than the Stiles’ 2° color- 
| mixture functions by this test. It also indicates that 
the CIE and Judd color-mixture functions both 
| correspond to very young observers in regard to 
halance point for the Granville grays. 
The rank correlation (Spearman's) between 10°/0, 
for red-green balance and age of the observers at the 
) 2° position is 0.60 with an uncertainty of plus or 
minus 0.22, the uncertainty being 4.9 times the 
probable error. Thus, it is suggested that there is a 
poor but significant correlation of lens and macular 


» 
oo. 
») 


In figure 
| (10°0,) of the 


{see 





| 
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pigmentations with age. Figure 4 shows the correla- 
tion of 10°/O, with age; the rank correlation ts better 
being 0.76 plus or minus 0.13. This correlation is 
attributable to the known increase in lens pigmenta- 
tion with Here again, there are only six 


age. 


$0 60 80 


OBSERVER AGE 


a 2° field by 39 observers to produce red-green balance between the 


Granville grays plotted against observer aqe 
{ } g f 


itures rhis tendency suggests that the amount of lens and macular pigmentations 


observers apparently less pigmented than Stiles’ 


| 10° color-mixture functions. From figure 5 it will 
be noted that there is no correlation of macular 


pigmentation with age, the rank correlation being 

().145 plus or minus 0.61. Also, there is no signifi- 
cant correlation of macular or lens pigmentations 
with either eye color or hair color. 

From table 5 and figures 3, 4, and 5, it is possible 
to classify the oberservers and the four color-mixture 
functions into five groups, depending on their 
amounts of ocular pigmentation. From table 5, it 
will be seen that certain observers describe the sim- 
plex gray as redder than the complex at both the 10° 
and 2° positions. These are the young observers 
and those tentatively considered to have the least 
ocular pigmentation and constitute group 1. Those 
in group 2 apparently have a little more pigmentation 
and so call the simplex redder at the 10° position but 
a match at the 2° position. Group 3, comprising 
the majority of the observers, contains those who 
call the simplex redder at the 10° position and greener 
at the 2° position. With more pigmentation, the 
observers in group 4 call the simplex a match at the 
10° position and greener at the 2° position. The 
most heavily pigmented observers, and these contain 
the oldest observers, call the simplex greener at 
both positions. Thus, these groups are called red-red, 
red-match, red-green, match-green, and green-green. 
The same grouping may be made with respect to the 
color temperature or reciprocal color temperature of 
the red-green balance points at the 10° and 2° 
positions. The three criteria for the five groups are 
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shown in table 6. According to this classification, 
all four color-mixture functions fall into the red-red 
group of very lightly pigmented observers with the 
CIE the most heavily pigmented and the Stiles’ 10 
the least. 


Granville grays 
one observer of 
determination 


temperature for 


of 


has been found to vary widely from 
normal color vision to another, and 
this value of reciprocal color 
any one observer serves to charac- 


terize his color vision in an approximate but useful 


5. Conclusions 


way. 


The value of reciprocal color temperature 


field for red-green balance of the 


required in a 10 
Granville grays is dependent on the amount of vellow 
pigmentation in the lens of the observer's eye, and 
1 measure of this pigmenta- 


The reciprocal color temperature of the source 
required to produce red-green balance between the 
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eful 
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ita- 


Classification of obser 


[ABLE © : 
the simplex gray re lative to the complex gray for source C, and alternatively by the color te mperature and reciprocal color tempera- 


ire of source required to make neither gray redder or greener than the other 


Obs may Simplex relative to complex 
10 2 

red-red red red 

ed-matel red match 

red-green red green 

matct ! mateh green 

reer ‘ green green 


‘ion. Similarly, the value of reciprocal color tem- 
perature required in a 2° field for red-green balance 
of the Granville grays is dependent both on lens 


wan) 


pigment and on macular pigment, and the difference 
in these two values (value for 2° field minus value 
for 10° field) is tentatively taken as a measure of 
the macular pigmentation of the observer. 

The tentative measure of lens pigmentation af- 
forded by the Granville grays correlates well with 
the age of the observer, but the tentative measure 
of macular pigmentation shows no correlation with 
ave. Neither measure shows significant correlation 
with eve color or hair color, but a tendency, signifi- 
cant at about the 15 percent level of confidence, was 
found for female observers to have less pigmented 
eves than male observers of the same age. 
~The reciprocal color temperatures of source re- 
quired by four sets of color-mixture functions for 
red-green balance of the Granville grays have been 
computed and found to be higher than the average 
value for 39 actual observers with normal color 
vision. Of the three sets of color-mixture functions 
intended to refer to 2°-field observation, the 1931 
CIE standard observer agrees These data 
thus afford no basis for supplanting the 1931 CIE 
standard observer either with the Judd ‘1’ eolor- 
mixture functions or the 1955 Stiles’ 2° color-mix- 
ture functions for a pilot group of 10 observers. 

With respect to these observations on the Gran- 
ville gravs, the Stiles’ 1955 data for 2° viewing does 
not fulfill the requirement recommended in 1955 by 


best. 


ers into five groups according to indicated ocular pigmentation as shown by their de scription of 


Color 


temperature of | Reciprocal color tempera- 
source at match point ture of source at match 
point 
0 0 6,750 148< 106 OQ» Woe 
0 6.750° =0 10 148< 10°/0,; 
0 6,750 0 10°/Oy 148< 10/0, 
a) 6,750 is) 1we/e 148 = 10°/0, 
6,750 0 0 w/e 10° OQ; 148 


the CIE [1] that any revision of the “standard ob- 
server for colorimetry should represent average 
normal vision, adjusted as for an observer 30 years 
of age’’, but the data for 10° viewing are fairly satis- 
factory in this regard. 
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Precise Evaluation of Surface Area With Indirectly 
Calculated Dead Space 


William V. Loebenstein 


The determination of surface area from nitrogen adsorption data obtained at low 
temperature by volumetric methods is subject to several errors. One error that becomes 
especially serious when the surface area is small is associated with the calibration of the dead 
space with helium, This error is large because the volume of adsorbed gas is calculated as 
a small difference between two comparatively large quantities. 

\ statistical method of successive approximations has been developed for calculating 
the dead space. The latter is obtained as that quantity necessary to vield the best agree- 
ment with the known isotherm equation. This procedure fixes the value of V,, (directly 
proportional to the surface area). The procedure also provides an estimate of the experi- 
mental error associated with the over-all determination, 


1. Introduction In routine determinations where the isotherm 
equation is known from previous experience, the 
Adsorption of vapors at temperatures near their | dead space need not be experimentally obtained. 
boiling points provides the most convenient means | The surface area can be found by assigning a value 
for determination of the surface areas of finely | for the dead space that will allow the experimental 
divided or porous solids. In volumetric apparatus, | measurements to fit the known isotherm. This paper 
which is most commonly used, the amount of ad- | describes a procedure of successive approximations 
sorbed gas is the difference between that admitted | by which this is accomplished. The conservative 
to the sample container and that remaining in the | technician may always prefer to determine the helium 
dead space. The latter is the void space between or | calibration directly. He can use the resultant dead 
within the solid particles plus the volume of the | space value as a starting point which may then be 
connecting tubing. “refined’”’ according to the technique described here. 
The precision with which a dead space calibration | It is unnecessary to continue the refinement beyond 
is normally made using helium is fairly well limited | the point in the calculations where the error of the 
by the equipment and does not depend upon the | dead space ceases to be the dominant factor contrib- 
surface area of the sample. The effect of this limi- | uting to the over-all error of the surface area deter- 
tation of precision upon the measurement of the | mination. 
volume of gas (nitrogen, for example) adsorbed is It is fortunate that the adsorption isotherms which 
especially serious when the surface area is relatively | apply to a large number of materials measured at or 
small. This is because the amount of gas adsorbed | near the boiling point of the adsorbate are expressible 
is calculated as the small difference between two | in one linear form or another. These are represented 
large quantities. Indeed, the uncertainty inthe | by the well known ‘Free Surface” equation of Bru- 
measurement resulting from the dead space calibra- | nauer, Emmett, and Teller [{1]' (BET) (eq 1), and 
tion is greater than that of all other factors which | by the equally familiar Langmuir equation (eq 2). 
contribute to the over-all error of the surface area 
determination. It is difficult to express this fact 
using rigorous mathematical language because the 


surface area evaluation results from the parameters PLV (po— p)I=MKE— P/V a) p/pot1/(Vnl) (1) 
of a least-square line assumed to represent all of the 
points. No simple relationship exists between the p!V=pol(VmB)+ p/ Vm (2) 
over-all error and the errors associated with the 
determination of the individual points, especially p=pressure attained when adsorption and 


since an error in the dead space actually destroys the 


as . desorption processes are at equilibrium, 
linearity of the plotted curve. 


Po= Vapor pressure, 


Anvone who has made surface area determinations V—volume of gas adsorbed STP 
in the conventional manner can calculate from his 2 Omscsmetenin (ander conditions of the experi- 
own data the percentage change in surface area re- : meant \ elated to the mene of adsorp- 
sulting from an error of, sav, one percent in the dead tiem euel = 
space. The author has found that this error exceeded \ a conatends equal to the volume of gas 
100 percent for a sample of glass beads whose area | = (STP) required to cover the surface to 
was of the order of 0.2 m*/g as measured with N, at the extent of one molecular laver. 


90° K. Normally, of course, the careful technician 
would in such a case make sufficient replicate de- 
terminations of the dead space to insure that its 
precision was appreciably better than one percent. Figures in brackets indicate the literature references at the end of this paper. 
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2. Determination of the Dead Space 
Factor Z 


The volume of adsorbed gas may be calculated 
from measurements as follows: 


\" 


~~ 


k [N—(fp/R) Z| 
where NV is the total number of moles of gas admitted 
to the adsorption tube and the quantity (fp/R)Z 
represents the moles which remain in the gaseous 
state in equilibrium with the adsorbed phase. The 
factor f corrects for the nonideal behavior of a real 
gas at low temperatures. It is usually expressed as 
(l+ap) where @ is a constant [2] depending upon 
the particular gas and its temperature. The quan- 
tity ? is the gas constant in appropriate units; Z is 
the ratio of the dead space volume to the absolute 
temperature; and the proportionality constant /£ con- 
verts the moles of gas adsorbed to milliliters per 
gram at STP. The problem then is to choose a 
value of Z by trial so that with the measured values 
of \ and p the best fit to the known isotherm equa- 
tion is obtained. 

It is customary in using the BET isotherm equa- 
tion to plot p [V(p p| as ) against PP) Po as NX. The 
resulting graph is linear and may be written as fol- 
lows: 

v=a+bXx } 
where the intercept a and the slope h depend only 
on Cand V,,. It should be noted that 1/V,,=a+6. 
The surface area is determined by multiplying V, 
by a packing factor characteristic of the gas adsorbed 
and of the temperature at which the isotherm was 
carried out. 

The best value of Z in eq (3) is obtained when the 
resulting values for the gas adsorbed obey eq (4). 
It is not possible to determine Z with sufficient pre- 
cision by graphical means alone, hence an analytical 
method was developed. This was achieved by con- 
sidering eq (4) as a special case of the more general 
empirical equation 

Y=at+bNXN+eX (5) 
The least square value for the added parameter ¢ can 
be determined and the value assumed for 7 adjusted 
to make this coefficient approach zero. 

The least square solution * for the coefficient ¢ is 


S28 27y— Sr Sry 
Sz*Sz* Sz’)? 





The denominator of eq (6) contains only terms * in r 

and, hence, is independent of Z and need not be 

determined. The numerator is designated by F, 
F iz) =S82?Sa°y— SiS, y (7 

As each successive choice of Z is made. the Value of 
Y is computed for each experimental point and F, 
evaluated from eq (7). If the new F,.,, is nearer 
zero than the last, the new choice has then resulted in 
an improvement in linearity.* 

The question may arise as to how far this pro- 
cedure need be carried before the results can be con- 
sidered acceptable. This can be answered in terms 
of the experimental error as measured by the random 
scatter of the points. This is discussed in the next 
section of this paper. A standard deviation s,,, 
attributable to random scatter is relatively insensi- 
tive to small changes in Z. On the other hand, the 
value of 1/V,,, may be extremely sensitive to changes 
in Z. Consider, therefore, two values of Z that 
produce corresponding values of /\,,, both greater 
than and less than zero, respectively. Two sets of 
values for ) may be determined and from these, two 
values of 1 Wii If the difference between 
latter two is small compared with s,))y, 


these 
, estimated 
from either set, no further improvement in Z can be 
justified by the data. 


3. Estimate of Experimental Error 


The preceding section was concerned with the 
contribution of the error of the dead space to the 
over-all variability of the results. When this is 
sufficiently reduced, the variability that remains is a 
measure of the scatter of the experimental points 

The errors associated with the regression param- 
eters a and b may he readily calculated. Further- 
more, an equation may be derived [8, 4, 5] for the 
standard deviation of the sum of a and 6. The result 
may be written 


2 [i/n+ (1 —NX)2/Sez’]. (8 


The term s2., is evaluated as follows *' 


(ny Ky! Sy- (Sry */ Sr.“ (Q 
Confidence limits in 1/V,, may be determined by 


inserting “student-t”’ values into the relationship 


(10) 


is the independent variable, is 
this is not strictly true, the 
small as compared 


, which is taken 
W hil 
since these errors are probably very 
is the volume of gas adsorbed, which 
i the ordinate Furthermore, the errors in Y are assumed to be com- 
pletely random and normally distributed about a mean of zero 

4‘ The quantity Fz) (as it passes through zero) has a positive slope (with respect 
to 7 An estimate of the magnitude of this slope from previous choices aids In 
the selection of 7 

’ The quantity (n—% 
number of degrees of freedom than (n—2 This i 
utilized for the determination of three parameters inst 

6 In computing sy.z from eq (9), rounding off by dropping of 
postponed until the end of the operation, otherwise serious errors 


troduced 


ire (or relative pressure 
i to be measured without error 
issible 


for derived quantities sucl 
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confidence interval surrounding V 


the Sep 
proportional to the surface area) is given as 


m 


Hence, | 
which 1s 


l . l 
/V.+l0/Vn)~'"* L(1/Vi») (11) 
The 95-percent confidence limits are implied if not 
aherwise stated. In layman’s language this means 
chat the interval determined by the stated limits 
yi], on the average, bracket the correct value of 
iV, 95 percent of the time. The appropriate 
‘student-t’’ values may be found in any elementary 
»xtbook on statistical methods. 

The method outlined will now be applied to a 
44188 g sample of Pitch Coke (Koppers Co.) which 
yas found to have a surface area of 1.6 m?/g. Table 1 
ymmarizes the adsorption measurements using 
vitrogen at 90° K for which the BET free-surface 
aquation was known to apply. The values in the 
cond column of table 1 were determined with use 
3) to express Vin terms of Z. Values of some 
917.52 


of eq 


of the quantities used in this example are k& 


> » 2000 ; 1 (OR)\-1 
ml mole; 2 = 6.236 10* (mm Hg) (ml) (mole)~! (°R) 
1+3.3 p * 107°; and 2,715 mm < po < 2,722 mm. 
TABLE | -ldsorpt on measurements of Ns; at 9U K on Pitch 
Coke Kopp rs ( o 
Relat f 
bere ! Vip p 
Y 1 n terms of Z by eq (3 
art 1700) 19. 367 42.670) Z 
14 1.7909) 112.469 62.366) Z 
17% 2.3364) '114.942—(78.231) Z 
1972 ” 8769) 116.2900 87.508 7) 
IsZ 3.0430 [17.893 06.916) Z) 
341 $3322) TIS.870 104.080) Z 
40 $7118) 120.375 113.057) Z 
79 3. ONS87) /[21.243-—(119.425) Z 


The geometry of the adsorption vessel furnishes a 
means to estimate the first values to be ascribed to 
Z. After a few trials the choice of Z was narrowed 
down to within 0.135 and 0.137 for which Fz, was 
less than zero and greater than zero, respectively. 
An estimate of the standard deviation was made ac- 
cording to eq (S) and the results are entered in 
table 2. 

It can from table 2 that the difference 
between the two values for 1/V,, for the first compari- 
sn amounts to more than three times the estimate 


be seen 


f "hy: Poar _ _ : 
Of Sar» This means that further improvement in 
TABLI 2 Res ‘its of successive appro. mations 

| 1 2d nparison 
/ 135 0.137 0. 136% 0. 1363 
} Wee Ol ny mr 
LV, ‘ist 2 Si o 7% 2 79] 
AU V, 4 0. 009 
Vim 0.047 0.044 0.049 
1204 1269 


Z should be made in order that the limiting factor in 
the precision of the experiment be the random scatter 
of the points. It is significant that if the dead space 
had been calibrated directly, and the results ob- 
tained were 12.15 and 12.33 ml for duplicate deter- 
minations, the agreement would have been con- 
sidered acceptable. Yet these values correspond to 
0.135 and 0.137 for Z which has just been shown to 
have too wide a spread to be acceptable. 

Upon further refinement of the computations by 
successive approximations, the spread in Z was de- 
creased to fall within 0.1362 and 0.1363 and the 
results are also summarized in table 2. The differ- 
ence between the values of 1/V,, has now been di- 
minished to an acceptable value of 0.009 since this 
is only about one-fifth the standard deviation of 


1/V,,. Hence, there is no need for further refine- 
ment. The average of the acceptable values of 


1/V,, is 2.786 with 95-percent confidence limits of 
+0.1267. This corresponds to a surface area of 


-. 7 0.080 . 
1.656 _. Mm 
0.072 


cr 
a 


4. Application to the Langmuir Equation 


The development is simpler in the event that the 
isotherm applicable to the data is the Langmuir 
equation (eq 2) where now 1/V,,—6. The quantity 
pie when plotted as a function of p results in a 
straight line providing the correct dead space is 
used. The experimental error is obtained by the 
same reasoning that was used in the previous section. 
Equation (9) is equally applicable for determining 
s?,. However, instead of eq (8), the corresponding 


relationship is 


(12) 


The remainder of the treatment for this case is ex- 
actly the same as that already described. 


5. Discussion 


The effect of an error in the dead space upon the 
determination of the volume of gas adsorbed be- 
comes especially serious whenever the amount of 
adsorbed gas is small compared with the unadsorbed 
gas. Indeed, the uncertainty in the measurement 
of the dead space can overshadow the unreliability 
of all other factors which contribute to the over-all 
error in V,,. This error can be decreased appreci- 
ably by a different choice of experimental conditions 
designed to increase the adsorption. The use of 
krypton as the adsorbate at liquid nitrogen temper- 
atures has become increasingly popular in this con- 
nection especially where the surface areas to be 
measured are extremely small (i. e. of the order of 
100 em*/g). 

The technique described in this paper is not pro- 
posed as an alternative to the use of krypton, but 
rather as a supplement. It can be applied to iso- 
therms obtained with krypton as well as with 
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nitrogen or any other gas for that matter, as long 
as the isotherm equations are applicable. For many 
practical reasons‘ it is desirable to use nitrogen in 
weference to krypton where at all possible. The 
imiting criterion should be whether sufficient pre- 
cision is obtainable using nitrogen based upon the 
scatter of the adsorption points alone. 

It is a common experience in many laboratories 
to repeat adsorption measurements because many of 
the points fail to conform to the BET isotherm 
equation. Small errors in the dead space could 
account for such failures, especially when the volume 
adsorbed is small. For such cases the technique of 
this paper is of considerable value. 

There may be reasons for discrepancies existing 
between measurements of dead space and their cor- 
rect values other than the usual limitations in the 
precision of the experimental observation. Helium 
may itself be adsorbed upon certain materials * at 
the temperature of the dead space calibration. Also, 
probably more commonly than realized, adsorbents 
mav have an appreciable volume within pores larger 
than that of the helium molecule, but smaller than 
nitrogen. This would give rise to a molecular-sieve 
phenomenon that might alter the dead space obtained 
from a helium calibration. 

The most reliable estimate of surface area is ob- 
tained by taking the arithmetical mean of several in- 
dependent determinations. In this way, chance 
variations even among samples can be averaged out. 
Such a procedure for as few as two determinations 
can seldom be justified economically. The only 
practical alternative is to make a single determina- 
tion. Hence, it is especially important to select a 
representative sample and to employ an analytical 


1. Highly sensitive differential pressure measuring techniques are required 


since all the adsorption points must fall within the pressure range of 0.2 to 0.6 mm 
where the BET isotherm is valid 

2. Long periods of time are required before equilibrium is attained at eact 
ulsorption pressure 

3. Thermal transpiration (thermo-molecular diffusion) effects become so pr 


nounced that special calibration is required to correct for this phenomenon 

* A polymer carbon described by F. H. Winslow, W. O. Baker, N. R. Pape, 
ind W. Matreyek, J. Polymer Sci. 16, 101 (1955) adsorbed 3 cm of helium per g 
even at 30°C. The surface area of this sample (nearly 600 m was determined 
from its nitrogen adsorption at 90° K. The absolute error of the dead space was 
quite large because of the appreciable adsorption of helium at that temperature 
The relative effect of this error upon the surface area was neglibible, however 
because of the very large adsorption of nitrogen as compare d with the unadsorbed 
gas in equilibrium with it 


technique which leads to results of the highest pre 
cision that the data are capable of vielding = 
estimate of this precision is also very desirable Tr 
is believed that the present work accomplishes these 
ends. 

Obviously, the calculations are somewhat 


, - - more 
time consuming with the proposed technique than 
with the conventional method of determining dead 
space. The time saved by elimination of the dead 


space calibration partially compensates for the 
additional time required for the calculations. 

The more experimental points taken within the 
applicable pressure range of the isotherm equation 
the better one is able to distinguish between experi- 
mental error and incorrect dead space. However. a 
large number of points is time consuming, and jp. 
creases the possibilits of introducing errors due to 
undetected changes in experimental conditions 
While it is bevond the scope of this paper to deter- 
mine rigorously the minimum number of points re- 
quired, it is felt from experience that about eight is a 
practical compromise. ' 

The helium calibration emploved in the conven- 
tional method for surface area determination some- 
times poses a problem not universally encountered 
The limited availability of helium in many areas 
could act as an effective deterrent for its use. partie- 
ularly where surface area measurements are of a 
rountine nature. For these localities a technique 
such as that described in this paper should be 
especially attractive. 
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Sample Calculations of Gamma-Ray Penetration 
Into Shelters: 


Contributions of Sky Shine and Roof Contamination ' 


Martin J. Berger and James C. Lamkin 


\n approximate method is presented for calculating the penetration of gamma radia- 


tion in shelters 


volts, are given for the following problems: 


Sample calculations, for an assumed source energy of 1 million electron 


Dose rate inside houses and underground 


shelters whose roofs are covered with radioactive fallout, and (2) dose rate in open holes 


due to reflected radiation (sky shine 


from fallout contamination on the surrounding ground. 


\ detailed examination is made of the dependence of the dose rate in a shelter on the shape 


of the shelter, and on the position of the detector within the shelter. 


of the calculations is + 30 percent. 


1. Introduction 


The shielding against gamma radiation provided 
by shelters, houses, and similar structures is of con- 
siderable practical interest. The experimental data 
on this problem are usually obtained under complex 
drecumstances, so that they are difficult to interpret. 
It is desirable, therefore, to obtain some theoretical 
understanding of the shielding problem. Calcula- 
tions with the required accuracy (say 30 to 50%) 
appear feasible, and are perhaps easier than experi- 
mental determinations. 

The well-developed semianalytic theory of the 
penetration and diffusion of gamma_ radiation 
1,3]? treats only infinite homogeneous media, and 
is thus not directly applicable to the problem under 
consideration which involves more complicated 
boundary conditions. Yet the theory provides the 
foundations on which the calculations can be carried 
out. It would be possible, but rather laborious, 
to take the boundary conditions into account by 
means of random sampling (Monte Carlo). An- 
other approach appears more promising—the adapta- 
tion of infinite-medium results to realistic conditions 
through a suitable schematization. 

It is the purpose of this paper to explore one such 
schematization that can be applied to calculations 
if shi¢ lding against fallout radiation. In this schemat- 
vation, the radiation dose inside a shelter is com- 
puted as a suitable weighted integral over the angular 
distribution of the radiation dose in an_ infinite 
medium. The present paper contains only pilot 
calculations. More extensive and systematic cal- 
culations are now in progress as part of a National 
Bureau of Standards program of shelter evaluation, 
carried out under the auspices of the Federal Civil 
Defense Administration. 
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2. Schematization 


The nature of the schematization can be most 
easily explained by applications to specific problems. 
Generalization to other situations is then straight- 
forward. Figure 1 illustrates the schematization 
to be employed for three typical situations that will 
form the basis for our discussion. 


2.1. Underground Shelter With Roof Covered by 
Fallout 


The shelter roof and the surrounding ground are 
assumed to be covered with radioactive material 
that constitutes a plane isotropic source of gamma 
radiation of constant strength per unit area, which 
emits gamma rays uniformily in all directions. 

Simplification 1. Different materials of low aver- 
age atomic number, such as air, concrete, earth, or 
water, have nearly the same gamma-ray mass 
attenuation coefficients. Provided that all distances 
from the source are expressed in units of mass per 
unit area, it is therefore a good approximation to 
replace the combination air-concrete (or earth) by a 
concrete medium that is infinite and homogeneous, 
except that it contains an air cavity (the inside of 
the shelter). 

Simplification 2. In order to reach a detector 
somewhere inside the shelter, radiation must pene- 
trate varving amounts of concrete or earth, depend- 
ing on its path. Radiation coming through the 
ceiling has the least amount of material to traverse 
and is therefore least attenuated. There are two 
types of paths on which photons could reach the 
detector through the floor. They can detour the 
shelter to come up from below, but in this case the 
attenuation is so strong that the contribution of the 
detoured photons is negligible. There is also the more 
important possibility for photons to enter the shel- 
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ter through the ceiling. to penetrate into” the 
floor (or side walls) and then to be backscattered. 


Some radiation will also penetrate to the detector 
through the side walls of the shelter. But because 
of the small size of a mean free path in concrete or 
earth (a few inches) only a small side-wall region 
near the ceiling can make a nonnegligible contribu- 
tion, resulting in an effective ceiling area somewhat 
larger than the true ceiling area. Both the back- 
scattering correction and the side-wall penetration 


correction are small. They are briefly discussed in 
sections 6.1 and 6.2 but are omitted from = the 
initial schematization. It will be a good first ap- 


proximation to consider only radiation entering 
through the ceiling. The mean free path in air is 
of the order of several hundred feet, and is very 
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large compared to the shelter dimensions. Th 
radiation thus travels in a straight line inside the 
shelter and ean reach a detector only if its direction 
of motion lies within the solid angle subtended by 
the ceiling with respect to the detector. This solid 
angle is indicated in figure 1 by the label ‘allowed 
directions.”” The notion of “allowed directions” will 
be formalized later through the introduction of an 
“angular response function.”’ 
Summa yoof the Calculation 
sists of the following steps: (1) Replacement of the 
air-shelter complex by an infinite concrete medium; 
(2) computation of the amount and angular distribu- 
tion of the radiation in such a medium, at a distance 
from an infinite plane isotropic source equal to the 
thickness of the shelter roof: (3 computation of the 


The calculation con- 
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function for a specific shelter- 


ygular response | . 
letector configuration; (4) computation of the total 


mount of radiation received by the detector, by an 
ptegral over the angular distribution weighted by 
he angular response function; and (5) corrections to 


count for additional radiation reaching the de- 
wetor through the side walls, and by the reflection 
fom interior surfaces of the shelter. 


2.2. House With Roof Covered by Fallout 


4s indicated in figure 1, the schematization for 
this problem is essentially the same as for the 
previous problem. A number of points require com- 
ment, however. 

1) The contamination now ts confined to the roof 
«that the plane-isotropic source is finite. Yet the 
situation is practically equivalent to that of an 
underground shelter underneath an infinite con- 
taminated plane. Note that in the underground 
shelter, fallout from regions other than that directly 
above the roof does not contribute appreciably to 
the dose, because of the large amount of intervening 
earth. 

The equivalence is, of course, only approximate. 
To illustrate the order of magnitude of the error in- 
volved, one may consider the following point. In 
the schematization, approximately 10 to 15 percent 
of the energy emitted by the source Is reflected back 
toward the roof. In the actual situation, with air 
above the roof, this reflected radiation would spread 
laterally, and most of it would come down on the 
ground surrounding the roof. The reflected radiation 
assumed incorrectly to be incident on the roof is low 
in energy and diffuse in direction, so that its pene- 
trating power is much less than that of the source 
radiation. Hence, the resulting dose-overestimate 
will only be a few percent. 

2) There is a somewhat better chance than in the 
ease of an underground shelter, that radiation origi- 
nating on the roof will reach the detector through the 
side walls. But this radiation consists mainly of 
photons which after repeated all scattering are Incl- 
dent on the side walls with low energies and therefore 
are not very penetrating 

3) If the ground surrounding the house is also 
covered with fallout, the penetration of radiation 
through the side walls becomes very important. The 
numerical details of this case are not worked out in 
the present paper, but the same schematization could 
be applied as in the case of a roof contamination, 
through the introduction of appropriate allowed 
directions and angular response functions defined 
with respect to the side walls, or perhaps the windows. 


2.3. Open Hole Surrounded by Fallout 


Again the schematization Is essentially the same 
isin 2.1. The detector is now in a position where it 
cannot ‘see’? the source, i. e., only seattered radiation 


ean reach it. This seattered radiation will consist 
mainly of photons that have gone up into the air and 
are reflected into the hole by scattering. This radia- 
won COM pe nent is customarily called shy sh ine. The 


radiation distributions that must now be used in steps 
) and (4) of the schematized calculation are those 
pertaining to the scattered flux at the source plane. 


(2 


3. Mathematical Formulation 


An isotropic source of gramma radiation is as- 
sumed to be located in the plane z=0. Let Jp(2, @) 
sin 6dé denote the dose delivered to a detector in an 
infinite homogeneous medium at a distance 2 from 
the source plane, by photons incident on the detector 
at angles between @ and 6+dé@ (with respect to the 
--axis). Dose means the quantity technically called 
“absorbed air dose”, i. e., energy dissipated in air, 
measured in rads (100 ergs/g) [2]. It will be con- 
venient to express the dose distribution in the form 


1 »(2,0) f(2g(2,9), (1) 


#) is an angular distribution normalized to 


| g(z,0) sin 6dé=1. (2) 


where g(2, 
unity: 


) 


The factor f(z) will represent the attenuation of 
radiation by the roof of thickness z in the problems 
of the house or underground shelter. For the open 
hole, one must replace f(z) and g(z,@) by quantities 


f,(2) and g,(2,@) pertaining to scattered radiation. 


Let the shelter-detector configuration be such that 
only a fraction ¥(@) of the photons with direction 6 
can reach the detector. This function is called the 
angular response function and is further discussed 
in section 5. The dose ))(z) received by the partially 
shielded detector in the shelter is expressible as the 
product of three factors: 

D(z)=Cf(z) G(z), 


where 


G (2) q (2,0) ¥(0) sin 6dé, 


and C is a correction factor taking into account the 
contribution of radiation penetrating through the 
side walls or reflected from the interior surfaces. 
The determination of the geometrical protection 
factor G is the main objective of this paper. 


4. Calculation of the Radiation Field 


The standard solution for an infinite homogeneous 
medium is provided by the moment method [1]. 
Tabulations on this method are available 
from which one can extract information about the 
dose f(z) from a plane isotropic source [3].° The 
calculation of the angular distribution g(z,@) has been 
described in {4}. 

The numerical computations of this paper all per- 
tain to a source emitting 1-Mev photons. [his 
source energy is a little higher than the average 
energy of fallout radiation, but close to the energy 


based 


As part of the work sponsored by FCDA, further tabulations, with applica- 
sto shelter « ng prepared it NBS, 


iluation, are now be 
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Fiat RE 2. dose from a 


1-Mev plane rsolropu source im an infinite medium 


Angular distribution of absorbed air 


rhe calculation is for a water medium, but is in ver 
applicable to concrete 


y good approximation also 


4) 9(2,0), at a distance wor=4 from the source plane 
C) @o(2,0), at the source plane yor =0 Phis distribution represents only the 
contribution of scattered photons 
rhe curve (b) represents an isotropic distributior 
of Co”-radiation, which is most suitable for con- 
‘ , vane ; : 
trolled experimentation. The choice was moti- 


vated by the ready availability of results for g(z,@). 
In any case, calculations with a 1-Mev source will 


be sufficient for the illustrative purposes of this 
paper. 
Let f(z) be the dose in an infinite concrete medium 


at a distance z from a plane isotropic source. The 
results of a calculation by the moment method can 
be summarized by the approximate formula 


with 


fo! Z kk, (flyz), ja) 


representing the contribution of unscattered radia- 
tion, and with 

t (2 1.04ke —°-?!*H0? (5b) 
representing the contribution of scattered radiation. 
The range of validity of (5b) is 0<p 28, F; is the 


exponential integral [5], and wo=0.0635 em*/g is the 
mass attenuation coefficient for gamma radiation in 
concrete at the source energy (1 Mev). The constant 
i determines the source normalization. If the source 
plane is covered with 1 curie of gamma emitter per 
square foot, and the dose rate is expressed in rads 
per hour, the normalization constant has the value 


I i fe a (5e) 
If the fallout is spread over a rough surface, the 
schematization of the source as a perfect geometric 
plane is not altogether suitable. A possible alterng- 
tive schematization consists in assuming that the 
radioactive material is mixed homogeneously with a 
thin top laver of the ground of thickness ft. In this 
case, the dose function f(z) must be replaced by 


| jie oe 
4) J (2") dz’. 


Three angular distributions are shown in figure 2: 
(a) g(z, 8) ata distance poz=4 (0.897 lb of concrete 
per square inch); (b) a completely isotropic distri- 
bution; and (c) the distribution ,(2,6) for seattered 
radiation only, at the source plane yoz=0. The 
distribution (a) peaks at 60° and is for use in the 
problems treated in sections 2.1 and 2.2. The dis- 
tribution (¢) peaks at @=90° and is for use in the 
problem in section 2.3. Both of these distributions 
pertain to an infinite plane-isotropic source and 
were calculated according to the method of reference 
[4]. The distribution (b) is flat (when plotted against 
cos 6) and has been included as an intermediate case. 
Together, the three cases Span a wide range of 
possible distributions. 


f* (z,t) (6) 


5. Calculation of the Response Function 


To each angle @ there corresponds a family of unit 
vectors (sinécosg, sin@sing, cosé), indicating possible 
photon directions. If these unit vectors are required 
by the schematization to lie within a specified solid 
angle, then not all values of the azimuth ¢ are 
permitted, but only azimuths in a region Ag <2. The 
response function is defined as the fraction (8) = Ag 27 


5.1. Cylindrical Geometry 


As a simple first example, we consider a shelter in 
the form of an upright evlinder of height HZ and 
radius 2, with the detector located at the center of 
the bottom of the shelter The response 


see fic. > 


FigurgE 3 Calculation of the netion for a cylindrica 


shelte 


response hu 


112 


fyncul 


We 
widtl 
loss { 
ton | 
of th 
shelt 
divid 


the s 


and | 


Tl 
illus' 
of th 


The 


pars 


K 
to 
det 
tha 
tan 
of 1 
hei 
bot 
COO 
hov 
ver 
pas 
per 





junction is in this case particularly simple: 

















(5¢e) 
H ‘ " 
©, the 0 <8<0,=cos i Re Ag=2r, ¥(0)=1 
Metric V dines (7) , oe So \ 7 | 
terna- \ 
t th <O<m , Ap=0, ¥(0)=0. \ | 
- le . \ TOP VIEW 
With a ‘ - \ | 
n thi 5.2. Parallelepiped Geometry ng ' \ 
7 . | mau | _ 
We consider a box-shaped shelter of length JZ, ee we —P 
width W, and height /7. It can be assumed without —p a ae 
(e) loss of generality that L>W. The response func- 
“ion depends on solid angles, and is thus a function ; 
of the shape but not of the absolute size of the Ff 
ure 2. | chelter. Therefore, we eliminate one parameter by 
crete | dividing all dimensions of the shelter by ZL, so that ; SIDE VIEW 
listri- the st aled length Is unIty, the sealed width is a 
ttered o 
The W/L=e<1, (8) | | . 
n the 
‘ dis- and the sealed height is 
i the ' all 
‘tions H/L=h. (9) 
and "7 ; , : 7 . Figure 4. Calculation of the response function for a box- 
rence The construction of the response function is shaped shelter. 
ainst illustrated in figure t for a detector in a bottom corner 
case, of the shelter. 
re of 
(Angular range Response function, ¥(6,€,h) 
on 
j 
unit O<0S%=cos™" ; L 
2+ 6 
sible . 
red Pe ee ee sin-1 (__* —) (10) 
solid \ pe | or h tan @ 
are 
The } - € 1 l 
ri < f «Qs! sina! o—= 2Oa~! 
Qn pian vVh?+€ | = (j tan ) as” (G tan a) 
“f <f 0 
r in . 
and The response function (10), evaluated with the | 
rof | parameters e=0.5 and h=0.5, is shown in figure 5. 25 
mse 
5.3. Combination Rule » 
Expression (10) can readily be generalized sO as 
to provide the response function for an arbitrary tes 
detector location. In the first place, it is to be noted 
that—within our schematization—the vertical dis- 
tance of the detector from the shelter ceiling (or top 10 
of the open hole) must be assumed as the effective 
height Il. Now consider a detector located at the 05 
bottom of the shelter, in a horizontal position ? with 
coordinates (@,Be), as illustrated in figure 6. We P 
now imagine that the shelter is partitioned by two 4 5 6 7 8 9 LO 
vertical walls (perpendicular to each other), which m9 
. Figure 5. Response function for a detector at a bottom corner 
ass , “a are ver ; , ; , 
cal +t | through the point r, und ure ve ry thin but oj a box-shaped shelter with scaled width (horizontal eccen- 


perfect, radiation shields. The detector can then tricity) e=0.5 and height h=0.5. 
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. As a rule, only the first few harmonies (=0. 1 
<7) are at one’s disposal, but they suffice to specify 
i ‘ y » j . “% 
| g(2,6). When (12) and (15) are substituted into (4) 
| 4. it follows from the orthogonality of the Legendre 
polynomials that 
EE - a _ 4 ’ : 
2/+-1 2n+1 
’ — a J y y . a J 
fae epee | 2 2 
5 ‘ a 
Ficure 6 Diagram illustrating the application of the com- Ss > > , 
bination rule by which response functions for off-corne: l (COS 6 P,,(cos @)sin ade (14) P 
positions are calculated _ 
L 9] l 
, 255 U9 
be thought of as being located at the corner of each i=0 < 
of four adjoining box-shaped shelters with the 
dimensions listed below. Thus, one can limit oneself to the calculation of the 
first few response function expansion coefficients 
- j a; | ¥ (0)P (cos Asin 6d@,/7=0,1, ..., LZ. (15 _ 
Jo 
min (a, Be } - — — — 
I With these coefficients, one is in a position to evaluate , FOUR 
max (a, Be max (a, Be m . = : : 
G(z) for any angular distribution g(z,@) through 
’ min (l—ae, Be simple superposition according to (14). <i 
’ max (l—a, Be max (l—e, Be) Numerical evaluations of G have been carried out | mre 
for detectors in various positions in box-shaped |“ 
3 min (a, ¢— pe shelters, on the basis of the angular distribution 
MAX (a, e— Pe max (a, o— functions shown in figure 2. These results apply as Fis 
' min (1 am, € De / follows: with 
ma ] a, € Be MisaN | ae Be) tains 
' 

— 
0<a<1, 0<s<!1 Angular Infinite medium and 
min (2, y)=smaller of « and y. Problem distribu- nee conti 
max (r, y)=larger of r and y. tion ee the | 

cures toot- . 
tect 

ris i 

The response function Wp at point P is Underground shelter’ or A f(z) =0.983 funet 
house with roof covered smal 
‘ by fallout, roof thick- rd | 
. ordae 
(A S ‘vu A.€ } y 11) ness poz 1(0.896 Ib of 
vr —¥ “ concrete per square cent 
inch) forw 
: : , , arge 
5.4. Legendre Polynomial Expansion Isotropic comparison dis- ty _ 
tribution ~~ 
In order to obtain the geometric protection factor : —_ 7 
Y/ = a . Re ¢ ° . Skv shine (seattered radia- Cc / 2 ded. lr" 
G(z), it is necessary to evaluate the integral in (4) for tion) into open hole 
. ° . . . ra ° ’ ILO res] 
each distribution g(z,@) of interest. When many dis- if 
tributions must be included, this procedure tends to ; re 
become laborious. It is then convenient to perform ; _— 
a set of standard integrations based on an expansion The results are presented as functions of the wvlit 
of the response function into a series of Legendre following dimensionless parameters: ie ~ 
polynomials: we 
: , plot 
- : a width ' 
wi 2l+1 P(e . Horizontal eccentricity € (16 the 
vy (@) = qa ait { (cos A). 12) leneth | resy 
U = 
shel 
The theory of radiation diffusion provides as the Pas ont yroof-area ye 17) 4 Prot 
: . oe Vertical eccentricity f. of t! 
most direct result, a similar Legendre polynomial height h 
expansion for the angular distribution: dem 
ren . . P ' hor 
I Phe range of parameters considered includes O0<e<l 
(> § wi2l+l Polo a and 0.2<7r<20. This appears sufficient for all 
q(z,6) dy cos @). L3 : aie f ‘ ' 
. . ») . . > . . 
= 2 practical applications 
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tri ily, for fired hor zontal ecce ntricily € r. 


ire based on the use of the corresponding 
and hollow circles 
i cylindrical 


The curves labeled 1), (b), and (¢ 
Jar distributions in figure 2 he triangles, solid circles, 


wpresent results obtained with the response function (7) for 


Figure 7 contains plots of @ versus r+, for a shelter 
with a square cross section (e=1). Figure 8 con- 
tains plots of G@ versus e, for different fixed values of 
; Results are shown for the detector in a corner, 
and in a center position (a=8=43). The figures 
confirm quantitatively what one would expect on 
the basis of a commonsense estimate. The pro- 
tection provided by the shelter increases rapidly as 
ris increased, but is generally a very slowly varying 
function of ¢« for fixed +r, except when ¢ has a very 
small value. Dose rates in the corner are on the 
order of four times smaller than dose rates in the 
center. The angular distribution (a) with its strong 
forward peak gives rise to G-values 2 to 3 times 
lager than those obtained with the diffuse sky- 
shine distribution (c), whereas there is little difference 
between cases (c) and (b). 

The question may be asked to what extent the 
rsults for the geometrical protection factor are 
affected by the choice of a rectangular shelter cross 
section. To obtain a partial answer, calculations 
were carried out for a detector at the center of a 
evlindrical shelter, with the of the response 
function As can be seen from these values 
plotted in figure 7, the results thus obtained are, for 
the same value of +r, practically identical with the 
results for a detector at the center of a box-shaped 
shelter with e=1. This means that the geometric 
protection factor depends very little on the shape 
of the cross section. Combined with the previously 
demonstrated insensitivity of G with respect to the 
horizontal eccentricity ¢, this has the following 
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tricity, for various constant values of the vertical eccentricity. 


pertain to the corresponding angular dis- 
detector in corner. 


(b), and (ce 
Detector in center; (B), 


The curves labeled (a), 
tributions of figure 2. A 


implications: (1) The geometric protection factor 
at the center of a shelter of any shape is, in good 
approximation, only a function of the vertical 
eccentricity 7; (2) under these circumstances, a 
calculation using the very simple response function 
(7) will be adequate. 

Center and corner values provide upper and lower 
limits of G. To indicate the transition between the 
two extremes, a systematic mapping of G is presented 
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in table 1 for numerous detector positions in the 
horizontal plane, the shelter being box-shaped with 


e=1 and r=2. 


TABLE 1. Geometric protection factor G as a function of the 
position in the horizontal plane, for a detector in a shelter 
with horizontal eccentricily € 1 and ertica eccentricity 


r= 2, calculated for the angular distribution (a) of figure 2 


The position indicators a and 8 are those defined in figure ¢ Note that the 
ransformations a—|1—«, 8-1—8, or (a8, Ba) leave G unchanged 
a 
oO 0. 125 0. 25 0. 375 O35 
8 
0. 0 O17 0. 22 0. 26 0. 27 0). 29 
. 125 22 3! 34 37 38 
25 26 34 1) 13 +4 
375 27 37 13 17 iS 
5 2% 38 14 tS 1) 


Corrections, Accuracy 


6.1. Penetration Through the Side Walls 


This effect has, so far, been investigated only in a 
limited manner to obtain the order of magnitude of 
the necessary correction. To simplify the caleula- 
tions, they were done for a detector at the center of 
the bottom of a cylindrical cavity. For the case of 
an underground shelter (the problem in section 2.1 
the response function (7) could be amended as follows: 


T A2+2’) g(z2+2’.6 
0,<O0<-, W(a)=" d 
at yg q a] 
7 

. < 

5< OS, ¥(0)—0 

where 
>'=H—R cot 6. 


The meaning of (7’) is that the surface elements of 
the side wall contribute to the dose in proportion to 
the intensity and angular distribution of the radi- 
ation at the distance ’ of the surface elements 
from the source plane. 

Sample calculations have been carried out for a 
shelter with roof thickness uoz=4, horizontal ececen- 
tricity e=1, and a depth of 8 ft. (The side-wall cor- 
rection differs from the uncorrected response fune- 
tion in that it depends on the absolute size of the 


116 


shelter.) It was found that for vertical eccentriciti 
‘ss 


r=1, 2, and 4, 
the increase of the dose is 7,3, and 1 percent, respec 
tively. 


6.2. Backscattering 


Some radiation will be scattered tow ard the deter. 
tor by reflection from the side walls, floor, and ceiling 
of the shelter. The order of magnitude of this effec 
can be estimated roughly on the basis of available 
Monte Carlo calculations of backscattering. Tho 
results of reference 6] suggest a dose increase by 
approximately 15 percent. 


6.3. Over-all Correction and Error Estimate 


In the final equation for the dose received by the 
partially shielded detector, , 
Di 2) (7 (i(- 2 

the factor (, incorporating side-wall and backseat- 
tering effects, is estimated to be 1.2. It is further 
estimated that the factors in (3) have the relative 


errors (6fG)/fa=+0.2 and (6C)/¢ +0.2. whieh 
leads to the following error estimate for D: 

6D ‘= ) i) n2 

D VG TN _ 18 


The authors are indebted to L. V. Spencer for 
valuable discussions, and to J. Hubbell for help with 
the calculations. 


7. References 


{1} L. V. Spencer and U. Fano, J. Research NBS 46, 446 
1951) RP2213 
[2] Report of the International Commission on Radiological 


Units and Measurements, NBS Handbook 62 (1956), 


[3] H. Goldstein and J. E. Wilkins, Jr., AEC Report NYO- 
3075 (1954 

[4] M. J. Berger, J App! Phys. 26, 1504 (1955 

15] Table of sine, cosine and exponential integrals, Federal 
Works Ageney (1940 

[6] M. J. Berger and J. A. Doggett, J. Research NBS 56, 89 


(1956) RP2653 


WasHiInecron, May 16, 1957. 


mal © 


Ther 
he un 
younds 
nan 
\lore I 
or pre 
he me 


he all 


I known 


oom | 


} resist 


ip to 
Since 
netall 
elativ 
ted 
eter 
In | 
neasu 


! opel 


The w 
venere 
vdric 
rom. | 
begins 


Wan 


The 


sump! 


| ; 
' Was 


is TEC 
+ prelin 
ater. 
POUT 


' weigh 





ricities 


‘espec- 


deter. 
reiling 
efter 
tilable 
The 


se by 


te 


\ the 


Scat- 
rthes 
lative 
vhich 


nai of Research of the National Bureau of Standards 


Vol. 60, No. 2, February 1958 Research Paper 2828 


Heat Content of Sodium Borohydride and of 
Potassium Borohydride from O° to 400° C 


Thomas B. Douglas and Ann W. Harman 


The heat content of sodium borohydride, NaBHy,, and of potassium borohydride, 


KX BHy, was measured, relative to that at 0 


‘, at 50-degree intervals from O° to 400° C. 


The so-called “drop” method was used with a precision Bunsen ice calorimeter and a silver- 


core furnace. Preliminary 


tests showed that spontaneous decomposition of these com- 
pounds is very slow at temperatures up to 400° C. 


Measurements were made on two samples 


of each compound, those of sodium borohydride having been obtained from different sources. 
Smoothed values of the heat-content function, heat capacity, entropy, and the Gibbs free- 


energy 


funetion were derived from the data. 


These values were increased by 0.5 percent 


to correct for the presence of an estimated | percent of the corresponding metaborate in each 


sample. 


The heat capacity of sodium borohydride increases monotonically with tempera- 


ture in the range investigated; that of potassium borohydride passes through a point of in- 
flection, which suggests that there may be a broad second-order transition in the neighbor- 


hood of room temperature 


1. Introduction 


There has been a wide interest for many years In 
je unusual properties and possible uses of com- 
wounds between boron and hydrogen, and a great 
anv of them have been prepared and_ studied. 
\lore recently, chemical methods have been devised 
x preparing metallic borohydrides, which contain 
ve monovalent group BH,;. The borohydrides of 
ve alkali metals are much more stable than the 


linown hydrides of boron, dissolving in water at 


. for 
with 


146 
gical 
O56), 
YO- 
leral 
6, 89 


oom temperature with but slow decomposition, and 
msisting spontaneous decomposition in the dry state 
to temperatures of several hundred degrees C. 
since X-ray patterns indicate a crystal lattice of 
vetallic and borohydride ions, it appears that the 
lative stability of these compounds can be attrib- 
ted directly to the high coulomb energy char- 
wterizing such typical ionic salts. 

In recent vears a number of investigators have 
neasured some of the fundamental thermodynamic 
roperties of the borohydrides of the alkali metals.' 
The work reported in this paper contributes to this 
gneral program. The heat content of the boro- 
ydrides of sodium and potassium was measured 
rom 0° C up to temperatures where decomposition 
egins, and from these values other important thermo- 
lvnamic properties were derived. 


2. Samples 


were on two 
Each sample 


The heat-content measurements 


amples of NaBH, and two of KBH,. 


'Wwas in the form of fine white crystals, and was used 


preliminary 


s received without further purification, except for 
heating and evacuation as described 
It appears that purification of these com- 
bounds, to a high degree (e. g., to a purity of 99.9 
veight percent), is difficult with present techniques. 


iter. 


a 


References to such work ace given later in this pape 
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This fact may be due to slight solvolysis of the salts 
when in solution, with subsequent decomposition of 
the free acid, for it has been reported that increasing 
the pH of the solution increases the stability of the 
salt. 

One of the samples of NaBHy,, hereafter designated 
as No. 1, was obtained from Metal Hydrides Ine., 
Beverly, Mass., who reported it to have a purity in 
excess of 99 percent. This sample was received in a 
glass-stoppered bottle, and was so kept until used. 
The other sample, designated as No. 2, was obtained 
from another source. Its suppliers had recrystallized 
it twice from liquid ammonia. Their subsequent 
analyses indicated the following percentages of the 
theoretical: Potassium, 99.9 percent; boron, 99.9 
percent; and hydrogen, 99.4 percent. (The hydro- 
gen reported was that obtained by hydrolysis.) 
However, considering the limitations of reliability of 
the analytical methods used, they believed that the 
purity of this sample (No. 2) could not be more 
precisely stated than as “safely above 99 percent”’. 
The sample was received in a sealed glass ampule, 
which was not opened until very shortly before the 
sumple was used. 

Both “samples” of the KBH,, distinguished 
subsequently as No. 1 and No. 2, were really speci- 
mens of the same sample, but specimens on which 
heat-content measurements were made 2 months 
apart. This “stock”? sample was obtained from 
Metal Hydrides Inc., who had found it to contain 
0.1 percent by weight of chlorine and 99.6 percent 
of the theoretical hydrogen (the latter figure resulting 
from gasometric analysis). On this basis they 
believed the sample to have a purity of approximately 
99.5 percent. The sample was kept in a vacuum 
desiccator until very shortly before the heat-contert 
measurements were begun. 

Assuming that at some stage during their purifica- 
tion these borohydrides were nearly 100-percent 
pure, contamination by the atmosphere very likely 
occurred by conversion of some of the borohydride 
to borates. It appears that whatever borate is thus 
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formed would maintain the 1:1 atomic ratio of alkali Alkali Metal, Gravimetrically. After decompogi.| ett 
“ele . ; . : . si-} ” 
metal and boron characterizing the parent boro- | tion of the borohydride with strong acid, the solution, © 2! 
hydride. Such a borate might be a metaborate or, | was evaporated with sulfuric acid. After evapora. compe 
before an opportunity for dehydration, more likely | tion of the boric acid and ignition of the residue | with ! 


an orthoborate (such as KH,BO,-H,O). All such | the alkali metal was weighed as sulfate. jater | 
borates other than the metaborate are known to de- The foregoing procedures were tested on potassium! the 


hydrate to form the corresponding metaborate at the | borohydride using an extra sample on which no heat.| meas 
higher temperatures of the present work, and at | content measurements were made. The results are rectiOl 
lower temperatures they would be expected to | given in table 1. Each value in the third colons apple 
possess heat capacities which are approximately | was calculated with the assumption that the partic.| 
additive combinations of the heat capacities of the ular element determined was present only in the rae 3. Pi 
constituent metaborate and chemically bound water. | of KBH,. Of these values, the last three. resulting 
For these reasons it seems simpler to interpret the | from the analyses for boron and potassium, agre: 
analyses of the samples as though the somewhat | closely among themselves but are all appreciabh It 
impure sample consisted wholly of borohydride, | higher than the first value, which resulted from the! 


metaborate, and water. Letting \V/ represent either | amount of hydrogen evolved. Averaging the lag | ioe 
Na or K, decomposition by water or oxygen could | three values and assuming that the sample contained valid 
then be formulated by the following reactions: only KBH,, KBO,, and H,O, it is calculated that the| obvio 
weight percentages of these components were re. | purity 
MBH,+ 2H,0= MBO, + 4H), (1) | spectively 98.8, 0.8, and 0.4. | in the 
durin 
MBH, + 20,= BO, + 2H,0. (2) | Tape 1 Test of analytical methods on an extra sample 6 samp 
potassium borohydride SOI 
(In time, the water formed according to eq (2) might hydri 
react with additional borohydride according to No. of de- Meanresult) Average | tera 
eq ( 1 ).) Element determined — as on —— that 
: : , ‘ « from mean 7S 
Several months after the completion of the thermal | purit 
measurements reported in this paper, specimens of Hydrogen (evolution) 3 98. 79 0.10 decot 
the unused portions of the stock samples were dn heey 3 90. 38 4 | evolu 
analyzed in the Chemistry Division of the Bureau. Potassium (gravimetrically ' 99. 30 Uo) —_ 
The analytical procedures were as follows. ; borol 
Hydrogen Evolved. A sample of 0.5 to 1.0 ¢ was melt 
transferred to a gas-evolution flask and dissolved in The methods for determining hydrogen evolved | °?! 
50 ml of water. Twenty milliliters of 1:1 hydrochloric | and boron present were extended to unused portions | the Ii 
acid was then slowly added dropwise. The hydrogen | of the samples on which heat-content measurements rh 
gas generated was swept by nitrogen through a con- | had been made, duplicate determinations being made weigh 
denser and an absorber containing Asearite (sodium in each case. The mean results, calculated as per- these 
hydroxide on asbestos) and anhydrous magnesium | centages of borohydride, metaborate, and water, are | ear 
perchlorate. The dry hydrogen was then burned in a | given in table 2. The KBH, stock sample analyzed | troll 
combustion tube packed with copper oxide heated | was the source of both heat-content samples No. 1/ ey 
OSL. 


to 700° C, the water formed being absorbed in a | and No. 2 of this substance, and it will be noted that 
weighed absorption tube containing anhydrone. | three sets of results are given for this material. The} “0! 


Each such determination required about 4 hr. The | first set was obtained on a specimen protected from | 2 
accuracy of the method was tested by using standard the atmosphere, and the second and third sets were | migh 
melting-point aluminum (NBS standard sample | obtained on two specimens which had deliberately welg 
44d). The hydrogen evolved in four such analyses | been exposed to the laboratory atmosphere for 5 and | ™E 
of aluminum corresponded to 99.80, 99.62, 99.69, and | 20 days respectively, just preceding analysis. } trace 
99.66 percent of aluminum respectively. Spectro- It will immediately be evident from table 2 that “™ 
graphic analysis had shown the aluminum to contain | these analyses, which were carried out after the heat- |‘ 
0.1 percent of impurities, but the content of oxygen, | content measurements, indicate decidedly less boro-) 2" 
carbon, and nitrogen was not known. hvdride content than the analyses that preceded the | '™ i 
Boron. Either an aliquot portion of the solution hvdi 
. . . ; ii 

of borohydride used in the preceding pre weedure Or sy ABLE 4 Analyse s of horoh yd ide sam ple S atte 7 the heat- smal 
a new sample of 0.25 g was used. After the solution seantends Gamuantinmsiiie ai 
had been adjusted in pH to the methyl-red endpoint, TT 
it was treated with mannitol and then titrated with seis taal arn bic 
0.1N NaOH to the phenolphthalein endpoint. The ener sabi, Tousen h re 
titration was continued until further addition of Weleht Weight © | Weight % me 
mannitol no longer destroyed the pink color. NaBH, sample No. 1 92.6, 3 Se | aa 
Alkali Metal, Volumetrically. A 0.25-g¢ sample was eet, aeemneeed. “ 98. 38 15s 06 ~ 
dissolved in water, treated with a few drops of BH. copeeed 30 days ap _— " Sof t 
methyl-red indicator, and titrated with 0.1.N HCI. i : aon 
Fig 
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beat-content measurements and which were referred 
“ at the beginning of this section. | The slow de- 
omposition of potassium borohydride in contact 
with moist air also is evident. An attempt is made 
ater (section 4) to estimate the chemical composition 
of the samples at the times when the heat-content 


OMposi-| 
Solution 
‘Vapora- 
residue | 


tassium| , 
oe measurements were made, so that appropriate cor- 
ies a -ections to the observed heat-content values may be 
S are; ° ° 
columy appliec ° 
partic. —_ . -): 
he form 3. Preliminary Tests of Chemical Stability 
‘sulting at High Temperatures 
agree 
eclabh; It was desirable to extend the heat-content meas- 


oM the} yrements to temperatures as high as possible without 
he last! jyeurring sufficient chemical decomposition to  in- 
itained| validate the thermal values. In such cases it is 
hat the! obviously important to know not only the initial 
purity of the samples but also any significant change 
|» the degree of purity resulting from decomposition 
during the measurements through exposure of the 
samples to high temperatures, 
Some qualitative information about the 
‘)vdrides of sodium and potassium is available in the 
literature. It has been reported [1],? for example, 
that when potassium borohydride of 98.8-percent 
| purity was heated in a fused silica tube in vacuum, 
decomposition commenced at about 500° C with the 
evolution of hydrogen and the simultaneous distilla- 
tion of potassium, Others have reported that sodium 
; borohydride decomposed only very slowly below its 
} melting point (about 500° C), with formation of free 
boron and free sodium, but that at the melting point 


Te re. 


mople 0 
boro- 


erage 
lation 
1 mean 


. 10 
07 
iL) 


09 


‘olved nal 
rtions | the liquid evolved hydrogen rapidly. 
nents, The present authors tested small specimens (each 


made| Weighing approximately 1 g) of their samples of 
; per-| these two compounds in open containers in a furnace 
r are} constantly being swept out with pure helium, con- 
lyzed | trolling temperature and time of exposure and ob- 


No, 1 / serving changes in weight and appearance after each 
that | test. Sometimes there were small gains in weight, 
The| amounting to a few hundredths of 1 percent, which 


from | naturally would mask any small losses in weight that 
were | might have occurred simultaneously. The gains in 
ately | Weight that occurred during heating at the higher 
sand| temperatures were probably due to absorption of 

‘traces of oxygen or water vapor from the furnace 
Sometimes a loss in weight that oe- 


that atmosphere. 

veat- Curred only on heating to a temperature of several 

yoro- hundred degrees C was followed by a comparable gain 

| the} weight (up to 0.1%) during standing overnight in 
a desiccator. This may be explained if owing to 
uvdrolysis the dried specimens actually contained 

heal-' small amounts of the corresponding metal hydroxides, 


which are hygroscopic. 
| The weight changes mentioned above will next be 
| described quantitatively in the case of first, sodium 
| borohydride, then potassium borohydride. —Im- 
mediately following the discussion of the behavior of 
each substance in the tests of a small specimen, will 
be given a description of the preliminary conditioning 
| of the heat-content specimens of that substance, 


}; ———_— 


Figures in brackets the literature references at the end of this paper 


indicat 
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| with the weight changes then observed. 

The test specimen of sodium borohydride, sample 
No. 1, lost 0.5 percent of its weight in 35 min at 
350° C; in a subsequent 35 min at 400° C the same 
specimen showed no further significant change in 
weight, but became somewhat grav. In view of this 
behavior, the specimens whose heat contents were to 
be measured were first subjected to preliminary 
heating to establish approximately constant weight. 
Sample No. 1 lost 0.1 percent in weight during several 
hours of drying in vacuum at 80° C, no more during 
several hours while being similarly treated at 130° C, 
and an additional 0.2 percent during 30 min of drying 
at 200° and 50 min at 350° C. Sample No. 2, when 
first heated for 30 min at 300° C, lost 0.3 percent in 
weight, but a repetition of this heating produced no 
further sensible weight change. 

The preliminary tests on a small specimen of the 
stock sample of potassium borohydride, with inter- 
mediate cooling and weighing, showed the following 
consecutive results: The only sensible weight change 
was a 0.2-percent loss during the initial heating, 
which was for 40 min at 300° C; there was still no 
change in appearance of the specimen after 35 min 
at 350° C; after 50 min at 400° C, some sintering was 
noted; then 50 min at 450° C led to a dirty-gray 
appearance; and a subsequent heating for 10 min at 
500° C produced somewhat further darkening. The 
only preliminary treatment of the specimens whose 
heat contents were to be measured was the heating 
of each in vacuum at 85° C for several hours, with a 
weight loss of less than 0.05 percent. During the 
subsequent heat-content measurements, each speci- 
men gradually underwent further losses in weight 
totaling about 0.1 percent. 

It may be pointed out that these specimens of the 
two borohydrides initially underwent appreciable 
weight losses, but no comparable losses in subsequent 
heatings even at higher temperatures. Barring con- 
ceivable catalytic behavior of volatile contaminants 
initially removed, it would thus seem that the initial 
loss was not through any sensible, spontaneous 
decomposition of the nominal compound, or else 
comparable or even more rapid weight losses would 
have continued to occur. The losses found are 
assumed to have occurred through the removal of 

| water. Whether this loss of water was more likely 
by simple evaporation or by chemical reaction with 
| borohydride is discussed in section 4. 


4. Correction for Impurity 


The immediate aim of the present work is to obtain 
the best values of the relative heat content of the 
pure borohydrides of sodium and potassium over the 
temperature range covered by the thermal measure- 
ments. As the precision of these measurements was 
several times the errors that would be caused by 
neglecting the presence of impurities in the samples 
measured, these errors will be estimated and 
corrected for. 

It is convenient to calculate first, the changes in 
heat content per unit mass of sample as given by the 
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thermal data, and then to estimate the corrections | heat content of the drying are small and can | 
; “ 


corresponding to replacement of the impurities by | easily allowed for. It is conceivable, however that 
an equal mass of the nominal substance. This cor- | heating caused the reaction of eq (1) to oceur and| 


rection for each impurity is obviously proportional | that the mass loss was due largely to escape of The 
to the weight percentage of the impurity and, | hydrogen, in which case the proper corrections se, used I 
approximately, to the difference between the average | the heat-content values would be many times as| | 
specific heat (heat capacity per unit mass), of the | great. However, the second alternative explanation | a pape 
impurity and that of the nominal substance. The | while it cannot be proved to be entirely false. does| Bl 
results of the chemical analyses given in section 2, | not fit the observed facts. It would require that the ie 
and the interpretation of the mass changes sum- | sample contain, before being heated, about five | (oF a 
marized in section 3, necessitate some arbitrariness | times as much water as the decrease in mass op the ki 
in deducing the compositions of the samples at the | heating; also, the large amount of borate formed | '° 


yorts | 





times when their heat contents were measured. | would lower the specific heat of the sample by ] per- | . = 
Fortunately, however, fairly adequate values are | cent or more. Excepting the final analyses of sae" 
available for the specific heats of the minor com- | sodium borohydride sample No. 1, discussed below by ‘ 
ponents (assumed to be water and NaBO, or KBO,, | the water contents, roughly deduced from the chemi- = 
as discussed in section 2). cal analyses, approximate these mass losses 01 port 
The specific heats of the crystalline mefaborates | heating. Furthermore, as the data in the next see. we 
. . . . . . . . . hl ] 
were calculated from equations given by Kelley [2]. | tion will show, a repetition of heat-content measure. | '"* 
These equations,® giving the heat capacity C, in cal | ments on potassium borohydride sample No. 1 was |e” 
mole deg K~' at K. are: carried out at 150° and 200 C: these later values | 
. na 
were lower than before the loss of 0.1 percent i! > 
- ; ‘ 2 ’ . Ny os fun 
Na BC ae : 12.9+ 0.0122 7 (289°-1239°K). (3) weight fat 200 to 300° ©) by O05 lo ().2 percent . 1) 
amounts which correspond approximately to 4 ; ) a 
’ , p _ a , 3] » pens ' , ater ws . i desir 
K BO, :C, = 12.6 +.0.0126 T (288°-1220°K (4 imple escape of water vapor ‘alibi 
The analy ses ol the samples after the heat- call 
Also using the results on the borohydride samples content: measurements all show less borohydride | @™ 
suis = S ; \ * Se , . anal vsesc fore (car D) s : <a ype, 
reported in this paper, it can readily be shown that — pict 3 Aca pr 6 ae Poh 
. . S i ys , ; five Cs i sO : a 
the heat content of the sample would be increased 7; , | = : — honed 
by verv nearly 0.5 percent for each 1 percent of compositions of the samples when their heat contents | P™ 
metaborate. in the Bs agg ee igen Pe by the | were measured. In this way it may be reasonably wom 
same mass of the fms. onding saadviedvidie "Thi assumed that sodium borohydride sample No. 2 and { a a 
is true for nv tempe ten se in the present range 0 each of the samples of potassium borohydride should |" 
=] " ? " i) ’ " “ePSe ' oe ‘ ° ‘ " . rm 
to 400 und in the Cust of either the odium Ol th be considered - ane of 90 weight LF aghgees of he 
— , "2 7 , 7S ‘a I nd « . . ry ‘| 
potassium compound the nominal substance and 1 percent of the corre- | &@& 
OLUSS . ( . ° rive 
; “ssi - ; sponding metaborate. e. 
A value for the effective specific heat of chemically- : \ similar interpolation of the analy yn sodium | "Me 
° . 4 » i ? i t 7 i SUS ¢ sot ; 
bound water near room temperature Was obtained : | “ds the f 
from the calculated losses in heat capacity when a borohydride sample No. 1 would indicate several t 
> CA r i ? SSCS i Ci ac i . . “ as U 
: ‘S oas uch Water as ; : 1g SS 
number of salts are dehydrated [2], the mean value nee oe ee h ter as the tot ul weight loss on did | 
Sci ales eek mate* dent" * Thie vale ‘ heating and enough NaBO, to give heat-content 
‘ng ? eg iis i e cor- : ; ? Pail te p ; Pa ’ snec 
responds to a specific heat for such water not far 2 goede” ee posed grey mre hie 
. . . . . . Ar»O. &- AC i , Sal © NO. 6 2ave i? lower heal- . 
from that of either sodium borohydride or potassium | oe mini 
borohydride. Only the pota ium ompound cor content values, but by only 0.25 percent at 100° and ther 
Be - ass c 0 Cc }- o . e - : m ‘ er 
é' : -* :, 0.15 percent at 200° and 300° C (section 5). It is } 
tained an appreciable amount of water during any belj I cl elec SealaRady aaa er f le No. } | men 
of the heat-content measurements, and the amounts apace r n we — a CS ee of tl 
present are believed to have been so small (not mor is responsible for the finding of so much impurity in not 
Se are Hheve : , en so sm: re . . ‘ , ’ ( 

: its final analysis, and that it too must have been 
than 0.16; )as to make the corrections in heat content - * tain 
for water negligible about 99-percent pure when its heat content was ree 

>. ° . . 
. . . in : easured., 

As described in section 3, the samples on which measured , : : | ven! 
heat-content measurement were mad wel first In accordance with the foregoing conclusions, each nil 
“at-c , asurements re made were firs ; ce 
heated to dry them Each of the imp! of sodium heat content of a borohydride sample plus cones befe 
> . > 4 zac * SS) cs s ° ° e . » 

Ye a oe < observed in this work was increased by an amount 
borohydride thereby decreased in weight by 0.5 per- l yz & ft) | _ | 7 ar | 
cent, and the samples of potassium borohydride by equa | to ye ge vies of the og re due to = = 
less than 0.05 percent (but by an additional 0.16% sample. Any mar yes of mass gate. Pe ——— 7 
when thev reached about 200° C during the heat- ment Was assumed to represent a Oss O water, had ie 
cuntess teumieains Tt thuee lnctes in oles | T treated as though the mass of the specimen hac 

' ‘asurements). ee inesea ig 
represent escape of water only, the effects on the decreased by that amount. But when @ mass i fl 
ie crease occurred, it was assumed to have oceurred } f : 
entirely through conversion, by a trace of water vapor bh 

Although these e« tions ¢ based on cld data, at least the equati f . : . : ws ; ) 
N 1BO: gives \ sues in the range 6 to HK) c in ~pproximate agreeme nt with the } OF OXN gen gas, of borohydride into metaborate by the - 
pom pen oot ease from more accurate results on other sodium oxysalts | reaction of eq (1) or eq (2). or both. An additional Pe 

‘ Insufficient data were found to justify assigning values at higher temperatures, | correction to the measured heat was applied ac- r 
where the amounts of water left in the borohydride samples were usually negli- | |), lj | taine 
gible anyhow ' coraingty. s fA 
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5. Measurement of Heat Content 

The calorimetric apparatus and “drop” method 
d have been described in detail in numerous re- 
| ports and papers from the Bureau, most recently in 
| a paper on the thermal properties of aluminum oxide 

3) Briefly, the method is as follows. | Phe speci- 
| men ina suitable container Is suspended in a furnace 
for a sufficient length of time to ensure Its reaching 
the known furnace temperature. It is then dropped 
into a Bunsen ice calorimeter, which, by absorbing 
a volume of mercury equivalent to the change in vol- 
| ume due to melting of ice, measures the heat evolved 


I us 


similar measure'pent en the empty container at the 
aire furnace temperature provides correction for the 
contribution of the container and for the heat lost dur- 
ing the fall into the calorimeter, and enables a computa- 
tion of the net change in heat content of the speci- 
men alone. Repetition at a sufficient number of 
furnace temperatures establishes the heat content as 
a function of temperature 

The temperature of the furnace, adjusted to some 


desired even value, Was measured by an accurately 


calibrated strain-free platinum-resistance thermome- 
ter and was kept constant to 0.01 deg Cc. A silver 
pipe, 10 in. long and of 1/2-in. wall thiekness, to- 
gether with end silver “guards,” reduced the tem- 
perature gradients to negligible amounts, thereby 
permpitting the sample to assume the same tempera- 
tureas that measured. That the time of the sample 
in the furnace was adequate was established by de- 
termining the relaxation time in a special drop run 
after a time in the furnace deliberately too short to 
give thermal equilibrium. The adequacy of the 
times used was tested also by varving the times in 
the furnace in each set of duplicate measurements so 
as to detect any systematic time-dependence, which 
did not appear. The furnace core, calorimeter, and 
specimen container were filled with pure helium, 
which flowed through during the measurements to 
minimize not only the times required for reaching 
thermal equilibrium but also the oxidation of speci- 
men and container. A evlindrical vessel constructed 
of the alloy SO Ni-20 Cr, closed approximately but 
not gas-tight by a serew lid of nickel, served to con- 
tain the specimen. The inside surface of the con- 
tamer had been thoroughly cleaned of oxide to pre- 
vent oxidation of the borohydride specimens. The 
container was evacuated and filled with helium just 
before each introduction into the furnace, but with 
ar after each run in order to place the weight on a 
known and reproducible buovaney basis.‘ 

The corrected individual measurements of the heat 
content of the cmpty container are recorded in abso- 
ute joules in table 3. In terms of absolute amount 
of heat, the agreement of the duplicate results at each 
furnace temperature is typical of that usually obtain- 
able with this apparatus in the absence of a material 


There was no reason other than convenience for choosing this particular con- 
tainer; preliminary tests had shown that it was not attacked by the specimens 
* All weights were corrected to vacuum 


TABLE 3. 


Individual heat-content measurements on the empty 
container 


Furnact Meas" red Furnace Measured Furnace Measured 
tempera- heat tempera- heat tempera- heat 
ture ture ture 
( ahs ( ahs ( ahs j 
0 } one. 0 250 } L, 331 , 150 f 2,491.5 
| 253. 6 | 1,330 ( 2, 489. 
f 514.3 f 1,614.0 7 f 2,789.4 
1? 
160 \ 513.8 a0 | 1,614.3 ond 1 2,780.0 
- f 780.7 f 1,902.2 
am 1 780.7 on \ 1,902.4 
= f 1,053.4 f 2,192.6 
. | 1.054 _ | 2193.4 
TABLE 4 Individual heat-content measurements on sodium 
borohydride ~ container 
Relative heat content of NaBH, 
Fun Corrected Hy 0% 
race Tl imne Mass in measred 
tem n fun excess of heat, 
per rhea original sampl Mean Observed 
ture, f +container Individ- observed minus 
ual run for set caleu lated, 
eq 5) 
Sample No. 1 (2.5715 ¢@ 
( m n ahs aha ia aba ia ahs ja 
$5 O.8 4,229.2 w4. 9 \ 
su) ‘3 . + 208 O03 8 ; 904 ‘ 0.4 
; 40 546, 2 At TY sma 
wo As 545. | naz if 9 . 
i 30 115.) 33. 7 | oe 
lit -~ , 
Yn 1 4 ) 1115.8 234.0 f 8 7 
. | 5 +13 1, 704.7 359.3 | . 
I 1 30 +1.4 1, 705.0 359 ; ome 
wu) 1s 2. 311.6 SKY. 0 | e 
—~ 1 $2.0 2 309. 6 w3 jf r.4 
. i sO +3.0 » 933.9 23. 4 | ~ 9 j 
. j 1) +8 2, 933. 1 23. 1 j 623. 2 ad 
_ i 30 4.7 3, 573. 2 761.8 ls j 
” 1 40 +4.0 3. 572.0 761.4 jf 61.6 ! 
4) 1 1,003.9 1,054.2 | 
iM) —-- a . 
Hw) +7 4, 806. 1 1,051.2 1, 053. 0 +5 
| a0 +12.0 1, wr 1. O53. 6 | 
Sample No. 2 (3.4104 g)' 
{ 30 0.4 4, 207. ¢ 7H. 4 l » 
ala j "w) ‘j 8 1. 207.6 70.4 { “60. 4 1.0 
i 30 s 1, 309 233. 2 l oa 
100 , 40 1.0 1. 300.8 233.3 | 243. 3 “9 
i 30 +1.3 2,718.1 488.0) \ . 9 
~~ i @ lt 2,718 1 jf oO — 
Sample No. 1 was obtained from Metal Hydrides, Inc 
Sample No. 2 was obtained from an independent source 


undergoing chemical or physical change with time. 
According to table 3, the deviation from the mean of 
two duplicates averaged 0.4 j, with a maximum of 
1.0j. It may be noted that these variations seem 
to be independent of furnace temperature. 

The results of the individual measurements on con- 
tainer with sample are given, in chronological order 
for each of the four samples, in table 4 for sodium 
borohydride and in table 5 for potassium boro- 
hydride. The third column of each table indicates 
the change in the mass of specimen plus container as 
compared with the mass before the heat-content 
measurements were begun. These changes were due 
almost entirely to changes in the mass of the specimen 
itself, as discussed in section 4. The values in the 
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TABLE 5. Individual heat-content measurements on potassium | 6. Thermodynamic Functions, 298.15° to! 


horohudride containe " Jatio 
. 700° K oe 
pound 
Relative t content of EBH 1 : ' mauion 
Fur Sictindaitad ee The observed values of heat content for each COM- peratu 
i ° ~ pele 
nace | Time | Mas measured pound were smoothed by using them to derive the spondi 
tem nm tur excess of heat “er . 7. - MH @ 
pera - uy original sumple Near Observed coefficients of an empirical equation Siving the rela ments 
ee a - pea te egg OF tive heat content as a function of temperature, Fo, natiO¥ 
qt each substance equal weight was given to the ayer! 3). bu 
7 . . . = 7) 
age of the mean values for each furnace temperature oMpo 
Sample No. 3 (3.6008 § given in the sixth column of table 4 or 5.7 For each, Yalu 
: : 
; substance, a form of equation was so chosen that the electe 
m x . . . 
a a LO. | fit. 3 Ny “A $ deviations between calculated and mean Observed jthou 
e- +1 + Rs c values (shown in the last column of table 4 or 5) are aodiun 
—- he 1. 148. INI. 9 ih rs comparable in magnitude to the differences between she we 
. if 30 1, 743.3 275.9 ™ . a. . . \ 
0 | 4 1 742 9 78 2 275.8 +. the individual values from duplicate runs at the more ¢ 
20K) l r - oe a : 360.0 » same temperature.” With the heat content at °C the v: 
wn «=f 2.7 2, 948. 9 413.7 63.8 o in excess of that at 0° C expressed in absolute joules here. « 
é 1 35 28 2. 949.7 3.9 i . . . : lene, 
{30 13 S774 2 ¢ per gram, the equations obtained for the crystalline of pot 
wo =) +0.3 ++ -—~$- a8. S r.2 salts from 0° to 400° C are as follows the sal 
. i 30 +1 1741.7 oO. 4 | = . 
1» 4 - aa a 2 ; { yere © 
1 40 +4 1. 740) 27 . 
m f 3 +.7 2. 340. ¢ hn ae { NaBH: Bureat 
- | i) 2 340 {HN 7 i 2 
r=) } a0 $ 1us. 5 Has. 1 - ' rhe 
, 1 40 1. 197.9 657.9 I, — Hoe = 2 AB84t 8.083 (100 '\¢ forego 
i 30 $845.1 TH I i) 
in) - —=s, ¢ * ’ J o> sym. 7 Wyre ~- 4 : 
\ 3 1H 4 750. 9 133.0 logy [(f+- 273.15) /273.15]. 5) K x 
TAO° K 
Sample No. 2 (3.3794 ¢ 4 
KBH,: ' 
\abB 
10) +0. 2 1, 127 Is! Is 
100 40 ( 1. 127.3 IN}. 4 v2 HH — Hex 4. OS9OSf—2.9030( 107°) 
j ") a $750.4 75a. 4 7350 
in) a - . } . ' - —s) ~ ms - 
7% + os me | — +~2.641(10~°)# — 1,482.7 logy [(f+ 273.15) /273.15 
100 | Hw 9 4 1127 s] } 
1 
* Samples l and 2 were actually specimens o ‘ me s rok ‘ o 
Metal Hydrides, In , . ' 
Saiaiiedl Gnas of weet wieetaten Equations (5) and (6) were then converted to th 
form most commonly used in chemical thermody- 


namics by expressing the heat content relative to 
that at O° K, im calories per mole, with the temper- 
ature in deg K.“ From the resulting equations the | 
cemmon thermodynamic functions (at the standard 
thermodynamic pressure of 1 atm) were derived 
using the following relations 


fourth column have been corrected to the basis of the 
fixed mass of pure compound given in the heading. 
Each value in the fifth column was obtained by sub- 


tracting from the corresponding “‘gross”’ value in the = as ABH, 
fourth column the mean value for the empty con- (H1° — Hy) | T= (HP + ys D+ (H 2s 
tainer (table 3) and then dividing by the mass of 7) (H 
sample. The “calculated” values of heat content 
referred to in the last column of tables 4 and 5 are C“—d(H He) /dT. (S 
the final smoothed values derived as deseribed in 
section 6. *7 ( 
P , S > Hy "\d' (9 
For some temperatures, two mean values of net “0s / I, 
. . . . . ”s 
heat content are listed in the sixth column of table 
5. If both these values were exactly correct, they (F H,)/T=S (H H°)/T. (10 
would naturally be the same. Their small disagree- 
' ameenn item | *otaaai barn. , 
+s Ae “ip a ny eee 1. r, ola tum sg where C, is heat capacity at constant pressure, S is 
. _—~ " or “G09 ‘ “ " “— ae - : A 
he rice aa NO ne loug 1 N mat specimen a entropy, F is Gibbs free energy, and the subseript 
t 1 Same s.vock sample as sampie NO rut measures ‘“()”’ refers to 0° K Equation (10 permits the cal- } 
2 months later, gave mean values of heat content 
lower than those of sample No. 1 by 0.5 percent at anes re 
100° and 0.1 percent at 400°C. The consistent de- a St, Shp Save Saath Tees Geet Se Tee 5 Oe ae Se SRS 
. F : . - - were each given half weight 
crease with time of values for potassium borohvdride * The logarithmic terms and those in ¢ and f# were introduced "= they 
. . . I e been found to account closely for the relative heat content at these ten 
indicates that, to at least the extent of these rather nesatenes of » number of oryetallin slits aoveld of parent transitions 3 
iscre 1e > eorrect} e for chs j In the case of potassium borohydride, an equation containing also a term in @ 
small dist r pare i 5, the ( Orres tion for ‘ hanging was hy whe ‘ion ver, the resulting equat on gave values almost identica The 
degree of purity (incorporated in the values of table | with those of eq @ he 
e\ Taking one defined caloric 4.1840 abs j,0° C =273.15° K,and the molecular rature 
©) are not completely adequate weights of NaBH, and K BH, as 37.843 and 53.952, respectively [4 sur 
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S° to ulation of the free energy of formation of the com- 
sound at Various temperatures when its heat of for- 
mation and its absolute entropy at some one tem- 

‘h com. perature are both Known provided that the corre- 

functions of the free ele- 

The heat of for- 

r, mation measured 

C aver. (5) but no corresponding value for the potassium 

PTature ompound was found im the literature 


IVe the 
le rela. 


ponding Thernloc Ny Ter 
ments also have been evaluated 


" borohvdande has been 


ot sodium 


reach, Values for the constants in eq (7) and (9) were 
lat the: selected from the work of other investigators. 
served Although the low-temperature heat capacity of 
5) are eodium borohydride had been measured earlier IG}, 
‘twee, the work of Johnston and Hallett \7 | appears to be 


at the’ more comprehensive and more accurate, leading to 
tt? C she values of (E15 us Hy) and Sys adopted 
joules here, 3,889.5 and 24.26 respectively.’ In the case 
talline of potassium borohydride, the values adopted for 
same constants, 3,948.8 and 25.48 respectively, 


he 
vere obtained in preliminary measurements at the 
Bureau [8]. 

' The thermodynamic functions derived in_ the 


oregoing manner are as follows (in cal mole~' deg 


» kat T°K and for the temperature range 298° to 
TH0°K _ 
\aBH,: 
H°—1T) T=$18.0607 + 7.3108007%) 7 
3.15 
1 .202.94(1/ T)log,, 7 +831.4/ 7, (11 
i) 
( 18.061 + 0.014622 T7—522/T (12 
o the 
Lodly- S°=41.586 logy T+ 0.014622 T+ 522/T 
ve ne S4 TEA (13 
nper- 
S the _ - ome on a 
Ae } I] / }] OSD log / 1. 2PO2.9(] / ) logy, 7 
rived tL 7 311(1073) T—309/T — 102.815. (14 


ABH,: 


7) H°—I1,) T—80.8100—0.06534038 T+ 3.4055 (107°) T? 
S IQIIMISCL 7) logy 7+ 32,070.2/T, (15 
( SO. 8$10—0. 1230681 7+ 1.0217(1074) 7? 
q 
8 .303/T. (16 
(10 ve a 
S 186.072 log, 7 —0.130681 7+ 5.108010 T° 
S is 8,303 /7—428.371. (lv 
‘ript 


H,)/T T+19,119.21/7) logy T 


186.072 log 


en 0.065340 7' + 1.703(1075) T?—23.767/7 
thes 509.18] 1S) 
te 
ns [3 - 
mer 
ntica I \ ' of “4 
Though th } ( ‘ red il 273° kK . 
PCuLaT iture ler h , 4 } ; } ra if 
surement and tl for ‘ I lerive I 1) I 


Values of the thermodynamic functions calculated 
for even temperatures from eq (11) to (18) are given 


in tables 6 and 7. 


Thermodynamic 


TABLE 6. 


functions of sodium borohydride * 


remper H iH ; S (* H ) 
ure 1 ‘ r 
l¢ al dea l dea cal deq 
A iole nole nole mole 
20s. 15 13. 04 20.67 24. 26 11.21 
sm) 13. O43 “71 24,39 11.29 
$20) 13. SSI 21. 11 25. 74 12.15 
$40) 14. O35 21.50 27. 03 12.00 
st) 14. 460 1.87 28. 27 13. 81 
480) 14. 860 22.24 4H} 14.00 
4M) 15, 238 > i) 30. 61 15. 37 
12 15. 507 ») $1. 72 16. 12 
0 15. G4 23. 31 32. 80 16. 86 
seu) 16, 26S 23. ¢ 33. 84 17. 57 
iN) 16, S82 3. 09 34.8 18. 27 
TLD 1H, SAS 24. 33 35, 34 1S. 05 
*Y 17. 178 4. 46. 80 14, 62 
0 17. 461 24.09 37. 74 20. 28 
wil) 17. 736 25. 32 38. H5 20. 92 
am) 1S. 008 25. 64 $0. 55 21. 54 
wn) 1S, 203 » ” 40. 42 22.16 
Ht Is. Sle 4}, 2S 41.28 22. 76 
“40 Is. 74 26. 00 42.12 23. 35 
twin) 14. O07 20. W2 12.04 23. US 
Hs 14. 244 27. 24 43.75 4.51 
TOO) 19. 477 27 4.54 25. OF 


Molecular weight = 37.843, 


raken from reference [7 
TABLE 7 Thermodynamic functions of potassium borohydride * 
remy H H f S ( k H ) 
ul 1 ‘ T 
cal dea cal dea cal deq cal deg 
hk nole mole mole mole 
Us. | 14. 244 13. 0S 25. 48 12. 23 
in) 13. 400 3. 12 25. 62 12. 32 
320 13. 431 3. 51 27.13 13. 20 
340 14. 508 23.74 Js 6 14. 06 
se) 15. 023 23. U4 24. G2 14.9) 
380) 15. 405 24.0, 31. 22 15. 73 
1K | as 24.13 42. 46 16. 53 
120) 16. 317 24. 18 $4. 4 17. 32 
40) 16. 67 4 2D 34. 76 18. 09 
Wi) 17. (4M 24.27 * 35. 84 18. 84 
int) 17. 308 24.33 36. 87 19. 57 
uM) 17. 500 24. 41 37 87 20). 28 
*Y 17. 854 24, 52 $8. 85 20.9% 
40) Is. 108 24.06 40, 75 21.05 
wr) Is. 340 24,84 1). 66 22. 32 
NS IS. 50S 25.07 41.53 2. 
ww) Is. TSU » $4 a 23. 00 
azn 10. ON) 25. ¢ 43. 22 24. 22 
Ho 14, 20 HO $4.04 24, 82 
TLD 10. 433 2H. 44 14.85 25. 42 
HN 14. 4s Hh US 15. 05 26. 00 
TH) i St 27. 5A 1). 44 AT 
Molecular weight = 53.952 
laken from reference [8 


7. Discussion of Results 


Although the thermal measurements reported in 
this paper are of relative heat content, it is more 


informative to examine 
capacity. 

in figure 1. 
calculated from eq (12) 


the resulting values of heat 


These are shown for the two compounds 
The curves represent the smooth values 


and (16). The points may 


be considered the unsmoothed values, as each one 
represents the difference between the mean values of 
heat content for two adjacent furnace temperatures 
as given in the sixth column of table 4 or 5.” 


The points so mputed 


, but the largest such corres 


wert 
curves { 


capacity 


wrected for the of the smooth 


curvature 


on amounted to only 0.2 percent of the heat 














Hleat 


borohua le, as determined by the drop ethod 


Piagure | 


capacity oft odin horot / ‘ and potas- 


, NaBHy, mean observed: ©, KBH a 1. Neill 
KBH 


It will be noted that in the temperature range in- 
vestigated, the heat capacity of sodium borohvdride 
Increases with temperature in the usual manner ob- 
served for many crystalline On the other 
hand, the heat-capacitv-temperature curve of po- 
tassium borohydride shows a distinet point of in- 
flection, which the existence of a 
second-order transition of some sort in the 


solids 


stiggvests broad 
nelghbor- 
hood of room temperature 

The authors’ smoothed values of heat capacity may 
be compared with the results of other investigators 
The present measurements on sodium borohydride 
overlap those of Johnston and Hallett [7] from 
301° K The agreement in mean heat Capacity 
over this temperature range is almost perfect, but 
their value at 298° Kk 0.25-percent highet 
The less precise values for sodium borohvaride of 
Boodman, Stegeman, and Mason [6] are somewhat 
lower in this temperature range than those of John- 
ston and Hallett Davis, 


)y-- 
mid 


Is about 


Mason, and Stegeman [5] 
obtained a value of 0.55 0.01 eal e' dee Co! for 
the mean heat capacity between 273° and 298° K; 


thes is 2-percent higher than the authors’ correspond- 


ing value. In the case of potassium borohydride, 
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Furukawa et al. [8], who use an adiabatic 
to make the preliminary 
temperatures to 370° K referred to In section P 
obtained somewhat higher values of heat canal 
than the authors. In the range from 298° to 379 k 
the discrepancy averaged 0.5 percent and reached 4 
maximum of 0.5 pereent.' 

The authors believe their 


that corrected 


» 4 ang 
smoothed values of the heat capacity of these two 
compounds are reliable to within | percent. Th, 


thermal measurements were made using a method 
and apparatus that appears to give consistently muc 
greater accuracy than this when the substance 
measured is pure, stable, and physically reprodueib 
It is therefore likely that in the present case thy 
uncertainty in the exact chemical compositions of 
the samples contributed the major uncertainty in th 
heat capacity, and, as indicated in section 4, by 
approximately the same percentage for each of thy 
two compounds over the whole range from 0 to 
100° C. If this is true, approximately the sam 


percentage error for each compound resulted in the: 


parts of the entropy and = heat-content funetio 
which are in excess of the respective values at 298° K 


This work was supported by the Bureau of Aero. 
nauties, UL S. Department of the Navy. Th 
authors thank R \ Paulson, of the National 
Bureau of Standards, for the chemical analvses of th 
samples following the heat-content measurements 
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Infrared Spectra of Crystalline Polyphenyls 


James E. Stewart ' and Max Hellmann 


af ‘ i . ~ . . 
Land Infrared absorption spectra of sixteen selected polyphenyls ranging from biphenyl to 
Se two m-octaphenyl are presented for the 5- to 38-micron region. The observed vibration bands 


Th are discussed in terms of a pseudo-symmetry approximation in which the individual rings 
nethog ure assumed to vibrate separately. Indications of interactions between rings are seen in 
i] 


the C-H out-of-plane bending vibrations of ortho- and metasubstituted rings. Evidence is 

 Muc} viven for the existence of two rotational isomers of m-quatraphenyl with a free energy differ- 

Stance ence of about 250 calories per mole. 

ucibh; 

~ 1. Introduction Randle and Whiffen [5] (referred to hereinafter as 
“mag RW) have analyzed infrared and Raman data from 
ri re Studies of the vibrational spectra of large molecules the literature for a large number of substituted single 
of the have been mostly restricted to empirical correlations | Te aromatic compounds —_—s the pseudosymmetry 
0°, of bands with the presence of specific functional | 2Pproximation and have arrived at some very satis- 
sam. groups. Recent publications, however, give evi- factory assignments. They have also calculated 


some out-of-plane vibration frequencies for aromatic 


in thee dence Of progress in the spectral analysis of large ! 
molecules using the force constants of benzene. 


netio molecules such as polvethvlene {1 substituted *. , Consvan' : 
yok | polyethylenes [2], and infinite chains of HgNH,) More recently, Shimanouchi, Kakiuti, and Gamo {6} 
™ | ta. (3). These substances have the property of have also reported calculations made in this way on 
possessing translational svmmetry and symmetric single ring compounds. cae 
repeating units, thus reducing the secular equations The polyphenyls were studied in the solid phase, 
Aero-| t9 a good approximation, from order infinity to | rather than in solution, because it was desired to 


The) some manageable finite order. The homologous observe, if possible, the contributions of the indi- 


tional) gries of polyphenvls, C,H,-(CyH,),-C,H;, have sym- | Vtdual units of a polyphenyl chain to the charac- 
of the, metric repeating units, but exeept for tne para- teristic absorption bands. In liquids, the variations 
nts ylyphenvls, they have no translational symmetry. 12 environment tend to broaden absorption bands 
| However, empirical speetra-structure correlations | Until the finer structure is lost. On the other hand, 
d analogies with benzene and the substituted | the analysis of crystal phase spectra is often com- 
enzenes [4] have provided information on the nature | Plicated by the occurrence of lattice vibrations and 
J. A fthe bands found in the spectra of aromatic com- the strengthening of combination Vibrations — and 
ounds and, furthermore, have shown that many of | Splitting of bands by intermolecular vibrational 
etical these vibrations are virtually independent of the | Interactions. 
“a nature of the substituent but do depend only on its 
U.S} position 2. Experimental Details 
19 In view of these considerations it seems reasonable ; 
Soa; | attempt to treat the vibrations of polvphenvls in. Infrared spectra of 16 polyphenyls are shown in 
4 terms of the svmmetry of the individual substituted | figures | through 8. The frequencies of some of the 
phenyl units. This is equivalent to replacing the bands and their assignments are given In table 1. 
substituents on a viven aromatic rine with identical A Perkin-Elmer model 21 spectrophotometer was 
‘ONI point masses. A) pseudosvmmetry can then be used for the measurements. The instrument was 


equipped with a sodium chloride prism for use in 


Oy assigned to the resulting structure and the vibrations P . . 
the 5- to 15-u region and a cesium bromide prism 


f the ring classified under the appropriate symmetry 


species in the usual way. This paper is a report on 
in Investigation into the feasibility of assigning in 
this Way the infrared absorption bands of a series of 
polyphenvls containing from 2 to 8 mono- and disub- 
ituted aromatic rings. These compounds have 
the following advantages: The ring-to-substituent 
lorce constants are probably nearly the same in all 
the compounds; the the substituents are 
juite large (at least 77 atomic mass units 
imber of phenyl units between a given ring and 
he last ring in the molecule is not likely to greatly 
fect its vibrations; and, confusion from overlapping 
minds of nonphenyl substituents is avoided. 


masses of 


so the 


for the 15- to 38-u region. Reststrahlen filter plates 
of lithium fluoride and calcium fluoride were used 
to control stray radiation in the 15- to 30-u and 
30- to 38-4 regions, respectively. Stray radiation 
did not exceed 2 percent. 

Samples were prepared in the form of alkali halide 
pellets using the die and procedure reported earlier 
7]. Potassium iodide was found to be most satis- 
factory in reducing scattered energy losses, but even 
so, scatter losses at shorter wavelengths were some- 
times excessive. No attempt was made to com- 
pensate for scatter or reflection losses. A satisfactory 
spectrum of biphenyl in pellet form was not obtained 
in the 5- to 15-u region. The spectrum shown is of 
a specimen of biphenyl melted and recrystallized 
between sodium chloride windows separated with a 
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M, 0, m, p refer to mono-, ortho-, meta-, and parasubstitute: 
H, and | C refer to vibrations which are primarily in-plane ¢ 
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were identical, as were spectra of crvstalline films 
deposited on rocksalt windows by cooling the melted 
material. After making due allowances for the effect 
of orientation [8] on the band intensities, it 
that the spectra of the ervstalline films most nearly 
resemble the pellet spectrum of m-quatraphenylI-Il1. 
A specimen of m-quatraphenyl was melted in vacuum 
and recrystallized. The spectrum of the resulting 
material was that of m-quatraphenyvl-Il. When the 
same substance was dissolved in benzene and re- 
crystallized by slow evaporation, a spectrum of 
m-quatraphenyl-L was obtained. The spectrum of 
m-quatraphenyvl-IT shown in figures 3 and 4 evidently 
contains some m-quatraphenyl-l. In some prepara- 
tions the bands at 784 and S80 em-—', for example, 
were greatly weakened The possibility that rota- 
tional isomerism might be responsible for the two 
forms was investigated by studying the effect of 
temperature on the spectrum of a melted film. A 
cell was constructed for use at temperatures as high 
250° CC. Spectra were observed between 2 and 
l5u with a Beekman IR-4 spectrophotometer equip- 
ped with sodium chloride prisms. The results 
this study are discussed in section 4 

The polyphenyvls measured were obtained from 
different Biphenyl, and o-, m-, and p-ter- 
phenvl were Eastman chemicals which had been 
purified by reerystallization. The higher m-polv- 
phenyls, viz, quatraphenvl-Il, quinquaphenyl, hexa- 
Phenvl, and octaphenvl were prepared in_ this 
laboratory by R. L. Alexander, Jr. [9]. The samples 
of o- and p-quatraphenvl were purchased from the 
Organic Specialties Co., and the sample marked m- 
quatraphenyvl-IL was bought from the Mathieson Co 
The remaining five compounds quatra- and quin- 
quaphenvls) were obtained from the laboratory of 
Prof. G. F. Woods at the University of Marvland 
Their preparation and properties have been reported 
in the literature [10.11.12] 

The spectrum of p-quatraphenyl contains weak 
bands at 740 and 714 em in the 
spectrum given by Cannon Sutherland [13] 
although the spectra seem. to agree in other respects 
These 
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3. Discussion 


Monosubstituted and ortho- and metadisubsti- 
tuted rings have pseudosvinmetry ¢ para- 
disubstituted rings have pseudosvmmetry Vy. The 
svinmetry species of the vibrations of the substituted 
rings correlate with shown in 
table 2. Vibrations that are forbidden to appear in 
the infrared spectra under the selection rules appro- 
priate for a given symiretry group are shown in 
parentheses. They are actually quite 
appear in the spectra, though probably with reduced 
intensity. The vibrations listed in table 2 and dis- 
cussed in the following paragraphs are arranged on 
the basis of an approximate description, e. ¢g., C —H 
stretebing, solely for convenience. It is recognized 
that more than one symmetry coordinate is involved 
in most normal vibrations. It should also be noted 


OY hile 


those of benzene as 


likels lo 


132 





sy 2 Correlation of summet sper } ; 
PaRLe 2 ! i sy | stituted rings» 
qhut I-}) ( ( l ( C be 1 
ben 
para (1 RB 1 / 1 B R 
het ( B 1 B t B i 
rtho (¢ i i B R R 
mono (¢ B 1 / i BR ‘ 
benzene (D B I / B 
Ou ( Hi 
B s / B B 1 Ry. & 
Bes i B B i,) Ss R 
1) s 1 B B i Ro 
Rhys i B B 1 R 
B I 1 | 
In-plane ¢ H 
I i Bi) S$ B / RB 8 
/ i, 8 BR I Bs 
/ 1s B 1 | R,8 
/ 1 B Bh { Rs 
1 I I I 
c— } 
1 I / B Is / 
i / BR { j 
1 1 I i I / 
1 1 I BR i | 
1 I B / 
( H ! 
i x { Rh B ] ~ / 
1,8 I B 1 1 Bs 
1,5 1 Rt l 1 R, 8 
is 1 B 1 1 I 
l I / f 
that the stretching and bending of  carbon-to- 
substituent) bonds in’ the substituted aromaties 


correlate with stretching and bending of carbon-to- 
hvdrogen bonds in benzene, requiring adjustments 
in the counting of vibrations 


3.1. CC Stretching Vibrations 


Characteristic aromatic rine stretching vibrations 
are found near 1,600 and 1,450 to 1.500 em Vi- 
brations of ortho-, meta-, and monosubstituted 
rings that correlate with the degenerate /., Raman 
line of benzene at 1,595 em infrared active, 
but the corresponding vibrations of VA parasubsti- 
tuted rings are imactive All of the polyphenyls 
have absorptions near 1,600 cm! which are resolved 
into at Following RW. the 
high-frequency component probably contains the 
A, absorptions, and the low-frequency component 
the B, absorptions of the mono-, ortho-, and meta- 
substituted rings. The 1,600 em~' bands are quite 
weak In p-terphenyl and even weaker in p-quatra- 
phenyl, as would be expected with a reduction in 
the mono- to pararing ratio. However, there seems 
to be no noticeable weakening of the 1,600 em 
bands in the other polyphenyls which contain para- 
substituted rings. All of the polyphenyl vibrations 
that are correlated with the /\, band of benzene at 
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* are 


sgh em infrared active, and the spectra of all 
‘the poly phe nyls have a comple x series of bands in 
RW have assigned the bigher fre queney 
porptions near 1,500 cm! to the A, in 
»ono-, ortho-, and metasubstituted rings and the 
g,, mode in pare asubstituted rings; and they assigned 
th lower frequency absorptions near 1,450 em"! 


this re2 710Nn. 
modes 


to 


e B; modes in mono-, ortho-, and metasubstituted 
ngs and the Bs, mode in parasubstituted rings. 


Those polyphenyls that contain meta- and para- 
gibstituted rings have bands near 1,400 em~! which 
re missing or Weak in compounds containing a high 
; proportion of orthosubstituted rings and ip biphenyl. 
These bands are probably C—C stretching modes 
shifted from the 1,450 to 1,500-cm~! region, possibly 
pcause Of the influence of neighboring — rings. 
These bands are discussed further in section 
The C-—-C vibrat'ons correlated with the A, 
Raman line of benzene at 992 em~! are sensitive to 
aibstituents and have not been assigned. Vibra 
jons corresponding the uncertain B,, benzene 
ibration have been placed at 1,324 and 1,315 em 
both weak) in and orthosubstituted ring 
ompounds by RW. Weak absorptions are observed 
year there in the infrared spectra of the polyphenyls. 


3.6. 


to 


Wono 


3.2. In-Plane Ring Bending Vibrations 


The doubly degenerate | 
fbenzene is found at 605 em~' in the Raman spec- 
rm. The and 2B, pair of vibrations in mono-, 
rtho-, and metasubstituted rings which correlate 
vith the 4, benzene mode re be active in the 
nfrared, but the corresponding «1, and B,, modes of 
parasubstituted rings should be inactive under V, 
vlection rules. All of the polyphenyls studied have 
iband at 614+4 em This band, which fairly 
except in the of p-terphenyl and p- 
watraphenyl, is assigned to the B, modes of mono-, 
wiho-, and metasubstituted rings. RW_ place the 
orresponding &,, mode of parasubstituted rings at 
em! in the by strong ab- 
sorption of the monosubstituted rings. <All but 
bree of the polyphenyls absorb weakly at 625+6 
m These bands are assigned to the A mode ol 
monosubstituted rings placed at 620--4 (61845 
min the Raman effect) by RW. In addition, all 
but two of the compounds containing metasubsti- 
tuted rings have 647+11 em~' which we 
ission to the «4, mode of metasubstituted rings. 

The 2, ring deformation of benzene has been eal- 
culated by Lord and Andrews [14] to lie at 1,008 
em. RW assign at 1001+4 and 99945 
‘m to the corresponding A, modes of mono- and 
metasubstituted rings respectively, All of the 
polyphenvls have bands of medium intensity near 
1,000 em which probably represent these modes as 
well as the corresponding B, modes in ortho- and 
by modes in parasubstituted rings. 


3.3. In-Plane C--H Bending Vibrations 
_ The in-plane C 
le in the fairly narrow 
1200 em-!. Of these, 


ring bending vibration 


trong, cause 


H99 — 6H recion obscured 


bands at 


bands 


H bending vibrations of benzene 
region from about 1,030 to 
a band is observed at 1,037 


em~' in the infrared and another at 1,178 em=! in 
the Raman spectrum, while the remaining two 


inactive modes have been calculated [14] to lie at 
1,190 and 1,145 em After removing the carbon- 
to-substituent bond bending modes there remain 
1 A, and 3 B,; 2 A, and 2 B,; 1 A, and 3 B,; and 1 


By, 1 By, 1 By, and 1 A, transitions in the mono-, 
ortho-, meta-, and parasubstituted rings, respec- 
tively. Only the B,, and Byu transitions are infrared 
active in parasubstituted rings under Vj, selection 
rules. 

RW have observed that numerous bands in the 
1,000- to 1,200-cem~' region are relatively invariant 
from molecule to molecule of a given isomeric type. 
They assign 1,029 +5, 1,031 +5, 1,045 +7, and 
1,018 10 em~' to A, mono, A, ortho, B,; meta, and 


by, paravibrations; 1,072 +7, 1,126 +12, 1,081 +10, 
and 1,125 + 10em~'to B, mono, B, ortho, A, meta, and 
By, paravibrations;and 1,156 +5, 1,159 +7,1,165 +6, 
and 1,179 11 em™! to B, mono, A, ortho, B, meta, 
and 1, paravibrations, respectively. These three 
groups of bands are readily associated with bands 
found near 1,040, 1,080, and 1,160 em™~', respectively, 
in all of the polyphenyl spectra. RW further assign 

band at 1,177 +6 em™! to an A, mode of mono- 
substituted benzenes. Bands found near 1,190 
em! in the polyphenyl spectra very likely correspond 
to this 4; monomode, and also to the B, ortho, B, 
meta, and the forbidden B, paravibrations as well. 

The regions of absorption, as given by RW, are 
seen to depend somewhat on isomeric substitution, 
and lead to some complexity in the polyphenyl 
bands. For example, the spectrum of p-terphenyvl 
has a pair of bands with peaks at 1,043 and 1,028 
em-', which can be assigned to the A, mono and 
B., paravibrations, respectively. These separations 
are more than the amount expected from RW’s 
tables, perhaps because of interactions of the Fermi 
tvpe 


3.4. Out-of-Plane C--H Bending Vibrations 


These are the vibrations used most frequently for 
the qualitative analvsis of substituted aromatic 
compounds. RW give the following regions of ab- 
sorption for the strong bands: Mono Bb, 15; 
ortho B, 751 +7: meta B, 782 +9, and B, 8 10; 
and para B,, 817 +13 em These values do not 
hold for strongly electrophylic substituents such as 
fluorine and nitro groups. In the spectra of poly- 
phenvls the metabands are found at 795 +13 and 
SO4 21 em™' and the paraband at 854 +6 em”. 
These values are somewhat higher than those given 
by RW. The orthobands are also found at higher 
frequencies, 776 +9 em™', rendering them distinet, 
or nearly so, from the monobands at 766 +22 em™!. 

The weaker C—-H out-of-plane bending vibrations 


76 


are given by RW as follow: Mono B, 908 +10, By 
and », 962 +6, and 982 +6, A, 837 +10; ortho 
A, 865 14 and 977 +9, B, 934 +11; meta A, 904 

13, B, 964 +10; and para B,, 840 (in the Raman 


effect), Ay or Byg 961 12 em~'. All of the poly- 
phenyls have one or more weak-to-medium bands 
near 850 em 7! which are easily mistaken for the 
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characteristic By, para-absorption. It is likely that 
these bands are contributed by vibrations involving 
the connecting links (see section 3.6). 

A more careful study of the strong characteristic 
C—-H bands reveals some curious progressions. The 
metabands are seen to be resolved into components 
with increasing complexity in compounds containing 
larger numbers of such rings. In those cases where 
only one metasubstituted ring is present the strongest 
peak in the 894 cm™'! region is generally found at a 
somewhat higher frequency than when more than 
one metaring present. In the 795-cm™! region 
(see fig. 9) the strongest peak shifts regularly to 
lower frequencies as the number of metarings in- 
creases, while an increasingly complex set of weaker 
peaks develops on the high-frequency side of the 
major peak. This complex set of weaker peaks 
develops on the high-frequency side of the major 
peak. This complex structure is less resolved in the 
spectra of some of the polyphenyls in the solution 
phase, but it is still recognizable, thus eliminating 
the possibility that lattice vibrations or intermoleec- 
ular vibrational interactions are involved. A plaus- 
ible explanation is that the C—H bending vibrations 
of different metasubstituted rings in the same molecule 
are interacting. This idea can be developed further 
through the use of quantum mechanical perturbation 
theory, which shows that if all interaction terms are 
sel equal to zero except those between nearest 
neighbors, which are assumed to be identical, the 
degeneracy of the vibrations is removed and_ the 
levels are split symmetrically. The inclusion of 
next-nearest neighbor interactions results in a regular 
change in the separation of adjacent components with 
the spacing of the bands becoming compressed at one 
side and expanded at the other. 

The magnitude of the interaction force constant 


Is 
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Wave number. cm™ 


Frat RE Y Bands in SOU-cm™~' re qion for compounds contain- 
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Numbers at right indicate number of metasubstituted rings in the molecul 


m-Heptapheny] was not available for study, but predicted bands are included in 
the figure. The line predicted for m-octapheny]! at 780 cm-~'!, was not observed 
perhaps because of insufficient resolution 





must be about 0.2 of the value of the diagonal te 
in order to explain the magnitude of the Observed 
splitting. The relative intensities of the components | 
of the bands were determined from the elements of 
the dipole moment transition matrices after first 
determining the correct sums of the product functions 
for each value of the perturbation energy. The 


rms 


theory predicts that the component of lowest fre. | 


queney will be the strongest, with the intensities of 
the higher frequency components decreasing more 
or less regularly in agreement with the observed 
spectra, provided that \, the interaction force con- 
stant, is taken to be negative. The physical signifi- 
cance of a negative value is that it is easier to distort 
the hydrogens on a ring out of their normal plane if 
the hydrogens on the neighboring ring are already ou 
of their plane. It is interesting to observe thay 
Kakiuti and Shimanouchi [6] caleulate 
values for the out-of-plane interaction force constants 
of neighboring hydrogens on the same ring, and that 
the magnitude of the interaction constant is about 
0.18 that of the C--H. out-of-plane bending force 
constant. 

Orr [15] has recently explained the broadening of 
the out-of-plane bending bands (960 em=!) of the 
olefinic hydrogens in orthosubstituted transstilbenes 
on the basis of damping accompanying the transfer 
of vibrational energy to a motion involving the sub- 
stituent groups which are quite close to the olefinic 
hvdrogens. This situation differs from that encoun- 
tered in the case of the metasubstituted polyphenyls 
because in the stilbenes energy is transferred bet ween 
vibrational states of considerably different energy 
However, Orr's process is very likely involved in the 
where neighboring rings are substituted in 
different ways. 

The bands characteristic of orthosubstitution near 
776 em”! and the 766-em~! band of monosubstituted 
rings in molecules containing orthosubstituted rings 
are also split into complex structures, but not in so 
systematic a fashion as that displaved by the meta- 
ring compounds. This is probably the result of 
interactions between the out-of-plane C—-H vibrations 
of ortho- and monosubstituted rings, as well as inter- 
actions between orthosubstituted rings. The num- 
ber of peaks resolved in the 760- to 780-em7! band 
For example, o-terphenyl has five 
peaks resolved, but contain only three rings. This 
Is possibly a consequence of the occurrence of more 
than one molecule per unit cell. o0-Terphenyl has 
four molecules in each unit cell [16]. The speetrum 
of o-terphenyl in solution in carbon disulfide does 
not display such pronounced splitting, for only three 
components are resolved. —0-Quatraphenyl in carbon 
disulfide solution, on the other hand, has five com- 
ponents in the 740- to 790-em™! region, just as in the 
solid phase 

No splitting is observed in the out-of-plane C—H 
bending vibrations of the parasubstituted rings, but 
the band is found at a somewhat lower frequency in 


Case 


seems eXCeSSiVe. 


p-quatraphenyl (826 em~'), which contains two 
parasubstituted rings, than in those molecules 
(843 +5 em™') containing only one parasubstituted 
ring. 
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3.5. Out-of-Plane Ring Bending Vibrations 


| The frequency of the B,, out-of-plane ring bending 


made of benzene has been calculated [6] to be 
agg cm RW have placed the corresponding B, 
odes of mono- and metasubstituted rings at 
7 +11 and 690 +15 em”, respectively. These 
hands are quite intense in the infrared spectra and 
ye used With the nearby CH bands in the qualita- 
ive analysis of substituted aromatic compounds. 
RW place the infrared forbidden 2, vibration of 
narasubstituted rings at 699 +6 em~! where weak 
hands are sometimes observed in’ the infrared 
spectra. The corresponding orthoband is of species 


4, and 1s also forbidden in the infrared. 

The 404 cm”! £, benzene vibration is doubly 
degenerate. Its A, component in monosubstituted 
ings, and its 1, component in parasubstituted rings, 
10 and 
according to R W, and sometimes occur 
ys weak absorptions in the infrared. All of the 
plyphenyls absorb weakly near 400 em~', which 
awgests that the «4, components of this vibration in 
ortho- and metasubstituted rings also absorb here. 
The B, component corresponding to the /, benzene 
vibration Is reported by RW to he in the 450- to 
500-em region and to be sensitive to the nature 
of the substituent. The polyphenyls have bands in 
this region which permit some general assignments, 
but until Raman data are available tor these com- 
pounds the situation cannot be described as satisfac- 
wrily as for the Il- to 14-u region. Biphenyl has 
160 em~', and the parasubstituted com- 
ounds have single strong bands at 471 11 «m 
vhich are readily assigned to the B, mono- and B,, 
gravibrations. The orthosubstituted compounds 
have a pattern.of three bands at 5083 + 10, 532 +10, 
and 564 7 em while the metasubstituted com 
pounds also have a pattern of bands in this region 
at454 +20, 514 +12, and sometimes 563 + 10 em 
Other bands usually occur also. In table 1 these 
bands are merely assigned to ortho- or metasubsti- 


44 —+10¢m 


tuted rings 
3.6. Connecting Link Vibrations and Unassigned 
Bands 

tach of the ring-connecting bonds contributes 1 
stretching, 1 torsion, 2 in-plane bending, and 2 
out-of-plane bending internal coordinates to the 
normal vibrations. If these internal coordinates 
acted independently, without simultaneous vibration 
of the other internal coordinates, their vibrational 
frequencies would probably be well below the region 
of observation available in this study. However, 
other internal are almost certainly 
volved and might in fact, predominate. Without 
actually performing a normal coordinate determina- 
tion it is oot possible to guess where these vibrations 
are to be found. 

There js little doubt that the para- and metasubsti- 
tuted ring systems are planar, or nearly so, while 
the orthosubstituted ring systems are nonplanar 
decause of steric hindrance. Hence the ring-connect- 
ing bonds in orthosubstituted systems must have less 
double bond character than in the other structures. 


coordinates 


If the vibrations loosely described as connecting-bond 
stretching actually involve large contributions from 
motions of the ring bonds, it would be reasonable for 
the frequencies to be not much lower than those 
of the “normal” ring vibrations, and they might 
indeed account for the bands observed near 1,400 
cm~' in the meta- and parasubstituted polyphenyls. 
The loss of double bond character in the ring- 
connecting bonds of the orthosubstituted systems 
would then explain the absence of these bands in 
the spectra of the orthosubstituted compounds. 

Five of the polyphenyls have bands of medium 
intensity at 727 +2em™'. A sixth compound has 
a pair of weak bands at 722 and 728 em™ and another 
compound has a weak shoulder at 714 em@'. All 
of these compounds contain at least 1 metasub- 
stituted ring, all but 1 contain 5 or more rings, and 
all but 1 have at least 1 metasubstituted ring which 
is not adjacent to a monosubstituted ring. The 
correlation of these bands with vibrations involving 
connecting links to metasubstituted rings seems 
reasonable. 

There are other bands in the infrared spectra of 
the polyphenyls which are possibly associated with 
the ring-connecting bonds. Among these are the 
bands at 846 +13 c¢m7! mentioned previously, and the 
multitude of absorptions observed in the cesium 
bromide region, particularly those correlated with 
ortho- and with metasubstituted rings. 


4. Rotational Isomerism in m-Quatraphenyl 

The connecting bond between phenyl rings in 
polvphenyl molecules is stabilized by about 7 
keal/mole [17] of resonance energy. This is sufficient 
double-bond character to insure planarity except 
in such molecules as o-terphenyl [16] where the steric 
repulsion is too great. On the other hand, a barrier 
of this magnitude is not too high to permit hindered 
rotation and the rapid establishment of equilibrium 
if more than one stable configuration exists. The 
structure of m-quatraphenyl has not been deter- 
mined, but it is reasonable to assume that the planar 
transform is the most stable one in view of the 
known planarity of biphenyl, p-terphenyl, and 
p-quatraphenyl. The planar cisform is very likely 
also a stable configuration with somewhat greater 
free energy than the transform because of the closer 
proximity of the hydrogen atoms. 

The intensities of all bands of liquid m-quatra- 
phenyl were found to decrease reversibly as the 
temperature was increased, partly as a result of 
thermal expansion, but some bands decreased more 
rapidly than others. Several bands were sufficiently 
well resolved to permit intensity measurements. 
The results, shown in figure 10, are readily inter- 
preted [18] as evidence for the existence of rotational 
isomers having a free energy difference of about 
250 cal/mole. A crude calculation of the additional 
electrostatic repulsion energy in the cisconfiguration 
vields a value of this order of magnitude. Included 
in figure 10 are data for two pairs of bands which 
behave similarly as the temperature is changed. 
The nearly zero slope (£ 20 cal/mole) for these 
pairs is evidence that possible perturbing influ- 
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ee 1 structure observed in the CH out-of-plane deforma- 
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On the Application of Steam-Driven Water Jets for 
Propulsion Purposes’ 


J. M. Burgers' and A. Ghaffari ' 


Calculations have been made concerning the momentum that can be given to a jet of 
water by mixing it with a jet of high-speed steam. Apart from an application of the equa- 
tions of momentum and enthalpy, this raises questions concerning the speed of condensation 
ind the acceleration of the water 


1. Introduction 


Basic equations for the chenge of temperature 
and of momentum to be expected when jets of steam 
and water are mixed have been viven by Chiarulli 
and Dressler | 1] 

The present xuthors extended these celeulations 
‘oa ease Where momentum must be communicated 
toa jet of water by mixing it with a jet of high- 
speed steam in such a way that complete condensa- 
tion of the steam will take place. It was thought 
ler constant Pressure would have 

rtain advantages. On the hand there is a 
much smaller loss of useful momentum, while on the 
other hand the thermodynamic aspect of the con- 

nsation becomes simpler It was thought 
that there would be less chance for the auppearaLce 
of shocks. = In under 
constant pressure, the chamber in which the condersa- 


that conde nseation Mie 
cone 


vlso 


order to obtain condensation 

tion takes place must have a sectional area decreasing 
the direction of the flow aus Is applied in the 
yector for steam locomotives necordimg to some 
w that must be determined 

lt was thought useful to have a concrete example 
mind with definite numerical values. It was 

sumed therefore that the pressure in the mixing 

hamber should be kept equal to OLS atm 

The water is supposed to enter tito the mixipg 

hamber in the axial with a velocity 
lOm sec and a temperature of 7’=20° C. The 


direction, 


athalpy of the water then is //’=20.1 keal kg, 
he density. p’=—998 kg m 
Data heave been taken from W. Koch and I. 


vhmidt, VDI-Wasserdampftafela (Springer-Verlag, 
52): since these data are presented in technical 

we must make use to a certain extent of the 
system in which the kilogram is takea as the unit 
fweight or foree. We shall write kg* im this case, 
nd reserve the notation kg when it refers to mass. 
The steam used ys a driving force is supp sed to 
« obtained from a boiler where the temperature is 
bout 165° C and the pressure 7.5 atm (saturation 
pressure for this temperature After coming from 
the boiler the steam is superheated to 7,=360° C, 


its, 


he pressure remaiming po=7.5 atm. The steam 
the preparation ef this paper wa ipported by the Office of Naval Research, 
Ss. Navy 
Guest Worker, National Bureau of Standards American University 
University of Maryland 
ds, and on leave in University 


National Bureau ef Standar 


from Teher 
gures in brackets t 


indicate tl terature references at the end of this paper 


then has an enthalpy //,=759.4 keal/kg. It will 
be assumed that the expansion of the steam from 
7.5 atm 0.8 atm without frietion 
and without change of entropy. If losses and 
entropy change cannot be avoided, we suppose 
that they are compensated for by a suitable in- 
crease Of 7) and po, so that we shall still have the 
sume data for the expanded state as are given 
below. After expansion the steam will then have 
the following properties: 


to occurs losses 


pressure p’’—0.8 atm=8,000 kg* m-* 
temperature 7°’ =99.0° C 
density p’ =0.462 kg m 
enthelpy /7’’=639.4 kerl ke 
velocity u’’=1,008 m see. 
It is assumed that also the velocity w’’ of the steam 
is directed axially. 

The entrance area of the mixing chamber will be 
denoted by A; of this the part 1 s is taken in by the 
water jet, and the part A (1—s) by the steam jet. 

We write 1, for the sectional area of the mixing 
chamber at the point where complete condensation 


has been reached; u,, for the speed of the water 
when all is condensed; pm for the density of this 


water, for which we shall use a mean value p,,=985 
kg m°; further 7,, for its temperature, /7/,, for its 
enthalpy. We note that the mechanical equivalent 
of heat in technical units is J=427 kg*m per keal; 
nid that g=acceleration of gravity =9.8 m sec’. 


2. Basic Equations for the Mixing Process 
hold 


for the mixing 


The 


process: 


following equations 


Kquation of contin uity: 


Asp’u’ A (1 S) p’’ iid Aw Pm Um: (1) 
When numbers are inserted, we solve for (A,,/A) u,, 
and obtain: 
(A,,/A) 4, 0.47 +-9.56 s, (2) 
Equation of momentum: 
Asp’u"?+A (1—s) p”’ wu’? =An pm Um’. (3) 
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Inserting numerical values we find: 


1,,/A) u 172 —372 s. (4) 
Fyuation of nthal py: 
Asp’ u’ H’ +u"?/2gJ) + AN —s) pu’ (AH +a’? /2¢J) 
Am Pm Um (Hn t+Um?/2gJ). (5) 
Inserting numerical values we arrive at: 
(Am/A) Un (FL, + Uy?/2geJ) = 258 — 156 s. 6 
From eqs (4) and (2) we can obtain u,,; from (6) and 


(2) we obtain /7/,,+-u,2/2¢J and then /7,, itself; from 
the latter quantity the temperature of the water can 
be found. It is also possible to find (A,,/.A 

The condensed water has the speed u,, and the 
pressure p,—=p’=8,000 kg*/m?. When it is to be 
used for propulsion, it may have to be driven out 
against a higher pressure p,. In such a case it will 
be appropriate after the mixing to make the water 
flow through an expanding duct, in which the velocits 
decreases from the value w,, toa value w,, 
that the pressure will increase from p,, to p,. In 
what follows we have assumed: ?,—2.0) atm 
20,000 kg*/m*. We then find: 


adjusted so 


~ 


U,- =U," —239 (m/sec)? 


Formulas for other values of p, can easily be derived. 

The equations can be worked out for various 
values of s. After some trials it was found that 
s=0.5 gave a convenient case. Larger values of s 
(meaning that more water is admitted tn the mixing 
chamber in comparison with the steam) make the 
exit speed rather low; smaller values of s (meaning 
that more steam is used in comparison with the 
water) make the final temperature higher and may 
lead to incomplete condensation. With s=0.5 the 
following results are obtained: 


Um = 54.5 m/sec 
ue 52.3 m/sec 
HT, = 52.9 keal/ke 
T..=52.9° C 

0.096. 


~ 


corresponds a saturation 
pressure of 1,450 kg*/m*. Since the pressure in the 
mixing chamber is 8,000 kg*/m?, it was thought that 
there would be no danger for incomplete condensa- 
tion provided the time for condensation is not made 
too short (see section 4). 


To the temperature 7, 


3. Momentum and Efficiency of the Issuing 
Jet 


The total mass flow of water in the issuing jet is: 


Um kg/sec= Am Ym Um/g technical units of 


A, p, 


mass per second. 8) 





Here y,, is used to represent the specific weight of t] 

water, which, when expressed in kg*/m’, is given be | = 

the same number 985 as for its density in ko * eli 
When the jet is applied to drive an underwater, ee 


object with a speed of | m sec. the usual arrange 

ment will be that the water is taken out of the oer. 
roundings, where it had the speed \ relative to the 
object, while it is reissued as a jet with the speed x. 
relative to this object. Hence the thrust to be ob. 


at the te 
; involv 
It is cle 
at too 1 
f eonsid 


iL 


tamed is given by: , 
. the velo 
= m sec LO 
thrust = (A,, 7,, u,./g) (u.—1\ * aint 
Ss { Ym , q be ) ke : q egse it W 
the wate 


The useful work done in unit of time is: 
mixing ¢ 


ondens? 


Because 


1 ! 


power= (1, ¥m Ung) (u ke*m see, (1 


presence 
sjent nul 
self, bu 
surface 


ondenss 


It may be possible to choose 
sion has a maximal value 


We then find: 


so that this eXxpres- 
: this occurs when ! Ly 


thrust ; A Ym Um f) u ko* Ya The S| 
a) Whe 
, a | 
eieilan ; An Ym My g) a ke*m see 10a e stead 
> surface \ 
Inserting numerical values with s—0.5, we arrive at: | order 
> the time 
seas approx 15,800 A kg* 9b herat ec 
ete 
360,000 1 ko* the mete! 
rs an 3 a we 
powel ipprox 560, «1 KZ*"Mm sec 10b latter 
| . simple ¢ 
- these formulas Ais expressed Wh wr, When the 1+ <a 
entrance area of the mixing chamber is taken as | J) tre 
100 em? =0.01 m*, the values become: a 
CORSE 
caren —— 13S ke* (Ye) ensatio 
« that 
power approx 3,600 kg*m see 100) | places i 


surface 


In order to obtain a rough idea concerning the dinar 


efficiency of the arrangement, we will consider only | “! !! 

the energy used for the steam production. We have |.” &@ 

assumed that the water needed for this purpose is lO a 

taken from the surroundings, at the temperature of ||. Ko 

20° C, so that the enthalpy is 20.1 keal/kg. It The « 

follows that per kg water to be turned into steam, Mension 
the enthalpy difference //,—/1’ = 759.4 —20.1=739.8 
keal is needed. Since the mass of the steam needed 
per second is given by A (1—s) pu’? =A (1—s) X 468 

ke sec, the energy required per second Is: 

vhere 7 

A (1—s) ©4683 739.3 «427 kg*m/see. af the 

the the 

With 0.5 this becomes: A 0.73 10° ke*m/sec, — density 

and for A=0.01 m? the energy required (at least) is | approx 

730,000 ke*m/see. This would give a maximal | egs up 

efficiency of L410 

3,600/730,000 = 0.0049. Wher 

’(=20 


» Condens 
siturati 


the upy 


The low value of the efficiency should not be used 
as a measure for the merits or demerits of the system, | 
since its merit is in the simplicity of operation. 
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4. Speed of Condensation 


For complete condensation not only a sufficient 
jiference Must be maintained between the pressure 
» the mixing chamber and the saturation pressure 
athe temperature after the condensation, but there 
involved also a question of time. 

itis desirable that the condensation should require 
ot too much time, since otherwise a mixing chamber 
of considerable dimensions will be needed in view of 
» velocity of the water which from u’=10 
msec to u 54.5 msec. In order to have a definite 
vase it Was thought that each element of volume of 
he water should stay not longer than 0.1 see in the 


rises 


mixing chamber. The problem is whether complete 
ondensation will be possi le under this condition. 
we have condensation of steam in the 


Jecause 
cana of water, the question of obtaining a suffi- 
ent number of condensation nuclei does not present 
celf. but the important matter is to have a sufficient 
surface the water in order to make rapid 
ondensation possible. 

The speed of conde msation de ‘pe ‘nds on two factors: 

Whether the difference between the pressure in 
»steam and the saturation pressure atl the water 
surface Will accelerate the steam sufficiently rapidly 
yn order that it can reach the water surface within 
the time limit set above; and (b) whether the heat 
berated by the condensation can flow off towards 
the interior of the water sufficiently rapidly. The 
tter problem depends on heat conduction and a 
simple estimate gave the impression that this was 
the really important question, 

In treating the problem of heat conduction into 
the interior a mathematical difficulty presents itself 

consequence of the circumstance that the con- 
ensation adds water to the already existing surface, 
» that the surface where condensation occurs dis- 
laces itself with time. For the of a plane 
surface a treatment will be viven (a) based on the 
inary solution of the one-dimensional heat equa- 
tion, In Which this displacement is neglected; and 
)}on a solution in which the displacement is taken 
method indicated by 


area of 


cause 


nto account according to a 
lI. Kolodner |2} 


The equation for the heat conduction in one di- 


lension Is, 


or oF , 
‘ { } 

Or Ov" 
vhere y is the distance measured from the surface 
f the water inwards. Here x= pC, where X is 
the thermal conductivity of the water; p is the 


lensity and (is the specific heat. For \ we can take 


ipproximately) 1.4 107° eal see? ema! °C; in 
-Y-s units p ] ana ( 1. so that K has the value 
4x10 Cin sec 

When the initial temperature of the water is 
r’(=20° C) and the surface on which the er 


Wndenses acquires a temperature 7, (=93 the 
‘wturation temperature for the pressure of 0.8 ae 
t the equation is, 


ut appropriate solution of 


—— . ¥ 
T=T’+a ( 1—erf — ) (12) 
Zyxt 
When we neglect the displace ment of the surface, we 
must have T=T, for y which gives, 
a= T,— T’ =73°C. (13) 


The heat flow to the interior then becomes (per em? 
and per SECC ) 


oT ——— K 
( d ( ) C(T,—T’) —° (14) 
/ OY /o , \ mt 
The heat flowing to the interior in the time ¢ is, 
2 
qi pC’ (T,—T’) y xt ~3.08ytecal/em*®. (15) 
\7r 
Taking t=0.1 see we find, 
q,=0.975 ecal/em’. (15a) 


The amount of heat to be taken up by the water 
is found as follows: There enters into the mixing 
chamber per second, A(1—s) p”u” kg steam.‘ 

This steam is condensed into water with a temper- 


ature of at most 53° C (for s=0.5). We take as 
mean temperature of the water: (20° + 53°)/2=36° 5. 
The steam has the enthalpy H”=639.4 keal/kg. 


The heat to be taken up by the water entering into 
the mixing chamber is given by the loss in enthalpy 
plus kinetic energy of the steam on condensation. 
Since the enthalpy plus kinetic energy of a given mass 
of steam remains the same as its value in the boiler, 
59.4 keal/kg, and since the kinetic energy of the 
the rate at which heat is taken 


Sa \ 7 
water is 0.1 ke al ‘kg, 
up by the water is 


AC s)p"u” (759.4—36.5—0.1) 


A™X 0.5 * 0.462 * 1.003 * 722.8 


~ 1.675 A keal/see = 1.675 A cal/sec. 
If we divide this by ¢ we find that per second a sur- 
face of 


210° A em 


1.675 & 10° * A/0.975 & 1. 


must be available. 
There enters the mass of water, 


wu’ ke see=5,000 A kg/sec=5X 10° A q/sec, 


As p’ 


which has the volume 5 10° 


=1,000 kg/m* and w’=10 m/sec). 


A cm? see (p’ 


It thus follows that per cm® of water we must have 
un surface area of 
1.72 10°/5 & 10°=34.4 em?/em’. 

When the water is dispersed into drops of radius R, 
the ratio of area to volume is 3/R. Hence we re- 
quire drops of a radius not larger than 0.104 em, 
that is, drops of a diameter of about 2 mm. 

"4 In this formulas 
sults since A finally 


1 is expressed as m However, this does not effect the re- 


drops out 





When the displace ment of the surface must be taken 





that with spherical drops the heat flow to the in-! The < 


into account, we assume (following Kolodner [2]) | terior is smaller than with a plane surface. Probably {the m 
that the surface at time ¢ finds itself at a position, the effect of the condensation upon the surface will | This mo 
agein give a further reduction by at mest 10 Percent, | js pessl! 
N=2byxt 16) Although the calculations cannot ezsily be made! geting 0 
precise, they make it pessible to obtain an estimate | goment 
We now must have 7=—7, for y A 2b yx for the total surface nrea needed = for the Water analysis 
Instead of (13) this cives entering into the mixing chamber per second, In Order } jpad to | 
= ' to allow the condensation of all the steam entering — peighbor 
T .. 7’ —a(i+erf b 17 in the same period, When it is assumed that each | that 
weter droplet will stay in the mixing chamber joy | gt the st 
A second cordition relating a end b is found by con- | CYS" 0.1 see Fine dispersion of the water is neces. Keepli 
sidering the balance between the heat liberated by | S#°Y. ' such a way that the drop size should be of Mu" ae 
the condensation and the heat flowing to the interior: about 2 mm in diameter. When this cannot be obtain fe 

3s can be written obtained, much more time must be allowed for thy 

condensation 
oT dX k ; ; 
q=—{ ay ) pl. 1, =pl Wa IS) 5. Equations of Motion for the Water Drop. | jiwm wh 
| lets 

“> af pe 603 ow kg—602 eal g Provisional equations have been developed for th 
ith (12) and (16) we obtain motion of the water droplets in the mixing chamber | To fo 
in order to find their speed as a function of thy ' ture of s 
Na exp h2 Lb K distance traversed, This is needed for the determi- | apy sect 
\ wel | — \ f pation of the cross-sectional area which the mixing | while th 


The time factor drops out, and with A= p(x we find, 


chamber must have at various distances from its 
entrance. 
For this purpose it is assumed that the water has 


already been dispersed into fine drops and it is 


( ‘omplet 
has beco 
The e 


a= y whe 19) supposed that condensation takes place on thes Aw 
drops, so that no new drops are formed. 
Substituting this into (17) we arrive at The speed of the water drops will be denoted by The e 
uv’; this speed increases with time, and we can also | pressure 
((T,.—T’ consider uw’ to be a function of the coordinate 
be” (1+erf b 20) measured along the axis of the mixing chamber. Wi Aww, 
Lyn then can write for the drops, dr— wu’ dt, and express 
their acceleration as, du’ dtu’ (du’ dy), The equi 
Now with the data used above The speed of the steam will be denoted by w”. | these eq 
- ; Also this speed can be a function of v. re of t 
CT l’) Ly =0.068; We assume that the mixture of steam and water) From 
thus retains a constant pressure, and that the steam 
b—0.064, and a=68.8° C retains ® constant temperature 7°’. The density of 
the steam, p’’, will then be constant. The tempera- 
From (18) we now obtain ture of the water droplets, on the other hand, With the 
increases In consequence of the condensation of 
J 2oLCy xt~2.90, ¢t cal em steam upon them. The density of the water drops 
p will be considered as constant 
Again taking f‘=0.1 we find We shall write A, for the sectional area of the , Differen 
mixing chamber at the position «. Let the number 
9:=0.916 cal/em?, of water drops per unit volume at 2 be 1; the number  Vo'(u” 
of drops crossing a section of the chamber per second 
a , , is then given by, 
instead of 0.975 cal em? as given by (15a ‘ V 1 nu’ 9) 
It appears from these calculations that the effect , m wh 
of the condensation of water upon the surface makes | As mentioned it is supposed that no new drops ar 
the heat flow to the interior somewhat smaller, | formed: when the state of motion is stationary th 
although the difference remains below 10 percent quantity No must be a constant (independent of 
Since it is necessary to disperse the water into fine | and +). ‘follow 
drops in order to guarantee a sufficient condensation The volume of a drop will be written, w=4rh* 3. | uly ift 


aren, the ecaleulatiors must be repeated for a spher- 
ical surface. This has been done only for the cuse 
where the displacement of the surface is neglected, 
in view of the mathematical difficulties involved in 


The mass of steam condensing on a single drop u 
unit time will be denoted by AZ, which may be @ 
function of s. We then have, 


pdu dt V/. 2) 


a more complete calculation. In this case it appears | 
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in-? The speed of the drops increeses mc Hsequence 
ably | of the momentum gained from the condensing steam. 
will | This momentum Is given by M/ u” per second, lt 
ent.| js ssible thet there is also y frictional force f 
lade | geting on the drops in surplus of the direct gain of 
Nate omentum it mies be. however, that a closer 
ater qgalysis of the phenomenon of condensation will 


b load to the result thet /—0, since the steam in the 


rdey ; 
»ighborhood of the drop will flow towards the drop, 


ring 
re » that there will be no appreciable shearing effect 
not | at the surface of the dre p 
one Keeping both the direct transfer of momentum 
eof Vu" and the frictional force f in the equations, we 
hy btain for a single drop, 
th : 
p’ d(aru’) dt Mu” }. pe 

Op | from which, with the aid of (22 

p wdu' dt=M (u”—wu’ f. 24 
th 
De] To form the conservation equations for the mix- 
the | ture of steam and drops we observe that the part of 
mi- any section si, taken in by water is given by.A,nw, 
ing | while the part taken in by steam is A,(1—nv). 
its | Complete condensation will have occurred when nw 


has become equal to unity. 


has The equation of continuity takes the form, 

Is 

es Apnwp' ul > A, nw) p“u" = constant = Ay. 25 
by | The equation of momentum becomes (since the 
Iso | pressure Is Constant 
WV Anwp'u’? > A, nw) pp?" coustant Ky. 26) 


The equation of enthalpy transport can be added to 
i’. these equations; it can be used to find the tempera- 
re of the drops. We shall omit it here. 


tel j From eqs (25 and (16) we deduce, 
im 
of A.nwp’u’ (ua —au’ Kyu” — kK, 26 
ri 
id, With the aid of (21) this can be written, 
ol 
ps Vp wu” —u’ Kyu’ —h 
he | Differentiation with respect to the time gives, 
0] 
er No’ (u” — ua’ (due dt Vp’ wildu™ ‘dt Np’ widu' dt 
" K,(du" /dt) 
1] 
m which, 
re 
lhe K,—Np'w)(du" dt Ni 28 
t follows thet the speed of the steam u” changes 

3. | wilvif there is an extra frictional force f. 

) 
| 

; 


(25), (26), to 


Equations and (21) can be used 
express A, and n as functions of the other variables, 


pu’ ou") K, (p'u’? pu”) Kk, 9 
S , ” , ” , ’ (2. ) 
pp uu" (u u’) 
Nop ; p ” " ” { u ”” Mu , ) 

i , , ” , , , , , ; Y 

pu p u")K, io & ° p’u”*)K, 
K,u" —K,)p"u" ” 
| ; , ” ” , , fo ’ , , (> ) 

wl(p'u’—p"u")K.—(p'u'*—p”u”)K, 





The central problem is situated in the following 
three equations, 


dw dt VM p 
du’ dt=|M(u" —u’)+f)\/ p’w 
du” dt {/(K,/N— p’w). 


If we attempt to neglect f, the last equation gives 
u”—constant. We then obtein from the other two 
equations, 

dur du’ 


wea 


the integral of which ts, 


wu” u’)=constant. 
the volume of the individual water 
0), 


, 
u 


If we know 
drops at the entrance of the mixing chamber (« 
and if we make a seale for uw’ (we know that 
increases from the value 10 m/see to the value 54.5 
m see with s=0.5) we can find the course of w as a 
function of wu’. 

Equation (30) then gives the corresponding value 
of n. When nw=1, the is completed. 


process 
Finally eq (29) would give the cross-sectional area 


A, as a function of wu’. 
What remains is to find the relation between wu’ 
and vw. This is determined by the equation 


du’ da=(du' dt) /u'’=M(u" —u’)/ p’u’w, 

where du’/dt is obtained from (23) with f=0. In 
order to integrate this equation it is necessary to 
know \/, and this must follow the analysis of the 


condensation process. 
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Mass Spectra and Relative Sensitivities of Some 
Polyphenyls 


Paul Bradt and Fred L. Mohler 


Mass spectra have been measured by evaporating polyphenyls from a tube furnace into 
he ionization chamber of a 60° mass spectrometer operated at constant ion accelerating 


voltage Mass spectra of pure tetraphenyl, 


synthetic mixture of terphenyl, tetraphenyl, 


exaphenyl, and octaphenyl are reported. A 
pentaphenyvl, hexaphenyl, and octaphenyl in 


known proportions was evaporated as the temperature was increased step by step from 125‘ 
to 200° © and time integrals of the five molecule ion currents were recorded over the time 
required to completely evaporate the sample. The relative sensitivity increases rapidly 


with increasing molecular weight. 


1. Introduction 


The Mass Spectrometry Section has published 
several papers on detection and analysis of heavy 
molecules evaporated from a vacuum-tube furnace 
directly into the ionization chamber of a mass 
spectrometer [1], 2, 3, 4].'. Mixtures of heavy mole- 
cules have been studied by recording the spectra as 
the temperature is increased step by step. 
This affords the advantage that the mixture ts 
fractionated in the process of evaporation and the 
complexity of the mass spectra is reduced. However, 
the results are qualitative as the abundance of the 
ions from a given compound will depend both on the 
amount of the compound in the sample and on the 
temperature of evaporation. A quantitative measure 
of the relative abundance of the compounds requires 
a time integral of the ion current over the time re- 
quired to completely evaporate the compound from 
the sample 

At the Bureau there has been interest in- the 
analysis of mixtures of polvphenvls formed when 
biphenyl is irradiated by neutrons and gamma rays 
ina high flux reactor. As a basis for analysis of 
these mixtures, R. L. Alexander of the Engine Fuel 
Section has supplied us with m-terphenyvl, m-tetra- 
phenvl, m-pentaphenyl, m-hexaphenyl, and m-oe- 
taphenvl. He also supplied us With a mixture con- 
taining the five polyphenyv!s in known proportions. 


2. Experimental Procedure 


Mass spectra were measured with a 60° mass 
spectrometer of 6-in. radius of curvature and the 
spectrum was scanned by varving the magnetic 
field with a constant electric field of 1,100 v. The 
repeller electrode of the ionization chamber was a 
line wire gauze that allowed free entry of the vapor 
rom a tube furnace. The svnthetie sample was in 
a copper tube which extended to within about 5 mm 
of the gauze and this was ina 6-mm Pyrex tube which 
fitted into the entrance port of the ionization cham- 
ber. A thermocouple attached to the copper tube 
recorded the temperature of the sample. 

The ionization chamber block attained a tempera- 
ture ofabout 190° C from heat from the cathode and 


— 


Figur } ket ] te thre 1 ture reference it the end of tl pape 


the sample tube reached a temperature of 131° 
without external heating. Heating the sample tube 
in turn raised the temperature of the ionization 
chamber by an amount that was not measured but 
results indicated that there was no condensation of 
material in the ionization chamber with the sample 
at 290°C. In other experiments the sample holder 
was further from the ionization chamber and the 
initial temperature was less than 131° C. 


3. Results 


Table 1 lists some of the larger peaks in the mass 
spectra of m-tetraphenyl, m-hexaphenyl, and@m- 
octaphenyl. A striking feature of the spectra is that 


TABLE 1 Vass spectra of tetraphenyl, heraphenyl, and octa- 
phen yl 


retraphenyl Hexaphenyl Octaphenyl 


( ( ( 
lt los ” 
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) 03 0.7 0.2 
12s 1.2 4 1 
1Av 2.5 11 
- s 13 
( 1.1 1.4 1 
17s os Oo. 6 4 
Isv 1.4 ’ $ 
' 3 i) 
l L165 1.t 2 
vet $ 
Js i 7 
jO8 20.7 
su 0.9 1.1 2 
252 1.2 1.3 i 
a] u 13 
27 $4.0) 2.2 
ISU 6.3 14 1.1 
sr 1+ 
alo 
ver a » 
2u 1.1 0. 22 
i) o9 $ 
$5 2.3 ] 
tsi) if 
H41 2. 3 4 
ow" } 
\7 3 
$2 y 4 
7Y 3 
1S 7 
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all fragment ion peaks are small and get progressively 
smaller as the molecular weight increases.  Frag- 
mentation involving breaking bonds between phenv| 
rings is negligible and most of the fragment ions 
involve breaking one of the rings. Spectra of 
terphenyl and pentaphenyl have been recorded with a 
ISO" mass spectrometer with a heated reservoir and 
inlet manifold [5]. All the spectra show that the 
molecule ion peaks can be used without any corree- 
tion for interference to analyze a mixture of 
polyphenvyls. 

The synthetic mixture contained 23.4-mole percent 
terphenyl, 21.6-percent tetraphenyvl, 17.6-pereent 
pentaphenyl, 22.7-percent hexaphenvl, and 14.7- 
percent octaphenyvl. A drop of molten material esti- 
mated to be about half a cubic millimeter was placed 
in the sample holder. This sample was completely 
evaporated into the ionization chamber over a period 
of 27 hr and the time integrals of the molecule ion 
currents were recorded. Only the molecule ton peaks 
were recorded and only two molecules contributed 
measurably to the spectrum at any one time. Read- 
ings on two peaks could be repeated within about 5 
min. The temperature of the sample was increased 
by small steps to maintain an easily measured rate 
of change. 

Figure | is a plot of the ion currents for the five 
molecule tons (designated @,, Dy. Ds, De. and @) as a 
funetion of time. Upward inflections of the curves 
mark the times when the temperature was increased 
The mass spectrometer filament was turned on at 
2:20 p.mM.and at 3 p.m. the terpehnyl and tetraphenyl 
ion currents began to rapidly rise to sharp maxima. 
Heating of the tube furnace began at 7:30 Pp. mM. and 
over a period of 6 hr the spectrum showed pre- 
dominantly pentaphenyl with temperatures ranging 
from 162° to 174° ©. The hexaphenyl evaporated 
from 187° to 214° © and the octaphenvl was _ pre- 
dominant from 250° to 287°C. There was no visible 
residue left in the sample holdes 
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The rate of evaporation of the terphenyl Was so 
rapid that there is some doubt as to the ace ur | 
the time integral. For the other 
rate of evaporation staved within a measurable! 
range. Table 2 gives the results. Column 3 FIV es 
the time integrals in units of current in arbitrars 
seale divisions and time in minutes. Column 4 giyés 
the composition of the synthetic mixture in mole 
pereent and column 5 is proportional to the time 
integrals of column 3 divided by the mole percent 
abundance of the compound of column 4. 


acy of 
compounds the 


TABI E 2 Re lative sensitivil i ot po sphe nu 


motecul 


{ONS 

\l Rela Mon: 

Moleeuk mie Integral of pet tive Isotope  isotop 

on current cen ! rree ensi 

tivit 

Perpheny 20) , OSxTO 4 03 4 10 1 Ww ) 
retrapheny $00) 14! PL 1. 74 7 LAT 
Pentapheny js oF 17.¢ 2. 58 1.41 2 ON 
Hexapheny 15s USS 2. 4 $ 27 1m) 1.02 
Octapheny “1 7.9 ve 1090 L 7 5 7 

The results show that the sensitivity increases 


rapidly with increasing molecular weight. The ey- 
periment was not designed to give sensitivity in 
absolute units. The long series of observations Was 
made by one observer and the experiment was not 
repeated because of the exceptional phir sical effort 
involved. Figure | indicates that the time schedul 
of increasing the temperature could be improved 


4. Discussion 


There are extensive data m the API Catalog of 
lwss spectral data on heav hydrocarbons obtained 
with IS0° Consolidated mass spectrometers using a 
heated reservoir and inlet svstem held ata tempera- 
ture sufficient to vaporize the sample completely 
Data commonly include sensitivity of the maximum 
peak of the compound relative to the molecule ion ol 
n-hexadecane. These data are not comparable with 
data obtained in this work and indeed it is difficult to 
intercompare the published spectra when there are 
large differences in the molecular weights of the base 
peaks. This is because the ( ‘onsolidated ass spec- 
trometer records the spectrum at a constant magnet 
field anda variable ion accelerating voltage, and there 
is a decrease in sensitivity with increasing moleculat 
weight. In this research the ion accelerating voltage 
was held at a constant value. 

The poly phens ls afford an excellent opportunity to 
study the variation in sensitivity with molecula 


weight, as the molecule ions are by far the largest | 


peaks and other characteristics of the spectra are 
similar in all the molecules. Unfortunately some 
limitations in the experimental data limit the theo- 
retical significance of the On theoretical 
grounds, with chemically similar molecules of this 
tvpe, would expect the total ionization al 


results 


constant pressure to be nearly proportional to the | 


molecular w eight 16,7] 
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The total ion current cannot be accurately evalu- 
ated. The sum of all the measured peaks in tetra- 
phenyl was 1.67 times the molecule ion peak, In 

' hexaphenv! the ratio was 1.76 and in octapheny] 
i328. These are lower limits as many little peaks 
were omitted and the sum of these is not negligible. 
In spectra obtained with a Consolidated mass 
epectrometer and heated inlet svstem the sum of 
the measured peaks in terphenyl, tetraphenyl, and 
hexapheny! was about 2 relative to the molecule ion 

In pentaphenyl the sum was 1.6 but the data 

are incomplete. Part of the difference between 
measurements with the two instruments comes from 
mass discrimination in the Consolidated mass spec- 
trometer It is concluded that the ratio of the 
molecule ion current to the total ion current remains 
the same with the possible exception of 


ped k 


roughly 
octaphenyl. 

There is a large isotope correction in these heavy 
molecules for molecules containing C™ and D. For 
terpheny] the abundance of Ci}ZH}, is 0.816 times the 
monoisotopic abundance of C)Hy, and for 
yhenvl, CyH,,. the factor is 0.582. The reciprocal 
of this fraction is the correction factor to be applied 
to the CH), peak to give the monoisotopic abund- 
The 6th and 7th columns of table 2 give the 
the monoisotopic relative 


octa- 


ance 
correction 
senSitIVItS 

The resolving power of the mass spectrometer is 
about 1/300 and ts insufficient to resolve unit mass 
for the hexaphenv! and octaphenyl ions. For octa- 
phenvl the peaks at 610 and 611 will be completely 
iresolved and there can be a large error in evalua- 
ting the 610. For hexaphenyl the 
error is much less 

It is noted above that theory indicates that the 
effective collision area for ionization of polyphenyls 
will be proportional to the molecular weight. The 
relative sensitivity increases much more rapidly than 
is Thus the ratio of the monolsotopice SCHSITIV ity 
of octaphenyvl to that of terphenvl is 5.76 while the 
2.66. The experimental 


factor and 


abundance of 


molecular welghit ratio ois 


more reliable for tetra- 
phenyl and octaphenyl. Here the sensitivity ratio 
is 3.1 and the molecular weight ratio is 2. The 
experimental uncertainties cannot account for this 
large increase in sensitivity with molecular weight. 
These preliminary results are published because 
the experimental work has been interrupted and 
there is no immediate prospect of getting more 
accurate data. Results are of some interest as thev 
afford the only experimental data on the sensitivity 
of mass spectrometers of this type for heavy mol- 
ecules. The change in sensitivity with molecular 
weight is in marked contrast to data obtained with 
instruments operating at constant magnetic field. 
The API Catalog of mass spectral data gives sen- 
sitivities for some fused ring aromatics [8] and in 
these spectra also the molecule ion is the only large 


ratio is considered to be 


peak. For isomers of C,H, (mass 202) the sen- 
sitivity is 56 times that for the molecule ion of 
hexadecane; for C\,H,. isomers (mass 228) sen- 


sitivity ratios range from 23 to 12; for C,H,» (mass 
252) the ratio is 26, and for C..H,, the ratio is 19. 


The sensitivity per micron shows the same un- 
systematic variation, 
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Infrared Spectra of Thermally Degraded 
Poly(Vinyl Chloride) * 


Robert R. Stromberg, Sidney Straus, and Bernard G. Achhammer 


} 


T 
in vacuum in the range 100 10° CC 
The principal changes occurring in the residue « 
tion of unsaturated structures and a change 
aromatic absorption. The data are used to s 
decomposition for poly vinyl chloride 


il to 


1. Introduction 


kineties of the thermal decomposition of 


The 


poly(vinvl chloride) in a vacuum have been pre- 


viously reported [1].!) Under pyrolytic conditions, 
poly(vinyl chloride) decomposes —predominantely 
’ into hvdrogen chloride and a colored. solid residue. 
The color of the residue has been attributed to a 
polyene structure existing after the dehydrochlorina- 
tion [2]. The role of oxvgen In the decomposition 
mechanism has not vel resolved although 
hydroperoxide mechanisms have been proposed :3]. 
Recently group theory has been applied to the 
nterpretation of the infrared spectra of poly(vinyl 
chloride) and the spectrum has been analvzed in a 
comprehensive manner [4]. The head-to-tail ar- 
angement of the monomer units [5], and the alter- 
ate arrangements of the chlorine atoms in the plane 
of the carbon chain [6], with the conclusion that the 
ervstallogra phic repeating unit of two 
monomer units, were further verified by the work of 
Krimm and Liang |4]. Polarized spectra were used 
vy Krimm and Liang to make band assignments. 
Other, more limited analyses of the infrared spectra 
of poly(vinyl chloride) have also been reported [7, 8]. 
The branching of poly(vinyl chloride) has been 
measured [9] by determining the methyl group con- 
centration of the reduced polymer. The application 
of infrared spectrophotometry to the explanation of 
degradation process occurring in poly(vinyl 
chloride) has, however, not been extensively reported 
the literature. Campbell and Rauscher [3] 
itilized changes in the infrared spectra of poly(vinyl 
chloride) in a study of the base-accelerated degrada- 
tion of this polymer. Others [10, 11] limited their 
studies to changes in only one portion of the spec- 
trum. 
This paper describes the results of an investigation 
Which infrared spectrophotometry was used to 
study changes in chemical structure occurring in the 
solid residue during pyrolysis of poly (vinyl chloride 
in Vacuum. 
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gure ul 


e changes in chemical strueture occuring in poly(vinyl chloride 
wer 


as a result of heating 
e studied using infrared spectrophotometry. 
luring pyrolysis in a vacuum were the forma- 
from an aliphatic spectrum to one showing 
upport a previously proposed mechanism of 


2. Experimental Procedure 


Three poly(vinyl chloride) polymers were studied. 
They were prepared from vinyl chloride monomer as 
follows: 


Polymer Preparation 


PV¢ + Initiator—y radiation from 0.3 Curie 
(‘o—-60 source 

PVC—bp Initiator—0.1 mole percent of benzoy! 
peroxide; 40° © 

PV C—azo Initiator—0.02 mole pereent of 2,2’ 
azobis (isobutyronitrile) ; 40° C 


All the polymer samples were powdered by grinding 
under liquid nitrogen. The white powders were 
sieved through a 325 mesh screen to a particle size 
less than 44 microns in diameter. 

Specimens of these polymers were pyrolyzed in a 
vacuum for 30 min at temperatures up to 400° C 
after outgassing at 110° C for approximately 1 hour. 
At temperatures above 180° C volatile products were 
evolved. The color of the residue changed to light 
brown at the lower pyrolysis temperatures and to 
dark brown at the more elevated temperatures. 

After the exposure to pyrolytic conditions was 
completed, the degraded powder samples were re- 
moved and prepared for infrared analysis. Untreated 
poly(vinyl chloride) is soluble in only a few solvents, 
which include eyclohexanone, tetrahydrofuran, di- 
methyl formamide and to a limited extent, ketones. 
Some of the degraded material had an even more 
limited solubility. Cyclohexanone and tetrahydro- 
furan tend to form hydroperoxides, which decompose 
easily and may initiate decomposition of the polymer. 
If these solvents are retained in the polymer film, 
they also can lead to erroneous conclusions about 
the structure of the polymer. This is especially true 
when hydroperoxide breaks down to give carbonyl 
and hydroxyl groups. Therefore, the solid phase 
pellet method was used to prepare the material for 
infrared analysis. Potassium bromide (Harshaw 
Chemical Co.) was used for the suspending medium 








because it has a refractive index, n, of 1.56 which obtaine 
approximates that of untreated poly(vinyl chloride shown 
(1.55). The polymer and the potassium bromide in A we 
the ratio of 1 to 50 were mixed in a dry box and the in the 


assigne 
yibratt 
formed 


pellets were prepared in the manner described pre- 
viously [12]. Although the ratio by weight of polymer 
to potassium bromide was the same for all samples, 





there was an increase in the molar ratio of degraded 1] of 1 
polymer to potassium bromide as a result of the loss The 
of hydrogen chloride exhibit 
A Perkin-Elmer Model 21 double-beam infrared the ot] 
spectrophotometer with a sodium chloride prism was \773 ¢1 
used for measurements in the 5000 to 670° em the cal 
region. For increased resolution a Perkin-Elmer 1733 ¢ 
Model 12—B instrument with a lithium fluoride prism vibrati 
was used for measurements in the 3000 to ?2HOCO0em interac’ 
region Benzo 
Figure 2 Pr on ¢ nfrere ectrun i Geor—101 PV" L000 1 
3. Infrared Spectra and Discussion ng lithium fluoride prisn from wl 
3.1. Untreated Polymers 3.2. SE 
The infrared spectra of the three poly(vinyl ehlo- ‘Figu 
ride) polymers studied are shown in figure 1. These | pigie 1. Assignments of absorption bands for auntreat.g | meet P 
patterns were obtained by placing a potassium bro- poly(vinyl ch ~ | temper 
mide pellet) containing the untreated poly(vinyl 8 8 fe 
chloride) in the sample beam and a blank potassium o—_* 
bromide pellet in the reference beam. A commercial a pric sanieaa References | of the 
polymer, Geon 1017, was studied with a lithium wane 
fluoride prism in the region of 3000 to 2820 em ww CCL stret graded 
This spectrum is given in figure 2. —— a oe | the po 
The assignments of the principle absorption bands amp | Semeibel cheat oneet treated 
of the basic polymer are given in table 1.) Krimm my CH wag, ¢ phiase approx 
and Liang [4] reported absorptions in the region of — 4 oe vl ge in the 
3000 em~! that varied slightly from those given in ar | Re ek. tes was th 
table 1 and figure 2 It is possible that this small N23 CH stret n phase 4 9) 2 
difference IS a result of the hature and perhaps et ri sanetets ade ms co of ¢ eapenes 
seatter of the pellet specimen. Krimm and Liang - n+ -~- peepee acime ed ealeinacausanidbaines ; ete 
"2 B. F. Goodrich Chemical ( that of 
| 
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pbtained their spectra from a film and the spectrum 
shown 10 figure 2 was obt: aine “d from a pe Hlet. 

A weak absorption band at 1590 em™~! is observed 
» the spectrum of the PVC-y sample. This 
assigned to a carbon—carbon ya. Pi hond strete ‘hing 
vibration [13] The double bond structure 
formed as a result of some dehyvdrochlorination 
i} of the polymer while exposed to the y-radiation. 
obtained from the PVC he sample 


is 


is 


The eurve 
exhibits several absorption bands not observed in 
the other polymers. The bands at 1795 em7' and 


assigned to the 


1773 em are carbonyl vibration of 
the catalyst, benzoyl peroxide [13]. The band at 
1733 em! is probably from the carbonyl stretching 


wr ketone [13], formed from 
with vinyl chloride 
an absorption [13] near 
at 990 em Is probably 


vibration of an aldehyde « 
teraction of the catalyst 
Benzovl peroxide shows 
1000 em and the band 
from unreacted catalyst 

3,2. Spectra of Poly(Vinyl Chloride) After Pyrolysis 
Figure 3 shows the spectra of PVC-bp_ residue 
after pyrolysis for 30 min in a vacuum at several 
temperatures The spectra are progressively altered 
as a result of increasing the pyrolysis temperature 
There is considerable radiation scatter for the pellets 
of the samples pyrolyzed at higher temperatures 
where the refractive index of the more highly 
eraded colored specimens did not match that 
the potassium bromide. The samples were pre- 
treated by heating at 110° C in a vacuum for 
approximately | hour The only changes observed 
in the spectrum of the material after this treatment 
was the disappearance the peroxide bands at 
1795 and em indicating removal of the 
mreacted catalyst by this pretreatment. At tem- 
peratures up through 250 the infrared pattern, 


de- 
of 


of 


1775 


the order of 300° C and higher, there is a marked 
over-all departure from the spectrum of poly(vinyl 
chloride). It was shown in_a previous study [1] 
that at approximately 300° C practically all of the 
chloride atoms have been removed from the polymer 
as HCI without any apparent appreciable break- 
down of the polymer chain. This removal of the 
highly electronegative halogen from the chain could 
cause frequency shifts in some absorption bands as 
well as other spectral changes resulting from the 
new structural configurations. Substantial break- 
down of the polymer chain occurs at 400° C [1]. 

In the region of the carbon—-hydrogen stretching 
absorption, 3000 cm~', a new band is formed at 
3030 em”! that may be associated with the carbon 
hvdrogen stretching vibration from an unsaturated 
carbon carbon linkage [13]. The intensity of this 


band increases with increasing temperatures of 
pyrolysis. Nonconjugated carbon—earbon double 
bonds are suggested by the absorption at 1670 
em-' [13]. The absorption band at 1610 em™' is 


assigned to the carbon—-carbon double bond streteh- 
ing vibration for conjugated bonds, either aromatic 
or aliphatie or both [13]. 

The shoulder at 1590 em™! is assigned to the car- 
bon carbon double bond vibration of an aromatic 
ring [13]. A resonance splitting of the absorption 
bands with the appearance of one band near 1600 
em-' and one near 1500 em! occurs when there are 
three unsaturated carbon-carbon vibrations in a 
ring [13]. The bands at 1610 and 1500 en— are there- 
fore assigned, at least in part, to an aromatic struc- 
ture. The broad absorption bands at 880, 820, 750, 
and 700 em! are present in the spectra of asphalts, 
heavy petroleum oil fractions, coal, ete., [14] and 
indicate the presence of aromatic structures, with 
the latter two indicating carbon—-hydrogen out-of- 
plane deformations of monosubstituted aromatic 


although somewhat altered, is easily identified as | groups [13]. The absence of aromatic structural 
that of poly(vinyl chloride). At temperatures of | bands in the untreated poly(vinyl chloride) and 
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Figure 3. Infrared spectra of PVC 


bp after 30 minutes’ pyrolysis. 
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their appearance and increase in intensity as a 
result of higher degradative temperatures can be 
readily observed. 

The absorption at 1720 em! indicates the pres- 
ence of carbonyl vibrations. This band did not 
appear to any appreciable extent when the pvrolyzed 
material was stored under an inert atmosphere. As 
this spectral change does not occur at room tempera- 
ture in untreated polymer, it appears that the 
double bonds in the degraded polymer are subject 
to attack by atmospheric oxygen, forming such 
oxygenated structures. 

The shifting of the carbon hydrogen vibration 
of CH, [3] from 1430 em™! in the spectrum of un- 
treated polymer to 1455 em™' for the more highly 
decomposed samples may be a result of a loss of 
the chlorine atoms or an increase in the CH, conecen- 
tration [13]. The band at 1380 em~! in the most 
highly degraded material may be assigned to the 
carbon hydrogen deformation frequency for the 
aliphatic C’—CH, bond [13]. 


The removal of the band at 690 em~', assigned to 


the CCl stretching vibration [8, 4], from the 
spectrum of the sample pyrolyzed at 300° C corre- 


lates with mass spectrometer and weight loss mea- 
surements [1] which showed that almost all the 
chlorine atoms were removed from the polymer at 
this temperature. 

The infrared spectra of the PVC~y polymer after 
pyrolysis at several temperatures are shown in figure 
4. The changes in the spectra as a result of increas- 
ing the pyrolysis temperature are similar to those 
observed for the PVC-bp material. The spectrum 
of the specimen degraded at 250° C can be readily 
identified as poly(vinyl chloride Increasing the 
temperature to 275° C gives a spectrum, from 1400 
to 650 em™! that is different from the spectrum of 
poly(vinyl chloride). This is probably caused by 
the almost complete removal of chlorine atoms as 
shown by the absence of a band at 690 em The 





presence of the conjugated aromatic structure 
revealed in the spectra of the i 


. | 
Is 
specimens pyrolyzed 
at 342° and 400° C by the bands at 880, 820, 739 | 


700 em~! and the shoulder at 1490 em There js 
also absorption caused by aromatic structure in the 
spectrum of the specimen pyrolyzed at 275° C | 
The spectra of the p\ roly zed PVC 
shown in figure 5. The same general observations 
made with the PVC bp and PVC—-y polymers also 
apply to this material. However, this sample jg 
observed to have undergone more degradation at . 
given temperature than the other two polvmers did 
For example, no C— Cl stretching band at 690 em 
remains in the specimen pyrolyzed at 253° C, while 
this band is relatively intense for the PVC bp poly. 
mer, pyrolyzed at the same temperature. This is jp 
agreement with the faster rate of decomposition, as 
measured by rate of weight loss for this material as 
previously reported [1]. This faster rate was attrib. 
uted to unreacted catalyst remaining in the samph 
The pyrolysis of poly(vinyl chloride) produces 
volatile products and a solid residue. At 400° ( 
some of these products are volatile at the ltempera- 
ture of pyrolysis but not at room temperature 
This portion is a brown, waxlike material at room 
temperature. The spectrum of this wax is similar 
to that of the more hight, degraded specimens ob- 
served earlier and probably consists of fragments of 
the dehvdrochlorinated polymer. The C—Cl band 
at 690 cm”! is absent, accounting for a flat spectrum | 
observed at the lower wave numbers. <A. broad 
carbonyl band at 1730 em~! is probably the result 
of the oxidation of double bond structures following 
the pyrolysis. The CH vibration [3] in CH, is 
shifted from 1430 em~! in the original material to 
1470 em A band at 1380 em assigned to the 
(-H deformation for the aliphatie C—-CH, bond 
[13], is present in this spectrum as well as in the 
spectrum of the residue of the material p\ roly zed at | 
100° CC. Although there are weak bands at 700, 
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7), and 755 em the relatively intense bands | the oxygenated structures observed. This action 


assigned to aromatic structures in the spectrum of 


‘the residue of the material pyrolvzed at 400° C 


are absent from this spectrum. The differences 
between the spectrum of the residue and the spec- 
trum of the waxlike material are attributed to the 


i diferent molecular weights of these portions. 


4. Conclusions 


Specimens pyrolyzed in a vacuum at temperatures 
up to 400° C yielded some spectra that were decidedly 
diferent from the original material and which varied 
fom each other depending upon the temperature. 
The volatile products obtained from the poly(vinyl 
consisted primarily of HCl and a few mole percent 
of benzene [1]. At 400° C, the volatile products 
produced from the dehydrochlorinated polymer 
consisted, in addition, of larger quantities of benzene 
and various hydrocarbons. The principal altera- 
in the infrared spectra of these degraded 
specimens were loss of chlorine atoms, with complete 
removal at temperatures near 300° C, establishment 
of double bond character, both aliphatic and aro- 
matic, and a possible increase in CH, concentra- 
tion. A variety of carbonyl structures is formed 
when the degraded material is allowed to stand 
exposed to the air for several days at room tem- 
perature. Within the limits of this technique the 
spectral changes are independent of polymerization 
initiators used. 

The mechanism of the decomposition appears to 
be a removal of chlorline as HCl in such a manner 
as to form a conjugated polvene structure, even 
when only small quantities of chlorine have been 
removed. No infrared evidence was found for 
iydroperoxides. It appears reasonable, therefore, 
that oxygen does not play a direct role in the 
dehvdrochlorination process of the systems studied 
here, but adds to the double bond structure, producing 


tions 


‘ 
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may further degrade the polymer by reducing the 
size of the chain. The infrared technique was, 
however, insensitive to measurement of chain 
length, although the apparent rise in methyl group 
concentration is an indication of an increased number 
of chain ends. 

The formation of aromatic structures in_ the 
spectra of the degraded aliphatic poly(vinyl chloride) 
correlates with the presence of benzene in_ the 
volatile products of such decomposition reactions [1]. 
This indicates that an aliphatic conjugated structure 
is present at chain ends. This structure can occur 

a result of dehydrochlorination at the original 
ends, either by initiation or termination, or by sission 
of an internal portion of a dehydrochlorinated chain. 
The dehydrochlorination reaction that has been 
previously proposed [1] would allow a free radical 
at the end of a polyene group in either instance. 
As a result of kinetic motion, various. cyclic con- 
figurations are formed. Because of resonance con- 
siderations, benzene the most stable structure, 
and as a result of this stabilization, the carbon 
carbon bond alpha to the benzene ring is weakened 
and breaks, evolving benzene as a product. The 
more complete removal of chlorine at the higher 
temperatures, together with the increased opportu- 
nity for chain cleavage at these elevated tempera- 
tures, vields a more aromatic structure. Moreover, 
the material thatis volatile at the temperature of 
pyrolysis but not at room temperature apparently 
consists of short aliphatic segments of the residue. 
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